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ABSTRACT 


Two-dimensional  and  three-dimensional  statistical  micromechanical  damage 
models  with  randomly  located  interacting  microcracks  arc  presented  to  investigate  the 
overall  nonlinear  mechanical  responses  of  microcrack- weakened  brittle  solids.  The 
macroscopic  stress-strain  relations  of  elastic  solids  with  interacting  micro-cracks  are 
micromechanically  derived  by  taking  the  ensemble  average  over  all  possible  realizations. 
Several  different  approximate  analytical  solutions  of  a  two-crack  interaction  model  are 
introduced  to  account  for  micro-crack  interaction  among  many  randomly  oriented  amd 
located  microcracks.  The  overall  elastic-damage  compliances  of  microcrack-weakencd 
brittle  solids  are  also  derived  by  further  taking  the  volume  average  of  the  ensemble- 
averaged  stress-strain  relations  over  the  entire  material  mesostmctural  domain.  Some 
special  examples  are  investigated  by  using  the  proposed  methods.  At  variance  with  exist¬ 
ing  phenomenological  damage  models,  the  proposed  framework  does  not  employ  any 
fitted  “material  parameters”.  “Cleavage  1”  microcrack  growth  and  “evolutionary 
damage  models”  within  the  proposed  context  are  also  presented. 
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PARTI 

Effective  Elastic  Moduli  of  Two-Dimensional 
Brittle  Solids  with  Interacting  Microcracks, 

I :  Basic  Formulations 

I.O.  Abstract 

Statistical  micromechanical  formulations  are  presented  to  investigate  effective  clastic  moduli 

% 

of  two-dimensional  brittle  solids  with  interacting  slit  microcracks.  The  macroscopic  stress-strain 
relations  of  elastic  solids  with  interacting  microcracks  are  micromcchanically  derived  by  taking 
the  ensemble  average  over  all  possible  realizations  which  feature  the  same  materia!  microstruc- 
tural  geometry,  characteristics,  and  loading  conditions.  Approximate  analytical  solutions  of  a 
two-microcrack  interaction  problem  are  introduced  to  account  for  microcrack  interaction  among 
many  randomly  oriented  and  located  microcracks.  The  overall  elastic -damage  compliances  of 
microcrack-weakened  brittle  solids  under  uniaxial  and  biaxial  loads  are  also  derived.  Therefore, 
stationary  statistical  micromechanical  formulation  is  completed.  Moreover,  some  special  cases  are 
investigated  by  using  the  proposed  framework.  At  variance  with  existing  phenomenological  con¬ 
tinuum  damage  models,  the  proposed  framework  does  not  employ  any  fitted  “material  parameters”. 
“Cleavage  1”  microcrack  growth  and  “evolutionary  damage  models”  within  the  proposed  context 
will  be  presented  in  Part  II  of  this  •"ries.  It  is  emphasized  that  microstructural  statistical  informa¬ 
tions  are  already  embedded  in  the  proposed  ensemble-averaged  equations  and,  therefore,  no  Monte 
Carlo  simulations  are  needed. 


I.l.  Introduction 


# 


The  nonlinear  mechanical  responses  of  damaged  solids  due  to  the  existence,  growth,  and  nu- 
cleation  of  microdefects  (such  as  microcracks  and  microvoids)  are  of  significant  importance  to 
engineers,  and  have  been  t!  e  subject  of  many  investigations.  See  Krajcinovic  (1989)  for  a  liter¬ 
ature  review  on  damage  mechanics.  For  brittle  materials  (e.g.,  concrete,  rocks  and  ceramics),  in 
particular,  microcracks  often  control  overall  deformation  and  failure  mechanisms.  To  date,  the 
only  exact  results  derived  for  microcrack-weakened  brittle  solids  are  for  dilute  microcrack  con¬ 
centrations,  where  microcrack  interactions  are  entirely  neglected  (i.e.,  “Taylor’s  models”).  On 
the  other  hand,  “effective  medium  methods”  were  proposed  in  the  literature  to  account  for  inter¬ 
action  effects  of  microcracks.  For  example,  the  “self-consistent  method”  [Hill  (1965)]  was  first 
applied  to  microcrack-weakened  solids  by  Budiansky  and  O’Connell  (1976)  with  special  attention 
directed  to  perfectly  randomly  distributed  and  weakly  interacting  microcracks.  The  self-consistent 
method  was  further  developed  by  Horii  and  Nemat-Nasser  (1983)  to  take  into  account  the  effects 
of  closed  microcracks  undergoing  frictional  sliding.  Christensen  and  Lo  (1979)  proposed  a  three- 
phase  “generalized  self-consistent  model”.  The  “differential  scheme”  was  investigated  by  Roscoc 
(1952,  1973),  McLaughlin  (1977),  and  Hashin  (1988).  Funher,  the  “Mon-Tanaka  method”  was 
developed  by  Mori  and  Tanaka  (1973),  Benveniste  (1986),  and  Zhao,  Tandon  and  Weng  (1989). 
Some  comparisons  and  assessments  for  the  self-consistent  method,  the  generalized  self-consistent 
method,  the  Mori-Tanaka  method,  and/or  the  differential  scheme  were  also  presented  by  Horii  and 
Sahasakmontri  (1990),  Laws  and  Dvorak  (1987),  Nemat-Nasser  and  Hori  (1990),  and  Christensen 
(1990,  for  pure  shear  load  only).  It  is  noted  that  the  foregoing  effective  medium  methods  are  only 
valid  for  low  microcrack  concentrations  since  they  do  not  depend  on  locations  of  microcracks. 

Ah  of  the  aforementioned  work  can  be  categorized  as  "stationary"  micromechanica!  models 
since  all  microcracks  are  assumed  to  be  stationary,  i.e.,  no  microcracks  are  allowed  to  grow  or 
nucleated  during  loading  histories.  For  a  constitutive  theory  to  possess  predictive  capability,  how¬ 
ever,  an  "evolutionary"  micromechanical  damage  model  is  warranted  to  account  for  “cleavage  1” 
(pre-existing)  microcrack  growth  and/or  “cleavage  2”  (new)  microcrack  nucleation.  In  the  current 
literature,  there  ar»  indeed  a  number  of  micromcchanical  “evolutionary”  damage  models  available. 
See,  e.g.,  Krajcinovic  and  Fanclla  (1986),  Fanella  and  Krajcinovic  (1988),  and  Ju  {1991b)  by  using 
the  “Taylor’s  model”;  as  well  as  Sumarac  and  Krajcinovic  (1987, 1989),  Krajcinovic  and  Sumarac 
(1989),  Ju  (1991a),  Ju  and  Lee  (1991),  and  Lee  and  Ju  (1991)  by  using  the  self-consistent  method. 
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hen  microcrack  concentrations  arc  higher  and  microcrack  spacings  are  closer,  strong  micro¬ 
crack  interactions  occur  and  effective  medium  theories  arc  no  longer  appropriate.  Emanating  from 
this  viewpoint,  excellent  strong  microcrack  interaction  models  were  proposed  by  Gross  (1982), 
Horii  and  Nemat-Nasser  (1985),  Hori  and  Nemat-Nasser  (1987)  for  two-dimensional  determin¬ 
istic  microcracks,  and  by  Kachanov  and  Montagui  (1986),  Kachanov  (1987),  Chudnovsky  ci  al. 
(1987a,b),  and  Kachanov  and  Laures  (1989)  for  two-  and  three-dimensional  deterministic  arbitrary 
micTocrack  arrays.  The  foregoing  work,  nevertheless,  are  only  stationary  strong  microcrack  inter¬ 
action  models.  Moreover,  the  work  due  to  Kachanov  (1987),  Kachanov  and  Laures  (1989)  rely  on 
Monte  Carlo  simulations  of  deterministic  microcrack  arrays,  and  depend  on  heavy  numerical  com- 
putations  of  stress  “transmission  factors”.  Therefore,  it  is  desirable  to  develop  simple  statistical 
micromechanical  theories  to  account  for  interactions  among  many  randomly  located  and  oriented 
microcracks.  Furthermore,  “cleavage  1”  microcrack  growth  models  are  needed  under  the  condition 
of  microcrack  interaction. 

The  purpose  of  the  present  work  (Part  I  and  Part  II)  is  to  establish  a  statistical  micromechanical 
framework  for  deriving  “evolutionary”  damage  models  and  corresponding  constitutive  equations 
for  brittle  solids  containing  many  interacting,  randomly  distributed  slit  microcracks.  The  proposed 
statistical  framework  considers  the  probability  and  conditional  probability  density  functions  of  mi¬ 
crocrack  locations,  orientations,  lengths,  and  relative  configurations.  See  Batchelor  (1970),  Batch¬ 
elor  and  Green  (1972),  Hinch  (1977),  Willis  and  Acton  (1976),  Chen  and  Acrivos  (1978a, b)  for 
references.  In  addition,  the  ensemble-volume  averages  of  stresses,  strains  and  compliances  are  sys¬ 
tematically  constructed  based  on  analytical  micromechanics  solutions  and  probability  functions. 
It  is  emphasized  that  the  proposed  method  is  very  different  from  that  proposed  by  Hudson  (1980, 
1981 , 1986).  Though  using  the  ensemble  average  approach,  Hudson’s  method  is  based  on  a  second- 
order  stiffness  theory  and  therefore  leads  to  irrational  behavior  for  solids  with  moderate  or  high 
microcrack  concentrations  [see  Sayers  and  Kachanov  (1991)].  The  proposed  approach  is  free  from 
this  anomaly. 

A  brief  outline  of  this  work  is  as  follows.  In  Section  2,  an  ensemble -average  approach  to  de¬ 
rive  damaged  stress-strain  relations  is  introduced.  Approximate  closed-form  analytical  solutions 
are  subsequently  presented  for  the  interaction  problem  of  two  arbitrarily  located  and  oriented  (but 
non-intersecting)  microcracks.  The  overall  elastic  moduli  of  a  statistically  representative  volume 
element  are  then  derived.  In  Section  3,  applications  are  made  to  a  number  of  special  cases.  In  par¬ 
ticular,  for  the  dilute  microcrack  concentration  case,  the  present  approach  recovers  the  well-known 


Taylor’s  model  by  neglecting  interactions  among  microcracks.  When  all  microcracks  are  open  and 
perfectly  randomly  distributed,  the  proposed  approach  shows  that  the  overall  compliance  matrix 
becomes  isotropic.  A  comparison  between  the  present  method  and  the  self-consistent  method  is 
also  discussed  in  this  isotropic  damage  case. 


1.2.  An  ensemble-average  approach  to  microcrack  interaction  and  effective  moduli 


1.2.1.  Background 


Consider  a  statistical  realization  within  a  statistical  RVE  (in  the  probability  space)  of  a  microcrack- 
weakened  solid.  The  volume  of  the  statistical  RVE  is  V,  and  the  exterior  surface  is  subjected  to 
prescribed  traction  tr°^.  Let  a  and  n  be  the  average  microcrack  radius  and  the  average  number  of 
microcracks  per  unit  volume,  respectively.  A  typical  avenue  to  describe  mechanical  responses  of  a 
statistical  RVE  is  to  relate  the  volume-averaged  strain  e  to  the  volume-averaged  stress  It  is  wcD 
known  that 


% 


i  =  S° 


^  (Cul  ®  n  +  n  0  Hull)  dS/ 


(1) 


where  S°  is  the  virgin  elastic  compliance  tensor  for  the  matrix  material;  JuJ  and  n  are  the  micro- 
crack  opening  displacement  vector  and  the  unit  normal  vector  on  discontinuity  surface,  respectively. 
Moreover,  5/  represents  the  union  of  all  discontinuity  (microcrack)  surfaces  (see  Fig.  1). 


Following  the  common  assumption  &  fa  cr°°y  the  microcrack-induced  strain  reads 

r  ss y*  i  (llul  0  n -t- n  0  ttul)  dSj 


(2) 


1.2.2.  Ensemble  average  of  microcrack-induced  strains 

In  this  section,  a  systematic  approach  of  forming  the  ensemble-averaged  strains  and  microcrack 
interaction-induced  local  stress  perturbations  are  presented.  The  basic  idea  behind  this  approach 
is  that  the  local  constitutive  relation  at  a  typical  point  within  a  statistical  RVE  of  a  microcrack- 
weakened  solid  should  be  obtained  by  averaging  over  the  ensemble  of  all  statistical  realizations, 
including  the  locations,  orientations,  lengths  and  relative  configurations  of  randomly  distributed  mi¬ 
crocracks.  This  approach  was  first  applied  to  the  study  of  fluid  suspensions;  see  Batchelor(1970), 
Batchelor  and  Green  (1972),  and  Hinch(1977).  The  ensemble  average  approach  was  later  applied 
to  interacting  inclusions  of  solid  composite  materials  by  Willis  and  Acton  (1976),  and  Chen  and 
Acrivos  (1978a,b).  It  is  emphasized  that  local  displacements,  strains  and  stresses  vary  with  po¬ 
sitions  within  a  RVE.  An  average  over  the  values  of  physical  quantities  occurring  in  a  very  large 
number  of  realizations  is  an  ensemble  average,  which  will  be  denoted  by  angle  brackets. 


Let  us  consider  a  two-phase  composite  statistical  RVE  (composed  of  a  linear  elastic  brittle 
matrix  and  inclusions)  subjected  to  external  load  The  local  constitutive  law  at  a  material  point 
X  may  be  expressed  as 


e(x)  =  S° :  <r(x)  +  e*(x,  C) 


(3) 


where  e  and  er  are  the  local  strain  and  stress,  respectively.  Funhermorc,  e*  is  a  perturbed  strain 
function  which  is  zero  if  x  is  a  point  in  the  matrix  and  is  non-zero  if  x  is  a  point  in  an  inclusion. 
Obviously,  c*  depends  on  the  full  configuration  of  all  inclusions  (denoted  by  C).  The  constitutive 
law  (3)  can  be  easily  statistically  (ensemble)  averaged; 

(6)(x)  =  S°:{tr)(x)-*-(e*)(x)  (4) 

In  addition,  the  ensemble-average  of  c*  at  x  (considered  to  be  non-zero  because  it  lies  in  an 
inclusion  centered  at  xi)  is  (assuming  inclusions  do  not  intersect  one  another): 

(c*)(x)=  /  (c*)(x|xi)/(x,)  dV,  (5) 

Aen. 

where  fl,  is  the  domain  of  a  single  inclusion  and  the  integral  is  performed  over  the  finite  domain 
such  that  X  can  lie  in  an  inclusion  centered  at  Xj.  Further,  f{\\)  is  the  probability  density  function 
(PDF)  for  a  single  inclusion  being  centered  at  xj,  and  (^'’(xjxi))  is  the  perturbed  strain  at  x  averaged 
over  the  subclass  of  realizations  having  an  inclusion  centered  at  Xj. 

At  this  stage,  it  is  reasonable  to  restrict  composite  solids  to  be  locally  homogeneous  [Hinch 
(1977)].  That  is,  all  PDF  do  not  vary  under  small  translation  on  a  macroscopic  length  scale.  With 
this  assumption,  /'xi)  in  (5)  may  be  regarded  as  a  constant  in  the  integration  and  equal  to  /(x). 
The  statistical  “local  homogeneity”  also  allows  a  small  translation  of  x  -  xi  in  the  two  arguments 
of  (f’(x|xi)).  Namely,  we  may  equate  {e*(xlxi))  to  (e'(x  -i-  (x  -  xOjx)).  Thus,  Eq.  (5)  becomes 

(e*}(x)  =  /(x)  /  (e*(x'|x))  dV'  (6) 

dx'en. 

where  x'  =  2x  —  Xi,  and  the  integral  extends  over  all  points  x'  within  an  inclusion  centered  at  the 
position  X.  By  divergence  theorem,  we  have 

(6*)(x)  =  /(x)  /  i(u®n-t-n®  u)(x'lx)dS  (7) 

Jdn,  ^ 

In  the  extreme  case  where  inclusions  become  line  microcracks,  Eq.  (7)  can  be  recast  as 

(0(x)  =  /(x)/  /  i(([u]l®n-»-n0|[uI)(x'|x,a)/(a)d^d5  (8) 

Jg  Jsj  4 

where  Q  =  (a,  n)  characterizes  the  microcrack  length  and  orientation  in  addition  to  the  information 
on  location  x\  f{Q)  is  the  probability  function  for  a  microcrack  with  geometry  Q.  In  what  follows, 
for  simplicity  of  demonstration,  attention  will  be  focused  on  two-dimensional  plane  strain  (or  plane 
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stress)  problems.  Accordingly,  G  defines  the  geometric  domain  (a,  9),  where  a  denotes  one  half 
of  the  microcrack  length  and  6  denotes  the  angle  between  the  global  (reference)  coordinate  and 
the  local  (microcrack)  coordinate  with  0  <  ^  <  tt;  see  Fig.  2.  It  is  well  known  that  for  a  line 
microcrack  in  an  infinite  linear  elastic  isotropic  solid,  the  normal  and  tangential  microcrack  opening 
displacements  take  the  form 

where  E  and  i/  are  the  Young’s  modulus  and  Poisson’s  ratio  of  the  virgin  elastic  solid,  respectively;  p 
and  q  ar^the  normal  and  shear  external  stresses  projected  on  microcrack  surface  in  local  coordinate 
system  (see  Fig.  2).  For  plane  stress  problems,  the  factor  (1  -  i/^)  is  removed  from  (9). 


By  substituting  (9)  into  (8)  and  using  the  Voigt’s  notation  for  strains,  we  arrive  at 


(e*)(x)  = 


=  (T)/(«,  dOda  (10) 


where  A  and  0  are  the  integration  domains  of  microcrack  lengths  and  orientations  (for  open  mi¬ 
crocracks),  respectively.  Further,  g  is  the  transformation  matrix  relating  the  global  and  local  coor¬ 
dinates,  and  T  is  local  stress  vector: 


[gl  = 


2  sin^  9  —  sin  29 

2  cos^  9  sin  29 
-2  sin  2^  2  cos  20  J 


(11) 


Following  the  same  arguments  in  deriving  (e*),  the  ensemble-average  stress  field  (in  the  local 
coordinates)  can  be  shown  to  be 

{T)  =  {p"}  +  ({p})hT»  +  (T)  (12) 

where  is  the  unperturbed  local  stress  field  due  to  remote  loading,  and  (T)  is  the  ensemble 
average  of  the  perturbation  in  local  stress  field  due  to  pairwise  microcrack  interactions: 

(f )(x,  a,9)  =  j  (t)(x,  a, 0|x, ,  ci,  0j)/(Xi,  aj,  0,  jx,  o, 0)  d9^daxdxx  (13) 

Here  (1')(x,a,0jxi,ai,0i)  is  the  stress  penurbation  of  a  microcrack  centered  at  x  with  (a,0)  av¬ 
eraged  over  the  subclass  of  realizations  which  have  a  microcrack  centered  at  Xi  with  (oj,  0i).  In 
addition,  /(xi,fli,0i|x,a,0)  is  the  conditional  probability  function  for  finding  a  micipocrack  cen¬ 
tered  at  Xi  with  (ai,0i)  given  one  microcrack  centered  at  x  with  (a,0).  Since  only  the  ensemble 
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average  is  considered,  the  integral  in  ( 1 3)  is  performed  over  the  domain  of  a  probabilistic  (not  phys¬ 
ical)  RVE  ( in  the  probability  space)  which  is  characterized  by  the  two-point  probability  function 
/(xi,ai,^i|x,a,0).  Thus,  the  shape  of  the  physical  specimen  plays  no  role  here.  The  active  (open) 
integration  domain  E  also  depends  on  loading  conditions  and  ranges  of  microcrack  interactions. 

Assuming  that  microcracks  do  not  intersect  one  another  and  that  reasonable  randomness  holds 
(Hinch  (1977)),  then  /(xi,ai,  ^i|x,a,d)  is  simplified  to  /(xi,ai,^i)-  Further,  by  the  localhomo- 
gtfnfi'O’ assumption, /(xi.G  1,^1 )  becomes  approximately  Therefore,  Eq.  (13)  can  be 

approximated  by  the  following  expression: 

(T)(x,a,^)  =  f(x,a,e)  Jjt)ix,a,6\xuau6x)  dOida^dXi  (14) 

The  assumption  that  microcracks  do  not  intersect  one  another  certainly  introduces  some  unknown 
errors.  However,  there  is  no  analytical  solutions  available  to  accommodate  interactions  among 
pairs  of  arbitrarily  intersected  microcracks.  Therefore,  following  usual  simplification  made  in  the 
literature  [Kachanov  (1987),  Kachanov  and  Laures  (1989)],  we  do  not  consider  the  cases  when 
microcracks  intersect. 

1.2.3.  Approximate  analytical  solutions  of  two^microcrack  interaction  problem 

Due  to  enormous  complexity,  it  is  practically  impossible  to  obtain  closed-form  analytical  so¬ 
lutions  of  strongly  interacting  many-microcrack  problems.  The  formulation  of  arbitrarily  located 
and  oriented  many-microcrack  interaction  problems  is  actually  quite  simple;  sec,  e.g.,  Kachanov 
(1987),  and  Kachanov  and  Laures  (1989).  Nevertheless,  the  formulation  involves  IN  (if  two- 
dimensional)  or  3A'  (if  three-dimensional)  linear  system  of  equations,  where  N  designates  the 
number  of  microcracks  in  an  RVE  (say,  N  =  100).  Clearly,  numerical  solutions  of  many-microcrack 
interaction  problems  are  suitable  for  the  Monte  Carlo  simulation  approach  mentioned  earlier  in  Sec¬ 
tion  1.  By  contrast,  the  present  work  intends  to  construct  closed-form  explicit  stress  solutions  (T) 
and  analytical  expressions  of  microcrack-induced  strains  (e*)  and  compliances  (S*)  for  interacting 
micTocracks  within  the  framework  of  the  ensemble-average  approach. 

In  order  to  construct  useful  explicit  analytical  solutions  and  to  gain  simple  physical  insight 
for  interacting  microcracks,  multiple-microcrack  stress  reflections  will  be  neglected  as  the  first 
approximation.  Namely,  we  will  only  consider  local  stress  perturbations  based  on  many  (arbitrary) 
pairwise  microcrack  interactions.  The  extension  to  the  higher-order  ensemble-average  formulation 
will  be  presented  in  Section  4  of  Part  II. 
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The  exact  analytical  solutions  of  the  boundary-value  problem  of  two  arbitrarily  located  and 
oriented  microcracks  embedded  in  an  infinite  linear  elastic  isotropic  solid  are  not  yet  available  (ex¬ 
cept  for  some  special  configurations  such  as  collinear  microcracks).  However,  approximate  analyt¬ 
ical  solutions  were  proposed  by  several  investigators  [e.g.,  Horii  and  Nemat-Nasser  (1985),  Hori 
and  Nemat-Nasser  (1987),  Chudnovsky  and  Kachanov  (1983),  Chudnovsky  et  al.  (1987a,b),  and 
Kachanov  (1987)].  In  the  present  work,  the  “pseudo-traction”  concept  is  adopted  to  find  approxi¬ 
mate  analytical  solutions  of  the  two-microcrack  interaction  problem.  For  mathematical  simplicity, 
only  ihc  first  term  of  Taylor’s  expansion  of  the  local  stress  field  v^ill  be  used  to  represent  the  average 

stress  across  the  microcrack  line. 

% 

The  local  coordinate  systems  1  and  2  employed  in  the  two-microcrack  interaction  problem  arc 
given  in  Figures  3(a)  and  3(b).  The  two  microcracks  have  lengths  2a\  and  2a2,  respectively.  The 
y\  —  and  ^2— directions  are  set  to  be  normal  to  the  microcrack  lines  C\  and  Ci.  The  original  two- 
microcrack  problem  is  decomposed  into  a  homogeneous  problem  and  two  sub-problems  1  and  2; 
see  Fig.  4.  In  the  homogeneous  problem,  an  infinite  solid  without  any  microcrack  is  subjected  to 
applied  stresses  at  infinity.  In  the  sub-problem  j  (j  =  1, 2),  an  infinitely  extended  solid  under  zero 
remote  stress  at  infinity  has  only  one  microcrack  j,  on  which  the  boundary  conditions  arc 

-P:-^pT-^P3=0  -q,+q^-^q,^0  onC,,;  =  l,2  (15) 


The  quantities  pj  and  pj  are  to  be  determined  in  such  a  way  that  all  boundary  conditions  of  the 
original  problem  are  satisfied.  In  the  subproblem  j,  the  “exact”  stresses  are  given  by  [Sneddon  and 
Lowengrub  (1969)] 

ai,  =  p,(- 1  +  -  F,)  +  q,(2G,  -  //;) 

+  +  (16) 
~  P]^]  9j(“l  +  Ej  —  Fj) 


where 


E, 

E, 

G: 


r  6j\  +  0.2 

COS(0jO  -  -  ■  >  "■) 


ra 


.  3. 


sin  0,0  sin -(0^1 0,2) 

— =  sin(0,o  -  -^r  -^) 

2 

ra^  3 

sin  0,0  cos -(0,1+ 0,2) 


(17) 
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See  Figure  3  for  the  definitions  of  r^i,  rj2,  Oji  and  9j2. 


In  order  to  satisfy  the  original  boundary  conditions  (15),  it  follows  that 


Pj  ~  ®2(i)  ■  ■  ®2(t) 

=  *2(0  ■  ■  *1(0 


^  ¥j  f  i-,j  =  1,2  (no  sum) 


(18) 


where  ande^(-)  are  the  two  unit  base  vectors  forthe  jth  local  c(X)rdinate  system  but  are  expressed 
in  terms  of  the  ith  Icxal  coordinate  components.  The  subscript  i  signifies  that  a  quantity  is  referred 
to  the  (zj,  y-)  coordinate  system.  Substitution  of  Eq.  (16)  into  Eq.  (18)  yields  the  local  stresses  for 
two  micfocracks: 


Pi 

91 
P2 

92 


r  (^ 

0 

Ol 

02  ■ 

0 

0 

03 

04 

QS 

06 

0 

0 

.07 

Og 

0 

0  J 

+ 


pf 

gr 

pf 

gf 


=  aT,_2+Tr-2 


(19) 


where 

=  —  \  +  E2  +  F2  cos  2(t>  +  H2  sin  2^  ;  a2  =  2G2  sin^  <i>-^  H2  cos 2<i>~  {\  —  E2  +  Ez)  sin  2<^ 


=  H2  cos  2<i>  -  Fi  sin  24>  ;  04  =  (G2  -  Hz)  sin  2<^  -  (1  -  £^2  +  £2)  cos  2(j!> 

=  - 1  +  £*1  +  Fi  cos  2?!'  -  Hi  sin  2<f>  •,  a(,  =  2Gj  sin^  (§  +  /fj  cos  2(^  +  (1  -  £1  +  Fj) sin 24> 


=  H\cos2(i>  + F\sin2<f>  ;  og  =  — (Gi  — /fi)sin2<^--(l  —  £1 +£i)cos2<;!> 


and  <f>  =  6\o  —  620.  It  is  noteworthy  that  the  “microcrack  interaction  matrix”  a  in  Eq.  (19)  actually 
corresponds  to  the  “transmission  factor”  matrix  A  in  Kachanov  (1987).  In  addition,  it  is  convenient 
to  define  Tj_2  =  (pi,9i,P2,92)^,  Tfl2  =  (pT,gf,pf,gff  andti_2  =  (pi,gi,  Since 

T]_2  =  T~2  ■♦■Ti-2,  one  can  solve  Ti_2  from  Eq.  (19): 


f,_2  =  a  (I-a)-’  Tr.j 


(21) 


The  comparison  of  the  stress  transmission  factors  obtained  by  using  the  simple  “first  term”  ap¬ 
proximation  (A)  and  the  exact  solutions  (Aei)  for  two  collinear  microcracks  is  given  in  Appendix  I 
for  various  normalized  microcrack-tip  distance  ratios.  When  d  -  0.25,  it  means  that  the  ratio  of 
the  two  microcrack-tip  spacing  to  the  microcrack  size  (2a)  is  0.25.  It  is  noted  that  A  =  1  when 
microcrack  interaction  is  totally  neglected.  From  Table  2  in  Appendix  I,  it  is  observed  that  the 
simple  “first  term”  approximation  is  quite  acceptable  in  general  except  when  two  microcrack  tips 
arc  really  close  to  each  other  (i.c.,  when  d  <  0.1).  It  is  emphasized  that  the  spatial  location  of  the 


second  microcrack  given  the  first  microcrack  is  random.  Therefore,  the  errors  (in  Table  2)  asso¬ 
ciated  with  the  proposed  approximate  analysis  should  be  statistically  averaged  over  all  possible 
realizations. 

Let  us  recall  that  T  =  and  define  Ki  as  the  first  two  rowj  of  the  matrix  fa  •  (l~a)~*]. 

Therefore,  from  Eq.  (21),  we  have 

f  =  K,TJl2  (22) 

We  can  now  substitute  Eq.  (22)  into  (14)  and  (12)  to  find  the  ensemble-average  stress  of  a  primary 
microcrack  located  at  x  over  all  possible  positions  and  orientations  (xi  and  ^i)  of  the  second  neigh¬ 
boring  microcrack.  In  addition,  the  local  and  global  remote  stresses  T“  2  and  t°°  are  related  by  the 
following  transformation  matrix  K2; 

sin^  6  cos^  6  —  sin  20  I 

sin 20  i  sin 20  cos 20  j  <^11  ) 

sin^(0  -I-  <j>)  cos^(0  +  4>)  -  sin  2(9  +  4))  \Z^\  ^ 

sin  2(0  +  4>)  1  sin  2(0  +  4>)  cos 2(0  +  ^)  J  ' 

Therefore,  Eq.  (22)  can  be  rephrased  as 

f  =  K,K2r«^  =  K.r'»  (24) 


where  K  =  Kj  •  K2.  The  detailed  exp/ic/r  components  of  matrix  K  are  given  in  Appendix  11. 

Remark  2.1.  It  appears  that  there  are  some  typo-errors  in  Sneddon  and  Lowengrub’s  (1969) 
solutions.  Eqs.  (16)  and  (17)  are  the  correct  formulas  after  the  problem  is  resolved.  | 

1.2.4.  Overall  moduli  of  brittle  solids  with  interacting  microcracks 

By  substituting  Eqs.  (24)  and  (22)  into  Eq.  (12)  and  assuming  that  x  and  (a,  0)  are  statistically 
independent,  one  obtains  the  following  expression  for  the  ensemble-average  stresses: 

(T)  =  (T“  -h  t>  =  (Ko  +  /(x)/(a,  0)(K))  •  (25) 


where  T°°  =  Kq  •  r°°,  and 


sin^  0  cos’  0  -  sin  20 

-  =  sin  20  1  sin  20  cos  20 

12  2  J 


K  rdr  d9\Q  ddio  dai 


(26) 

(27) 
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Substitution  of  (25)  into  (10)  then  yields  the  ensemble-average  damage-induced  strain  vector: 


(e*)(x)={(S*')  +  (S-^)}.r~  =  (S-)-r- 


(28) 


where 

(S-i)  =  g  •  Ko/(a,  9)  dO  da  (29) 

(S-2)  =  (K)/"(a, 6)  de  da  (30) 

Clearly,  (S*)  is  the  ensemble-average  damage-induced  compliance  matrix,  /(x)  and  f{a,9)  de¬ 
pend  on  material  microstructures  and  loading  conditions,  and  should  be  specified  by  experimental 
observations  (e.g.,  by  scanning  electron  microscopy  or  computerized  tomography). 


In  what  follows,  as  an  illustration  for  deriving  an  explicit  form  of  overall  moduli  of  micro- 
cracked  solids,  we  consider  the  special  case  in  which:  (1)  f{a,6)  =  f(a)fi9),  (2)  microcrack 
orientation  is  perfectly  random  (i.e.,  f(9)  =  1  /tt),  (3)  the  initial  microcrack  radius  a  is  a  constant, 
and  (4)  /(x)  satisfies 


^  fix)  dA 


[fmdA]  lA 


Moreover,  it  is  realized  that 


/(X)  dA  ^  N 


(31) 


(32) 


where  N  is  the  total  number  of  microcracks  within  an  RVE  (with  area  A).  Since  statistical  ho¬ 
mogeneity  (in  the  probability  space)  is  assumed,  the  volume-average  coincides  with  the  ensemble- 
average.  The  volume-average  (denoted  by  an  overline)  of  (28)  over  an  RVE  thus  takes  the  form 

(1  - 


(e-)  =  ju' g  ■  Ko  id  +  jf  g  ■  (K)  rfsj 


where 


rdrd9}od62Q 


(33) 


(34) 


and  E  is  the  active  (open)  probabilistic  integration  domain  for  variables  r,  Oio,  and  Bjo-  The  effective 
radius  r  has  been  normalized  with  respect  to  the  microcrack  radius  a.  Further,  u?  =  ^  =  na^  = 
CjIv  is  a  measure  of  microcrack  concentrations.  It  is  emphasized  that  our  definition  of  microcrack 
concentration  u  differs  from  the  commonly  used  definition  (w)  by  a  multiplier  tt  [see  Budiansky 
and  O’Connell  (1976)].  The  integration  domains  0  and  E  also  depend  on  loading  conditions  and 
range  of  microcrack  interactions. 
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By  substituting  Eq.  (33)  into  the  volume-average  of  Eq.  (4),  we  finally  obtain  the  following 
ensemble-volume  averaged  macroscopic  constitutive  relation  un  Voigt’s  notation): 


(e)  =  (S) .  r 


(35) 


where 


(S)  =  S° -»■  (S*>) -»■  (S-2) 


That  is,  (S)  is  the  ensemble-volume  averaged  overall  effective  compliance  matrix  of  a  microcracked 
brittle  solid.  Specifically,  under  the  plane  strain  assumption,  we  write 


(!+«/) 


1  —  1/  —1/  0 

—  1/  \  —  V  0 

0  0  2 


(S*»)  = 


(1  - 


(S*2)  = 


(1  -  1/2) 


-  1/^)  f 

Hu;  /  g.Kofl 

^uP-  I  g  •  (K) 
Je 


It  is  noted  that  S°  denotes  the  elastic  compliance,  (S*‘)  defines  the  first-order  microcrack  effects 
without  interaction,  and  (S*2)  defines  the  second-order  (in  microcrack  interaction  effects. 


Remark  2.2.  It  is  interesting  to  notice  that  Eq.  (28)  reveals  that  physical  nonlocal  effects 
are  brought  into  constitutive  equations  through  the  process  of  microcrack  interaction  and  ensemble 
averaging  (not  volume  averaging).  That  is,  in  Eq.  (28),  the  stress-strain  law  at  a  material  point 
within  a  RVE  depends  not  only  on  the  constitutive  behavior  at  the  point,  but  also  on  the  constitutive 
behavior  of  neighboring  points.  | 
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1.3.  Some  special  stationary  examples 

In  this  section,  we  apply  the  proposed  ensemble-volume  average  approach  to  solve  three  special 
examples.  These  examples  are  all  classified  as  stationary  micromechanical  models.  In  particular, 
the  first  example  shows  that  the  proposed  method  recovers  the  well-known  Taylor’s  model  when 
microcrack  interaction  is  entirely  neglected  in  Eq.  (36).  In  the  second  example,  we  examine  the 
randomly  located  aligned  (parallel)  microcracks  problem  —  for  transverse  cracking  phenomenon 
in  laminated  composites,  for  instance.  The  last  example  shows  that  the  overall  compliance  ma¬ 
trix  becomes  isotropic  when  all  microcracks  are  open.  Microcrack  locations  and  orientations  are 
assumed  to  be  perfectly  random.  In  what  follows,  for  simplicity,  we  further  assume  that  closed 
microcrack  contributions  are  neglected. 

I.3.I.  Dilute  microcracks  —  Taylor’s  model 

The  open  microcrack  angle  domain  0  =  [61,62]  can  be  found  by  solving 

sin^  ^<7^  +  cos^  6(7^  —  sin  26a^  >  0  (40) 


where  tensile  normal  stress  is  taken  as  positive.  After  dropping  the  interaction  term  in  (36),  the 
overall  complia.ice  becomes 

^  =  (41) 


where  [by  integrating  Eq.  (38)] 

■s;l  0 

0  S22 

c»i  c*i  c*l 

-  -61  -^62  -^65 

Specifically,  we  have 

Sj*,'  =62  —  61  —  ;^(sin  262  —  sin  26\)  ;  5^2  =  ^2  -  +  r(sin  262  —  sin  2^0 

2(02-01)  ;  5[^=sin2  0,-sin2  02  =  56;  =5^  =  56* 


(42) 


(43) 


For  various  biaxial  loading  conditions,  we  obtain  the  following  results; 

(1)  Uniaxial  or  biaxial  tension.  All  microcracks  are  open  (0i  =  0, 02  =  ’r)  and  isotropic  damage 
is  recovered.  The  orientation  domain  is  0  =  [0,  ▼]  The  non-zero  components  in  (43)  are 

■5u=52’=7r;  S^=27r  (44) 
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(2)  Biaxial  tension/compression.  We  have  the  axial  stress  >  0,  the  lateral  stress  cr*,  <  0, 
and  shear  stress  =  0.  The  open  microcrack  angle  domain  is  definea  by  0  =  (0, 0*]  U  [tt  - 
where 


0*  =tan-' 


'22 


’ll 


(45) 


The  non-zero  anisotropic  components  in  (43)  are 


S;l  =  2ff‘  -  sin  26»*  ;  Sjj  =  20'  +  sin  20’  ;  5;^  =  40'  (46) 


(3)  Biai^ai  tension  and  shear.  We  have  <7^’  =  <^22  >  0  and  cr^  >  ffg.  The  open  domain  is 
defined  by  0  =  [f,7r]  u  [0.0’]  u  (|  -  where 


The  non-zero  components  in  (43)  are 

5*,'  =  1  +  20*  =  Sjj  ;  =  TT  +  40' 

s\l^cos20'^s;l^s'^  =  s;l 


(47) 


(48) 


1.3.2.  Aligned  microcracks 


Attention  is  now  focused  on  an  array  of  randomly  located  aligned  (parallel)  open  microcracks. 
This  class  of  problems  is  relevant  to  the  transverse  matrix  cracking  problem  in  laminated  compos¬ 
ites;  see,  for  example,  Laws,  Dvorak  and  Hejazi  (1983),  and  Laws  and  Dvorak  (1987).  Since  all 
microcracks  have  same  orientation,  Eq.  (17)  can  be  simplified  as  follows: 


E,  =  £2  = 

Gi  -  G2  = 

where 


'■  cos(«,„  -  ;  F,  =  = 


r 

\/r^ 


TQ^  3 

- - ^  sin  010  sin -(011  +^12) 


sin(«.„  -  ;  ff,  =  W,  = 


^  3 

sin 010 cos -(011 -(-012) 
2 


rj  =  (r^  _  2arcos0io -(- ;  r2  =  (r^  +  2ar  cos0io  + 

^  _i  /^rcos0io  -  a\  ^  _i  /'rcos0io  +  a\ 
011  =  cos  ’  ( - - - \  ;  012  =  cos  *  I - - - 1 


and  0  =  0  =  4)  =  0io  —  020. 0io  =  02o- 


(49) 


(50) 


The  numerical  integration  of  Eq.  (34)  will  be  addressed  in  Section  3  of  Part  II  of  this  series. 
The  minimum  radius  of  integration  in  r  (center  distance  between  microcracks)  is  taken  as  2a.  The 
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Table  1.  Convergence  behavior  of  '  vs.  r, 


(A'la) 

(A'23) 

4a 

0.429843+00 

0.173965+00 

8a 

0.550972+00 

0.225418+00 

16a 

0.582654+00 

0.238972+00 

20a 

0.586498+00 

0.240618+00 

40a 

0.591638+00 

0.242820+00 

400a 

0.591894+00 

0.243244+00 

convergence  .avior  for  the  integration  in  (34)  (i.e.,  {Kij))  with  respect  to  the  choice  of  the  max¬ 
imum  integration  radius  rm,,  is  given  in  Table  1.  It  is  observed  that  use  of  r^tx  =  20a  is  qi  'te 
acceptable.  After  numerical  Gauss  integration  of  (34)  (by  using  40  Gauss  points),  we  obtain 


{K)  = 


0  0.5865 
0  0 


0 

0.2406 


Accordingly,  the  transversely  isotropic  compliance  becomes 

r  i  0  1 

7^  1  0 


(S)  = 


(1  - 


0 


0  ifrj 


+  ir 


0  0  0 

0  2w  +  1.173u;^  0 

0  0  2w  +  0.481a>^ 


(51) 


(52) 


Clearly,  the  first-order  terms  uj  in  Eq.  (52)  come  from  Taylor’s  model,  and  the  second-order 
corrections  in  (52)  come  from  the  pairwise  microcrack  interactions. 

1.3.3.  Isotropic  damage 

Under  uniaxial  or  biaxial  tension,  all  microcracks  are  typically  assumed  open.  It  is  also  as¬ 
sumed  that  microcracks  are  randomly  located,  randomly  oriented,  and  of  unifcMm  size.  As  a  result, 
the  overall  behavior  is  isotropic.  The  compliance  contribution  due  to  microcrack  interaction  can 
be  evaluated  by  performing  numerical  Gauss  integration  of  Eq.  (39).  The  results  are: 


(S-2)  = 


(1  - 1'^) 


■TTu; 


0.373  0.041 
0.041  0.373 
0  0 


0 

0 

0.664 


(53) 


(S)  = 


Therefore,  the  (plane  strain)  overall  compliance  becomes 

1  e  0 1 

1  0  +  TT 

0 


(1  -  */^) 


0 


2 

1-1.  J 


w-K  0.373^2  0.04 1  u>2 

0.04  la;2  w  +  0.373u.’2 

0  0 


0 

0 

2u;  +  0.664w^ 


(54) 
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It  is  noted  that  the  shear  moduli  (in  1-2  and  2-1  directions)  are  changed  (non-zero)  due  to  micro- 
crack  interaction. 


It  is  very  interesting  to  compare  the  foregoing  overall  compliance  matrix  with  that  derived  by 
the  self-ccmsistent  model  [e.g.,  Horii  and  Nemat-Nasser  (1983)].  The  overall  compliance  given  by 
the  self-consistent  method  in  this  isotropic  damage  case  is: 


S’"  = 


(1  - 


I— irw 
—1/ 
1-v 


0 


0 


0  • 
0 

2  I  -"yt/w 
I  — t/  i — nw  - 


(55) 


Obviously,  the  overall  compliance  given  by  the  self-consistent  method  goes  to  infinity  atw  =  l/rr. 
It  should  be  noted  that  any  effective  medium  theory  is  inherently  only  applicable  for  low  microcrack 
concentrations  since  it  does  not  depend  on  microcrack  locations. 


For  low  values  of  u;  <  <  1  /rr,  the  Taylor’s  expansion  of  Eq.  (55)  renders 


S""  = 


(1  -  1/2) 


i 

■  1 

—  If 

\-tf 

0  ■ 

u  +  7ru.’2 

0 

0 

—  If 

\-v 

1 

0 

+  TT 

0 

u  +  7ru;2 

0 

\ 

0 

0 

2 

l-u  J 

0 

0 

2(u>  7ru;2) 

(56) 


Apparently,  the  compliance  calculated  by  the  self-consistent  method  is  higher  than  that  of  the 
present  model.  In  addition,  the  self-consistent  method  predicts  no  change  in  shear  moduli  in  the 
1-2  and  2-1  directions.  Finally,  it  is  err.phasized  that  there  is  no  singularity  in  the  present  method 
irrespective  of  the  u;  value. 


From  Eq.  (54),  it  is  clear  that  the  relative  weight  of  the  second-order  terms  (in  due  to 
microcrack  interaction)  to  the  first-order  terms  (in  w)  depends  on  the  values  of  u;.  For  example, 
the  relative  weight  is  0.373w  for  the  11  and  22  terms  in  the  compliance  matrix.  If  the  microcrack 
density  is  very  small  (dilute),  the  second-order  terms  are  negligible.  On  the  other  hand,  if  the 
microcrack  density  is  high  (e.g.,  u;  =  0.5  -  corresponding  to  u;  =  O.Stt  defined  by  Budiansky  and 
O’Connell  (1976)],  then  contributions  from  the  second-order  terms  arc  quite  significant.  It  is  noted 
that,  however,  when  microcrack  densities  arc  extremely  high,  the  proposed  second-order  method 
may  not  be  accurate  enough  since  higher-order  terms  should  be  included;  see  Sec.  4  of  Part  11  for 
extension  to  the  higher-order  formulation. 


The  proposed  statistical  micromechanical  framework  is  not  restricted  to  widely  spaced  micro- 
crack  arrays.  From  Table  2  in  Appendix  I,  it  is  seen  that  the  proposed  two-microcrack  interaction 
solutions  are  quite  acceptable  even  when  the  ratio  d  is  only  0.1  (i.e.,  when  the  “microcrack-tip 
spacing”  is  only  one-tenth  of  the  “microcrack  size”). 


1.4.  Conclusions 


A  two-dimensional  statistical  micromechanical  theory  for  microcrack-weakened  brittle  solids 
is  presented  based  on  the  concepts  of  ensemble-average  and  pairwise  microcrack  interaction.  The 
overall  compliances  are  derived  by  performing  ensemble-volume  averaged  integration  over  the 
domain  of  a  statistically  representative  volume  clement  (in  the  probability  space).  Unlike  the  self- 
consistent  method  which  introduces  a  singularity  at  u?  =  1  /  jt,  our  method  predicts  a  smooth  increase 
of  effective  compliance  as  u  increases.  Further,  the  proposed  microcrack  interaction  framework 
docs  not  require  the  use  of  Monte  Carlo  simulations.  It  is  emphasized  that  the  proposed  second- 
order  teivns  are  second-order  in  not  second-order  in  geometric  quantities  [such  as  (a/r)^].  In 
addition,  pairwise  interaction  effects  do  not  vanish  in  the  event  of  perfectly  randomly  distributed 
microcracks;  seeEq.  (53).  In  fact,  effects  of  pairwise  microcrack  interactions  are  equally  distributed 
in  all  directions  and  therefore  the  overall  behavior  is  isotropic. 

In  Part  II  of  this  series,  we  will  further  investigate  issues  pertaining  to  computational  aspects  of 
the  proposed  approach,  extension  to  “evolutionary”  micromechanical  damage  models  (with  “cleav¬ 
age  1”  microcrack  growth),  and  a  higher-order  microcrack  interaction  model.  Extension  to  a  three- 
dimensional  statistical  micromechanical  theory  of  brittle  solids  with  randomly  located,  interacting 
penny-shaped  microcracks  is  presented  in  a  separate  pap<ir  [Ju  and  Tseng  (1992)].  Numerical  com¬ 
parisons  with  the  Taylor’s  model,  the  self-consistent  model,  and  the  differential  scheme  will  also 
be  presented  in  Ju  and  Tseng  (1992). 
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1.6.  Appendix  I.  Average  stress  field  of  two  coUinear  microcracks 


The  problem  of  two  interacting  collinear  microcracks  of  equal  length  in  an  infinite  plate  sub¬ 
jected  to  uniform  traction  at  infinity  has  an  exact  analytical  solution  (Sneddon  and  Lowengnib 
(1969)].  The  normalized  (nondimensional)  distance  between  the  two  microcrack  tips  is  denoted  by 
2k  (with  A:  <  1),  and  the  normalized  length  for  each  microcrack  is  designated  as  1  -  For  this 
special  configuration,  we  have 


^10  =  ^20  =  ^ij  =  0 ;  t,  j  =  1, 2 


(57) 


and  <f>  =  Sxo-  (?2o  =  0 


,  ,  \  +3k  3  +  fc 

r  =  l+  k-,  ri  =  -  —  =  ri,  =  r2i  ;  r2  =  =  rj2  =  r22 


Therefore,  the  only  non-zero  components  in  (17)  are 


A  =  £1  =  £2  = 


2(1  -t-  k) 

Vil  +  3fc)(3  +  k) 


(58) 


(59) 


Substitution  of  Eq.  (59)  into  (20)  and  solution  of  Eq.  (19)  then  render  the  approximate  average 
normal  stress  p  as  follows 

(60) 

The  exact  solution  for  this  average  normal  stress  takes  the  same  format  as  in  (60),  with  A  replaced 
by  Aei  [exact;  see,  e.g.,  Kachanov  (1987)] 


_  y/lTlk 

1  +  v/X 


(61) 


Table  2  shows  the  comparison  between  A  and  Aer  for  different  values  of  d  which  is  defined  as 
d  ~  2k I (1  —  k).  It  is  clear  that  the  error  associated  with  the  approximate  solution  increases  as  d 
decreases.  This  is  due  to  the  stress  field  singularity  at  microcrack  tips. 
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Table  2.  Comparison  between  A  and  Acx 


d 

exact 

approx 

error(%) 

10 

1.001037 

1.001035 

0.0002 

5 

1.00352 

1.00349 

0.0025 

1 

1.03528 

1.03280 

0.2396 

0.5 

1.07047 

1.06066 

0.9161 

0.25 

1.118034 

1.091089 

2.41 

0.1 

1.18821 

1.122683 

5.51 

0.05 

1.23801 

1.13721 

8.1419 

0.005 

1.34863 

1.15280 

14.52 

0.0005 

1.39238 

1.15451 

17.08 

i 


(62) 


(63) 


(64) 


i 


» 


i 


I 


I 
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1.7.  Appendix  11.  Components  of  matrix  K  in  Eq.  (24) 


where 


and 


Ku  =  k\sin2d  +  kzsin^ 0 +  ki  ;  K\2  = -kisinlO +  k2COS^ 9 +  k^ 

K\2  =  —kz  sin  26  —  2k\  cos  26  ;  Kzi  -  ^4  sin  26  +  ks  sin^  6  +  k(, 

Kzz  =  -ki  sin  26  +  ks  cos^  ^  +  ^6  ;  K23  ==  -ks  sin  26  -  2^4  cos  26 

ki  =  ^(/2  sin  2(f>  -  f\ cos  2(f>  -  fs)  ;  kz  =  ^(/z cos  24>  +  /j  sin  24>  +  /») 
*^3  =  ^(/2  sin^  sin  2()))  ;  ^4  =  sin  2<i>  -  fs  cos  24>  -  fj) 

ks  =  ^(/6  cos  2<i>  +  fs  sin  24>  +  f%)  ;  kt  =  ^(/s  sin^  <(>  -  jfs  sin  2<i>) 

fl  ~  (0104  -  QzQ2)Q6  +  0-2  ;  /z  =  — (O1Q4  -  01203)0%  +  O] 
fs  =OtiQ6+  0-208  ;  /4  =  Osfz  +  O7/1 

/s  =  -(0104  -  0203)05  +  04  ;  /6  =  (Q:iC.-4  -  0203)07  +  03 
/?  =  og/s  +  Ofi/fi  ;  /g  =  03O5  +  0407  ;  A  =  <iet(I  -  a) 
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L8.  Figure  captions 


I  Figure  1.  Superposition  of  microcrack  problem. 

Figure  2.  The  local  and  global  coordinates  systems. 

I  Figure  3(a)(b).  Definitions  of  geometric  quantities  for  two  neighboring  microcracks. 

Figure  4.  Decomposition  of  two-microcrack  interaction  problem. 

I 
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Figure  2.  The  local  and  global  coordinates  systems . 
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Figure  3(a)(b).  Definitions  of  geometric  quantities  for  two  neighboring  miciocracks 
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PART  II 

Effective  Elastic  Moduli  of  Two-Dimensional 
Brittle  Solids  with  Interacting  Microcracks. 
II :  Evolutionary  Damage  Models 


n.O^  Abstract 

In  Part  I  of  this  series,  basic  formulations  of  stationary  micromechanical  theory  and  overall 
responses  are  presented  for  two-dimensional  brittle  solids  with  randomly  dispersed  microcracks. 
The  basic  formulations  hinge  on  an  ensemble  average  approach  which  includes  pairwise  micro¬ 
crack  interactions.  In  this  paper,  statistical  micromechanical  evolutionary  models  arc  proposed  to 
account  for  “cleavage  1”  growth  of  randomly  oriented  and  located  microcracks  under  microcrack 
interaction  effects.  Biaxial  tension/compression  loadings  arc  also  considered  to  take  into  account 
mixed  microcrack  opening  and  closure  effects.  Efficient  numerical  integration  algmithms  for  the 
proposed  ensemble  averaged  constitutive  equations  arc  subsequently  given.  Further,  uniaxial  and 
biaxial  tests  are  presented  to  illustrate  the  proposed  models  and  procedures.  Finally,  a  higher-order 
mictocrack  interaction  model  within  the  proposed  micromechanical  framework  is  discussed. 


n.l.  Introduction 


A  statistical  micromechanical  ensemble-volume  average  approach  to  explicitly  derive  overall 
constitutive  equations  for  microcrack-weakened  brittle  solids  has  been  presented  in  Part  I  of  this 
series.  The  proposed  explicit  micromechanical  approach  takes  into  account  interactions  among 
microcracks,  random  microcrack  orientations,  locations  and  densities,  and  microcrack  opening  dis¬ 
placements.  Damage-induced  overall  anisotropy  can  be  predicted  by  the  proposed  framework.  The 
proposed  micromechanical  approach  is  at  variance  with  existing  “effective  medium  theories”  such 
as  the  self-consistent  method  [Hill  (1965),  Budiansky  and  O’Ckjnnell  (1976),  Horii  and  Nemat- 
Nasser  ^1983)],  the  generalized  self-consistent  method  (Christensen  and  Lo  (1979)),  the  differen¬ 
tial  scheme  [Roscoe  (1952,  1973),  McLaughlin  (1977)],  and  the  Mori-Tanaka  method  [Mori  and 
Tanaka  (1973),  Benveniste  (1986)],  etc..  In  addition,  the  proposed  framework  is  different  from 
those  micromechanical  interaction  analyses  based  on  Monte  Carlo  simulations  and  numerical  com¬ 
putations  of  stress  “transmission  factors”  [Kachanov  (1987),  Kachanov  and  Laures  (1989)).  The 
simple  approximate  analytical  solutions  for  the  two-microcrack  interaction  problem  render  closed- 
form  explicit  expressions  for  the  evaluation  of  ensemble-average  integrals. 

Nevertheless,  during  loading  histories  of  brittle  solids,  pre-existing  microcracks  may  become 
unstable  and  grow  along  certain  preferred  orientations,  depending  on  loading  levels,  microcrack 
configurations  and  interactions,  etc..  Therefore,  it  is  important  to  formulate  micromechanical 
“cleavage  1”  [Ashby  (1979)]  evolutionary  damage  models  to  account  for  both  microcrack  kinetics 
and  microcrack  interactions  under  specified  loads.  Some  micromechanical  “evolutionary”  dam¬ 
age  models  were  proposed  in  the  literature  for  “effective  medium  theories”;  see,  e.g.,  Krajcinovic 
and  Fanella  (1986),  Fanella  and  Krajcinovic  (1988),  Sumarac  and  Krajcinovic  (1987),  Ju  (1991), 
Ju  and  Lee  (1991),  and  Lee  and  Ju  (1991),  etc..  However,  these  micromechanical  evolutionary 
damage  models  are  only  valid  for  weak  microcrack  interactions  under  low  or  at  most  moderate 
microcrack  concentrations.  On  the  other  hand,  microcrack  interaction  models  due  to  Kachanov 
(1987),  Kachanov  and  Laures  (1989)  do  not  consider  microcrack  growth  or  nucleation  kinetics. 

The  primary  objective  of  this  paper  is  to  extend  the  framework  proposed  in  Part  I  to  accommo¬ 
date  “cleavage  1”  microcrack  growth  process  under  microcrack  interactions.  “Cleavage  1”  process 
implies  that  only  the  growth  of  pre-existing  microcracks  is  considered.  Therefore,  microcrack 
number  remains  constant  during  loading  histories.  Of  course,  one  may  also  consider  “cleavage 
2”  process  [Ashby  (1979)]  to  accommodate  nucleation  of  new  microcracks.  However,  that  is  be¬ 
yond  the  scope  of  the  present  paper.  In  Section  2,  micromechanical  kinetic  equations  characteriz¬ 
ing  the  “process  domains”  of  active  (open)  microcracks  are  introduced.  These  “process  domains” 


together  with  “open  microcrack  domains”  completely  define  the  integration  domains  of  the  en¬ 
semble  averaged  constitutive  equations  relating  macro-strain  and  macro-stress.  Moreover,  various 
tension/compression  loadings  are  presented  in  Section  2  to  illustrate  the  proposed  approach.  Effi¬ 
cient  numerical  integration  algorithms  are  presented  in  Section  3  to  implement  the  proposed  dam¬ 
age  models.  A  number  of  uniaxial  and  biaxial  tests  are  also  performed.  Finally,  discussions  on  a 
higher-order  microcrack  interaction  model  within  the  proposed  framework  is  given  in  Section  4. 


n.2.  Evolutionary  models  with  microcrack  interaction 


Let  us  consider  an  ensemble  of  randomly  oriented  and  located  microcracks  within  a  statistically 
representative  volume  element.  The  microcrack  number  density  is  denoted  by  n  and  the  micro- 
crack  length  is  denoted  by  2a.  The  representative  element  is  subjected  to  two-dimensional  external 
loading  r'”  =  erg,  In  accord  with  Part  I  of  this  work,  the  overall  compliance  matrix 
(S)  can  be  written  as  (plane  strain  with  linear  elastic  matrix) 


(S)  =S°-h(S*>)  +  (S*2) 


(1) 


where  [Sfee  Eqs.  (37)-(39)  in  Part  I] 


and  [see  Eq.  (34)  in  Part  I] 


(1  -H/)  f  ^  ^  ® 

=  ^  ^  ^  1/  1  1/  0 

^  L  0  0  2. 

(2) 

(S«1)  =  n  J^a]g.Kod9 

(3) 

(4) 

(k}=  (Krdrd9iod92o 

(5) 

Note  that  a\  is  the  half-length  of  a  primary  microcrack  and  definitions  of  g,  K  and  Ko  are  given 
in  Eqs.  (11),  (24)  and  (26)  of  Part  I,  respectively.  The  effective  radius  r  has  been  normalized  with 
respect  to  the  microcrack  radius  oi.  Further,  E  is  the  current  active  integration  domain,  including 
contributions  from  open  stationary  or  unstable  microcracks.  Contributions  due  to  closed  micro¬ 
cracks  are  not  considered  in  this  paper. 


For  biaxial  loadings,  the  open  microcrack  domain  can  be  approximately  defined  by  finding  the 
domain  of  tensile  (positive)  normal  stresses  on  microcrack  planes  as  follows 


org  sin^  0  +  <T^  cos^  9-0^  sin  29  >0 


(6) 


It  is  noted  that  Eq,  (6)  only  considers  the  far-ficld  stresses  tr°°.  In  some  cases,  however,  the 
open  microcrack  domain  can  be  affected  by  microcrack  interaction.  That  is,  some  originally  open 
microcracks  may  become  closed  due  to  strong  microcrack  interaction.  When  this  occurs,  the  local 
stresses  tr  should  be  used  in  (6);  i.e.,  one  should  check  the  following  “local  (nonnal)  opening 
condition”  to  ensure  that  a  microcrack  is  indeed  open:  {p)  >  0. 


n.2.1.  Kinetic  equations  for  microcrack  growth 

In  continuum  damage  mechanics,  kinetic  equations  characterizing  damage  processes  arc  typi¬ 
cally  postulated  on  the  basis  of  hypothetical  phenomenological  arguments  and  macroscopic  obser¬ 
vations.  However,  as  shown  by  Krajcinovic  and  Fanella  (1986),  microcrack  kinetic  equations  can 
be  actually  derived  based  on  micromechanics,  probabilistic  configurations  of  microcracks  along 
certain  microstructural  weak  planes,  and  fracture  energy  barriers  on  the  mesoscale.  When  the  local 
stress  intensity  factor  (or  strain  energy  release  rate)  reaches  a  certam  critical  value  for  the  weak 
plane,  a  pre-existing  microcrack  is  assumed  to  become  unstable  and  grows  along  the  weak  plane 
in  an  unstable  manner  until  reaching  a  certain  characteristic  size  such  as  a  grain  or  facet  size  [Kra¬ 
jcinovic  and  Fanella  (1986),  Sumarac  and  Krajcinovic  (1987)].  This  idealization  is  reasonable  for 
pre-existing  intergranular  microcracks  or  interfacial  microcracks  (between  the  matrix  and  inclu¬ 
sions)  in  certain  brittle  materials  such  as  concrete,  ceramics,  rocks,  ice,  and  some  brittle  compos¬ 
ites,  etc..  For  simplicity,  we  do  not  consider  micixKrack  kinking  (into  the  matrix  or  neighboring 
weak  planes)  nor  nucleation  of  new  microcracks;  sec  Fanella  and  Krajcinovic  (1988),  Ju  and  Lee 
(1991),  and  Lee  and  Ju  (1991)  for  some  recent  treatments.  In  other  words,  it  is  assumed  that  there 
exists  a  higher  energy  barrier  serving  as  a  microcrack  trapping  mechanism;  e.g.,  a  higher  fracture 
toughness  due  to  the  matrix  or  due  to  intergranular  kinking.  Therefore,  the  proposed  evolutionary 
damage  models  are  only  suitable  for  pre-peak  (not  post-peak)  behavior  of  a  class  of  brittle  materials 
aforementioned.  “Run-away  instability”  or  localized  macrocrack  growth  is  not  considered. 

Within  the  context  of  microcrack  interaction,  the  local  stresses  tr  (not  er°°)  should  be  employed 
to  compute  the  stress  intensity  factors  (or  strain  energy  release  rates)  on  pre-existing  microcrack 
tips.  This  treatment  is  completely  different  from  existing  (“effective  medium”)  micromechani- 
cal  evolutionary  models  which  employ  far  field  stresses  <t°°.  Following  Krajcinovic  and  Fanella 
(1986),  we  make  the  following  two-stage  simplification.  It  is  assumed  that  the  length  of  a  pre¬ 
existing  microcrack  can  either  be  Zqq  if  the  microcrack  is  stationary,  or  2<i  /  if  the  microcrack  is 
activated.  Note  that  2ao  denotes  the  average  length  of  an  initial  microcrack  and  2a /  signifies  the 
average  length  of  a  propagated  microcrack  (e.g.,  the  grain  or  facet  size).  The  relation  between  the 
initial  and  final  microcrack  lengths  is  assumed  to  be 


ao  =  Po  Of  (7) 

where  po  is  a  scalar  between  0  and  1.  The  microcrack  growth  from  2ao  to  2a/  is  assumed  to  be 
instantaneous  once  the  local  fracture  energy  is  reached. 
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The  actual  fracture  criterion  of  a  material  depends  on  its  underlying  physics  of  fracture.  This 
paper  does  not  intend  to  investigate  fracture  criteria  of  all  materials.  Instead,  we  attempt  to  demon¬ 
strate  possible  procedures  involved  in  a  micromechanical  statistical  evolutionary  damage  formula¬ 
tion.  For  illustration  purpose,  the  following  mixed  mode  fracture  criterion  is  employed  [see  Kan- 
ninen  and  Popelar  (1985,  p.  50-51)]: 


(8) 


where  Kj  and  A'//  are  the  Mode  I  and  Mode  II  stress  intensity  factors  of  an  initial  microcrack, 
respectively.  Similarly,  Kic  and  Kju,  respectively,  are  the  Mode  I  and  Mode  II  critical  stress  in¬ 
tensity  factors  of  the  weak  plane  where  a  pre-existing  microcrack  is  located.  The  foregoing  fracture 
criterion  may  be  suitable  for  some  materials  but  modifications  are  necessary  when  application  to  a 
specific  material  is  sought.  In  the  spirit  of  statistical  ensemble -average  stresses,  it  is  reasonable  to 
write  the  following  simple  approximations  (under  plane  strain): 


A’/ 


(9) 


where  (p)  and  {q)  =  the  ensemble-averaged  local  normal  and  shear  stresses  projected  on  microcrack 
surfaces  in  local  coordinates.  Substitution  of  (9)  into  (8)  leads  to 


{p?  +  =  Fl 


(10) 


where 


(11) 


For  uniformly  distributed  and  oriented  microcracks,  one  has  /{&)  =  1/ff  and  fix)  =  n.  From 
Eqs.  (25>-(27)  of  Part  I,  one  can  write 

=  (12) 

For  convenience,  let  us  define 

M  =  Ko  +  -(K)  (13) 

TT 

It  is  emphasized  that  (K)  (and  thus  M)  includes  microcrack  interaction  effects.  With  (12)  and  (13) 
at  hand,  Eq.  (10)  can  be  rephrased  as  follows 


1  ^  xoo  .  ^  oo  —  1 

+  74<7ii  +  75<722*^12)  =  * 


(14) 


(15) 


where  =  o-^/Fc  and 

-fo  =  Mn  +  k^Ml  ;  7i  =  A/„Mi2  +  ;  72  =  A/,\  + 

73  =  A/f3  +  A:^A/|3  ;  74  =  AfnA/js  +  l:^A/2i  Af23  i  Is  =  MnMu  ■*■  k^M22M2i 

Eq.  (14)  can  be  used  to  solve  for  the  “unstable  angle  domains”  (9^)  defining  the  “cleavage  1” 
unstable  microcrack  growth  domains  under  a  specified  far  field  loading  <r°°.  Conversely,  for  a 
specified  “unstable  growth  angle  domain”,  Eq.  (14)  can  be  used  to  obtain  the  external  loading. 

n.2.2.  Dilute  case  —  Taylor’s  model 

In  tliis  case,  the  overall  compliance  of  a  microcrack-weakencd  solid  takes  the  form 


(S)  =S°  +  (S*'‘)  +  (S''p) 


(16) 


where  (S*^‘)  and  (S’'p)  denote  the  compliance  contributions  from  initial  microcracks  and  micro¬ 
crack  growth,  respectively.  (S*'p)  can  be  expressed  as  [cf.  (42>-(43)  of  Pan  I]: 

(1  -  i/^) 


(S*>p) 


W0(/50^-  1) 


where  uq  = 


C**p 

0 

C**P' 

■^16 

0 

c*lp 
.3  22 

c»ip 

■^26 

c*ip 

^61 

c»tp 

■362 

C»*P 

■^66  J 

it 

~9\^ 

jfsin  2^2  “  sin  2^1’) 

(17) 


5^"  =  2{9l  -  ^P)  ;  =  sin^  e\  -  sin"  9l  =  s;,'”  =  5^"  = 

In  Eq.  (18),  O''  =  [9^,6^]  is  the  evolutionary  angle  domain  defined  by  Eq.  (14). 


(18) 


(1)  Uniaxial  tension:  ag  >  0,  aff  =  alf  =0.  All  microcracks  are  open  but  overall  responses 
are  anisotropic  due  to  preferred  orientations  of  microcrack  growth.  Specifically,  the  process 
domain  is  0p  =  [0,  0p]  U  [tt  -  tt]  (see  Fig.  1),  where  can  be  solved  from  (14)  for  given 
ag  (or  vice  versa): 

72  =  cos^  9^  +  fc"  cos"  9^  sin"  9^  -  ag"  (19) 

From  (18),  we  have 


5*;'’  =  29^  -  sin  29^  ;  S*]"  =  29^  +  sin  29^  ;  5^"  =  49^  ;  others  =  0  (20) 


(2)  Biaxial  tension/compression:  >  0,  ajt  <  0,  =  0.  The  bound  of  unstable  growth 

angle  domain,  9^,  can  be  computed  from  Eq.  (14).  The  process  domain  is  therefore  O'"  = 


[0, 0^1  u  [tt  -  ff].  Conversely,  one  can  obtain  the  axial  external  tensile  stress  for  given 
lateral  external  stress  and  prescribed  angle  bound  6^: 


The  non-zero  components  of  (S**'’)  are 

s;/''  =  20^  -  sin  20^  ;  Sjj'’  =  20^  +  sin  20^  ;  Sl\^  =  W  (22) 

Since  ^ii  <  0,  it  follows  that  is  negative  and  proportional  to  the  density  cjq-  This  is  due  to  the 
fact  that  Mo(k  II  contribution  to  (contractive)  lateral  strain  is  greater  than  Mode  I  contribution 
fOT  biaxial  tension/compression  loads.  To  wit,  let  us  consider  an  arbitrary  stadonary  microcrack 
and  compare  the  contributions  of  Mode  I  and  Mode  II  to  tj 

("I  =  2— — *1  sin^  0{co^  0a^  +  sin^  (Mode  I)  (23) 

h 

C*(^  =  JJ)  (24) 

^  ti 

The  sum  of  Mode  I  and  Mode  II  contributions  is 

e’u  =  sin^  <  0  (25) 

II.2.3.  Uniaxial  and  biaxial  tension  loadings 

Within  the  proposed  framework,  the  computauon  of  overall  compliances  require,  ’nformarion 
regarding  integration  domains  for  microcrack  opening  and  growth.  For  uniaxial  tension,  all  micro¬ 
cracks  are  initially  open  and  Eq.  (14)  becomes  (19)  in  which  72  is  evaluated  by  (15).  The  length  of 
a  microcrack  is  either  oj  =  qq  i(  stationary  or  oj  =  aj  if  unstable.  For  a  prescribed  configuration 
of  microcrack  growth  domain,  Qp  =  [0,  0p]  u  [tt  -  0^,n],  one  can  systematically  explore  the  status 
and  length  (02)  of  a  second  neighboring  random  microcrack.  The  details  (27  possible  cases)  are  as 
follows;  see  also  Figure  2. 


02 


Category  1.  If  0  <  ^  <  0'’  and  a\  =  aj:9  possible  cases  concerning  the  length  of  02. 
0  <  5  +  ^>,0  <  02  6^  <0  +  010  <  o"  _ 


(1) lf0<02O<0  +  0\o  +  then 

0/ 

(4)  I  f  0  <  020  ^  ®  +  ^10  ~  0^^  then 

ao 

(2)  If  9  +  ^10  +  <  ^20  ^  ^  +  ^10  +  0^ 

ao 

(5)If  0  +  010-0^  <  020  <0  +  010  +  0” 

(3)  I  f  6  +  610  +  0^  <  020  ^  ^ 

(6)  If  0  +  ^10  +  0”  <  020  <  ^ 

Go 

<  0  +  0\Q  <  2n 

02 

(7)  IfO<02Q<0*-  9io  —  oP,  then 

(8)  If  0  +  010  —  <  020  ^  ^  +  ^10  ~ 

oo 

(9)If  0  +  010-0^  <  020  <Tr 

a/ 

Category  II.  If  0^  <  0  <  and  aj  =  oq: 

9  possible  cases  concerning  the  length  of  02. 

0^  <0  +  010  <  a” 

02 

a”  <0  +  010  <  0” 

G2 

(1)  // 0  <  ^20  <  ^  +  010  -  then 

ao 

(A)IfO<02o<0  +  010  -a”,  then 

a/ 

(2)  I  f  0  +  010  —  0^  <  020  ^  ^  +  ^10 

Oj 

(5)  If0  +  010  -a”<02o<0-^  010  -  0^ 

Go 

(3)1  f  0  +  010  +  0”  <02Q<t^ 

ao 

(6)1  f 0  +  010-0”  <02O<ir 

«/ 

II”  <^>  +  ^10  <  2n- 

at 

(7)  IfO<02o<0  +  010  -  </ien 

ao 

(%)If  0  +  010- 0”  <  020  <0  +  0io-ct” 

«/ 

(9)  I  f  0  +  010  —  a”  <  020  <  Tt 

Go 

Category  ID.  If  <0  <  w  and  oi  =  a j.  9  possible  cases  concerning  the  length  of  02- 


o'’  <  0  +  0JO  <  0” 

02 

<0  +  010  < 

'-■K  +  a” 

02 

(1/  ^/O  <  020  ^  0  010  —  oi”i  then 

a/ 

(4)7/0  <^20 

<0  +  010  —  0”i  then 

Go 

(2)  If  0  +  010  -a”  <02Q<0  +  010  -  0” 

GO 

(5)  If  0  +  010 - 

-0”  <02O<0  +  010  -  a” 

«/ 

(3)If  0  +  010-0”  <  «20  <  tt 

a/ 

(6)1  f  0  +  010- 

-  a”  <  02O<  ^ 

Go 

Ol 


V  q’’’  <  &  +  $10  <  2it 

(7)  //O  <  ^20  ^  ^  ®io  —  —  TT,  then 

(8)  If  0  +  ^JO  —  —  TT  <  ^20  <  ^  +  ^10  —  OQ 

(9) Ife  +  9io-0P<e2Q<7t  aj 


We  have  employed  the  definitions  =  ir  —  6^  and  0^  =  z  +  6^.  With  the  foregoing  analysis 
at  hand,  one  can  evaluate  the  integrals  in  Eqs.  (4)-(5)  and  the  external  stress  from  (19).  For 
biaxial  tension,  the  foregoing  integration  domains  remain  the  same  and  the  only  difference  is  the 
evaluaddn  of  6^  in  Eq.  (14). 


II.2.4.  Biaxial  tension/compression  loadings 


In  the  case  of  biaxial  tension/compression  loadings,  some  microcracks  are  open  while  others 
are  closed.  Microcrack  process  and  open  domains  are  coupled  and  affect  each  other.  For  prescribed 
lateral  stress  and  “process  angle”  0^,  one  has  to  compute  the  corresponding  normal  stress 
(evaluated  at  ^  =  6^)  and  the  open/closed  microcrack  angle  boundary'  0"  by  Eqs.  (14)  and  (6)  in 
order  to  predict  progressive  stress-strain  responses.  In  particular,  Eqs.  (14)  and  (6)  can  be  recast  as 


^22  - 


7i  +  V (7i  -  7270)^1?^  72 
72 


As  mentioned  earlier  in  Section  2,  Eq.  (6)  (and  therefore  (26)2]  is  based  on  far-ficld  stresses 
and  IS  therefore  an  approximation  of  the  local  stress  formula  (in  biaxial  tension/compression): 


(p)  =  M„aJ?+M,2ag>0  (27) 

From  a  mathematical  viewpoint,  use  of  (26)  corresponds  to  a  second-order  (in  uj^)  accurate  interac¬ 
tive  microcrack  theory.  On  the  other  hand,  use  of  (26)  1  and  (27)  produces  third-  and  higher-order 
terms  (in  u;^).  Since  our  aim  is  to  develop  a  second-wder  theory  (in  a;^),  Eqs.  (26)  can  be  regarded 
as  a  good  approximation.  Nevenheless,  Eq.  (27)  should  always  be  checked  against  any  possible 
microcrack  orientation  to  ensure  that  a  microcrack  is  indeed  open. 

The  microstructural  integration  domains  can  be  systematically  summarized  as  follows.  Note 
that  there  are  76  possible  cases  in  total  and  a‘  =  z  -  9'  and  0*  =  z +  9'. 


Category  I  ’f  0  <  0  and  a?  =  ar.  19  possible  eases  concerning  the  length  of  a:. 


0  <  ^  +  ^10  < 

02 

<  6  +  6,0  <  6* 

Q2 

(1)7/0  <  Old  <  ^  +  +  then 

«/ 

(5)7/0  <  620  <  6  +  6,0  -  6P,  t/ien 

00 

(2)  I  f  6  6\q  +  <  $20  ^  ^  ^10 

Oo 

(6)  7/  6  +  6,0  -  6P  <  620  <  6  +  6,0  +  6^ 

«/ 

(3)  I  f  6  6*10  —  o*  <  Old  <  ^  +  ^10  “ 

oo 

(7)  7/  6  +  6,0  +  6^  <  620  ^  B  +  6,0  +  6* 

00 

(4)  7/  ^  +  ^10  —  <  ^20  <  Tf 

a/ 

(8)  I  f  B  +  6,0  +  O’*  <  620  ^  Tf 

Oo 

O'*  <.  $  •¥  9\d  ^  C(^ 

02 

<  6  +  6,0  < 

02 

(9)  If  0  <  620  <  ^  +  ^10  -  o',  then 

Od 

(13)  7/0  ^  620  ^  6  +  6,0  —  0^,  then 

a/ 

(10)  I  f  $  +  B\q  —  0“  <  Bid  ^  ^  +  ^10  ~ 

Od 

(14)  I  f  B  ■¥  6,0  —  O'**  <  620  ^  B  6,0  —  0“ 

Oo 

(11)  7/^  +  0\d  —  6^  <  02d  ^  B  +  B]d  +  B’^ 

«/ 

(15)7/6  +  6,0-6*  <62o<6  +  6,o-6p 

Oo 

(12)  I  f  B  +  610  +  B^  <  620  < 

00 

(16)7/6  +  6,0  -  6”  <  620  <  r 

a/ 

«*  <  fl  +  010  <  o* 

02 

(17)  I f  B  +  B\d  —  6"  <  620  ^  B  +  Bid  ~  B^ 

00 

(18)  7  f  B  +  Bid  —  B^  <  620  ^  B  -i-  Bid  +  B^‘ 

(19)  1  f  B  ■¥  610  +  6^  <  620  ^  B  +  ^10  +  B' 

Oo 

42 


Category  11.  If  0  <  0*  and  ai  =  qq:  19  possible  cases  concerning  the  length  of  02 


9^  <e-^  Bio  <  B' 

<12 

q’  <  0  +  0,0  < 

02 

(1)  I f  0  <  620  <  B  B\o  —  B^,  then 

oq 

(5)  //O  <  020  $  0  +  0,0  ~  f^en 

no 

(2)  I  f  B  •¥  B\()  —  B^  <.  620  ^  9  +  6\o  +  B^ 

(6)  //  0  +  0,0  —  0*  <  020  ^  0  010  “  0'’ 

no 

(3)  I f  B  +  B\o  +  0’’  <  B20  9  +  9io  +  B" 

no 

(7)  7/0 +  0,0-  0'’  <  020  <0  +  010  +  0'’ 

«/ 

(4)  I  f  9  ■¥  ^10  +  Q*  <  020  ^  ^ 

no 

(8)  7/  0  +  0,0  +  0'’  <  020  < 

no 

qP  <0  +  0,0  < 

02 

<  0  +  010  <  2ff 

02 

(9)  I f  0  <  020  ^  0  +  010  ~  0^,  then 

a/ 

(13)  7/0  <  020  <  0  +  010  -  tften 

no 

(10)  I  f  B  +  0,0  ~  Ct^  <  020  ^  0  +  010  ~  0* 

no 

(14)7/0  +  0,0-  /r  <  020  <0  +  010-0'’ 

«/ 

(ll)//0  +  0,O~0*  <  020  ^  0  +  0,0  0^ 

no 

(15)  7/  0  +  0,0  —  Ct^  <  020  <  0  +  010  —  n* 

no 

(12)  1  f  B  +  0,0  —  0^  <  020  ^  ^ 

«/ 

(16)  I  f  B  +  0,0  —  Q*  <  020  <  ^ 

no 

0*  <  0  +  0,0  ^  0* 

02 

(17)  I f  B  +  0,0  —  0*  ^  020  ^  0  +  010  —  0^ 

no 

(18)  1  f  B  +  0,0  —  0^  <  020  0  +  010  0^ 

«/ 

(19)  I  f  B  0,0  +  0^  <  020  <  0  +  0,0  +  0* 

no 

Category  ID.  If  q*  <  0  <  and  gi  =  oq:  19  possible  cases  concerning  the  length  of  02. 


q’  <e  +  exo<  ap 

02 

Qp  <  0  +  010  < 

02 

(1)  If  0  <  ^20  ^  ^  +  ^10  “  0*,  then 

ao 

(5)  7/0  <  020  ^  0  +  010  ~  Q**!  then 

(2)  If  B  +  Bio  —  0"  <  $20  ^  $  +  $10  —  $^ 

Qo 

(6)  7/  0  +  0,0  —  <  $20  ^  0  0io  ~  0* 

00 

(3)  If  $  +  $xo  —  B^  <.  $20  ^  ^  +  ^10 

«/ 

(7)  I  f  B  +  0,0  —  0*  <  020  ^  0  +  010  ~  0^ 

Oo 

(4)  I f  $  +  $xo  +  <  020  ^  ^ 

ao 

(8)  I  f  $  +  0,0  —  0^  <  020  ^ 

^  <B-^0xo< 

02 

TT  +  0*  <  0  +  0,0  ^  27r 

02 

i9)IfO<B2Q<$  +  $10  -  t/icn 

no 

(13)7/0  <  020  <  0  + 010  —  0*  —  ’T,  then 

Oo 

(10)  //  $  +  010  —  0^  <  $20  ^  0  ^10  ~  0^ 

«/ 

(14)  7/  0  +  0,0  —  0"  <  020  ^  0  +  010  ~ 

00 

(11)//  0  +  010  —  Qr^  <  020  ^  $  "^  $10  ~  0* 

Oo 

(15)  7/ 0  +  0,0  —  0^  <  $20  ^  0  010  ~ 

0/ 

(12)  I  f  $  +  010  —  0*  <  020  <  TT 

ao 

(16)  7/  0  +  0,0  —  <  020  ^ 

Oo 

0'  <  0  +  010  <  TT  +  a* 

02 

'mif  $+$10- 0'  <$2o<e+6io- 0^ 

oo 

(18)  I  f  $  +  010  —  0^  <  $20  <  0  +  010  —  0^ 

(19)7/0  +  0,0 -aP  <  020  <  0  +  0,0 -a* 

oo 

02 


Category  IV.  If  <  ^  <  :r  and  at  =  ay :  19  possible  cases  concerning  the  length  of  ai- 
<  e  +  Oio  <  0^  02  0^  <0-i-  010  <  0‘ 


(1)  7/0  <  ^20  ^  ^  +  ^10  ~  0^,  then 

«/ 

(5)  7/0  <  ^20  <  ^  +  0\Q  -  then 

ao 

(2)  I  f  6  ■¥■  6\o  —  <  $20  ^  +  ^10  ~ 

ao 

(6)  7/  0  +  ^jo  —  /?’’<  ^20  ^  0  010  ~  o’* 

0/ 

(3)  I f  6  +  6\q  —  6’  <  620  ^  0  010  —  0^ 

ao 

(7)  1  f  0  ■¥  0\o  ~  o’"  <  020  ^  0  010  ~  0* 

Oo 

(4)  I  f  6  +  $10  —  0^  <.  020  ^  rc 

(8)  I  f  0  ■¥  0\o  —  0“  <  020  ^  Tf 

Oq 

TT  +  Q*  <  0  +  ^10  <  TT  + 

02 

JT  +  o’*  <  ^  +  ^10  <  27r 

02 

(9)  //  0  <  ^20  <  ^  +  ^10  —  ck*  —  TT,  then 

Oo 

(13)  I f  0  <  02D  <  0  +  010  —  then 

0/ 

(10)  1  f  ^  +  0\Q  —  0*  <  020  ^  0  +  0\o  —  0^ 

Oo 

(14)  1  f  0  +  0\Q  —  o’*  —  TT  <  020  $  ^  +  010  —  q"  —  If 

ao 

(11)7/  ^  +  010  —  0^  <  020  ^  ^  +  ^10  ~  0^ 

a/ 

(15)  1  f  0  ■¥  0)0  —  0*  <  020  ^  0  ®10  —  0’’ 

Oo 

(12)  I  f  0  +  010  —  0^  <  020  ^  ^ 

Oo 

(16)  J  f  0  0)0  —/?’’<  020  ^  ^ 

0/ 

0“  <  0  +  010  ^  TT  +  O'* 

02 

(17)  I  f  0  +  0JO  —  0‘  ^  020  ^  ^  +  ^10  ~  0^ 

Oq 

(18)  I  f  0  +  0\o  —  0^  <■  02D  ^  +  ^10  ~  o’’ 

(19)  I  f  0  +  0\o  —  Ot^  <  020  ^  ^  +  010  ~  0* 

Oo 

The  foregoing  integration  domains  completely  define  the  integrals  in  Eqs.  (3H5),  and  there¬ 
fore  effective  compliances  can  be  computed.  In  a  strain-driven  algorithm,  the  axial  tensile  stress 
can  also  be  computed  accordingly.  It  is  emphasized  that  microcrack  interactions  affect  the 
values  of  na  and  thus  ag;  see  (26)i.  In  turn,  ag  affects  the  “opening  angle”  0’\  see  (26)2- 

11.3.  Numerical  algorithms  and  applications 

II.3.1.  Numerical  integration  algorithms 

As  previously  mentioned,  one  needs  to  numerically  evaluate  the  integrals  involved  in  Eqs.  (3)- 
(5).  Furthermore,  one  needs  to  develop  efficient  numerical  algorithm  to  solve  Eq.  (26)  involving 
both  process  and  opening  domains.  First,  the  numerical  integrations  of  (3)-{5)  can  be  efficiently 
achieved  by  using  the  “Gauss  quadrature”  scheme.  In  carrying  out  the  Gauss  quadrature  numerical 
scheme,  we  have  chosen  the  following  integration  ranges:  I2a,  20a]  for  the  distance  r,  and  [0,  «■]  for 
the  angle  0iq.  The  integration  bounds  for  ^20  and  6  arc  subsequently  determined  from  the  systematic 
analyses  presented  in  Sections  2.Z  and  2.4. 

Another  issue  is  the  number  of  Gauss  points  needed  for  each  integration  zone  aforementioned 
in  Sections  2.3  and  2.4.  Table  1  gives  the  convergence  behavior  of  Gauss  quadrature  with  respect 
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to  the  number  of  Gauss  points  used  for  each  integration  zone  for  the  "aligned  microcracks”  problem 
previously  considered  in  Section  3.2  of  Part  I.  Clearly,  use  of  20  Gauss  points  in  each  integration 
zone  is  sufficiently  accurate.  An  efficient,  iterative  numerical  integration  algorithm  to  solve  Eq.  (26) 
is  summarized  in  Box  1. 

Remark  3.1.  For  uniaxial  and  biaxial  tensile  loadings,  Eq.  (6)  is  automatically  satisfied  un¬ 
less  the  axial  stress  goes  to  infinity;  i.e.,  6'  =  7r/2.  Therefore,  the  iterative  process  for  finding  the 
opening  angle  bound  in  Step  (S)  of  Box  1  is  unnecessary.  | 

II3.2.  Uniaxial  tension  loading  tests 
% 

In  this  section,  we  consider  two  (plane  strain)  uniaxial  tension  tests  with  different  initial  mi¬ 
crocrack  ctmeentrations:  uq  =  naQ  =  0.1  (moderate  density)  and  0.4  (high  density).  It  should  be 
realized  that  these  two  concentrations  are  equal  to  O.Itt  and  0.47r  in  terms  of  the  commonly  used 
definition  (w  =  mral).  Moreover,  the  k  value  in  Eq.  (11)  is  taken  as  0.5,  initial  Poisson’s  ratio 
=  0.2,  and  po  =  0.6  in  (7).  For  convenience,  all  strains  are  normalized  (dimensionless)  as  follows: 


The  normalized  stress-strain  curves  are  displayed  in  Figure  3.  All  co.mpliances  have  also  been 
normalized  by  dividing  the  common  factor  (1  -  v^)IE.  We  observe  that  microcrack  interactions 
effectively  lower  the  onset  (threshold  stress)  of  microcrack  propagation.  The  existence  of  “comer 
points”  on  the  axial  stress-strain  curves  in  Fig.  3  is  a  consequence  of  the  deterministic  fracture 
criterion  (14)  and  the  uniform  initial  microcrack  size  oq.  If  we  employ  a  probabilistic  fracture 
criterion  or  nonuniform  initial  microcrack  sizes,  the  “comer  points”  will  disappear  and  be  replaced 
by  smooth  transitions  [c.f.  Ju  and  Lee  (1991),  and  Lee  and  Ju  (1991)].  It  is  noted  that  the  axial 
stress/lateral  strain  curves  are  the  same  for  wo  =  0. 1  and  0.4  for  Taylor’s  model  due  to  the  fact 
that  microcracks  have  no  contribution  to  cn-  In  the  case  of  the  proposed  microcrack  interaction 
theory,  a  significant  reduction  in  the  magnitude  of  lateral  strains  is  observed.  That  is,  microcrack 
interactions  enhance  Mode  I  displacements  and  reduce  Mode  D  displacements,  and  therefore  reduce 
the  magnitude  of  lateral  contraction. 

Figures  4  and  5  show  the  normalized  overall  compliance  components  predicted  by  both  Tay¬ 
lor’s  and  the  proposed  models.  It  is  noted  that  Taylor’s  model  does  not  change  either  {Sn)  or  (Sji) 
components  despite  damage.  By  contrast,  the  proposed  interaction  model  shows  that  the  shear  com¬ 
ponents  (S12)  and  {S21)  change.  It  is  also  observed  that  (S12)  and  (S2))  become  slightly  different 


as  the  microcrack  concentration  increases.  As  the  normalized  axial  stress  goes  to  infinity,  one 
may  expect  that  the  overall  compliance  matrix  asymptotically  becomes  isotropic  because  all  pre¬ 
existing  random  microcracks  grow  to  same  length  a j.  This  is  indeed  the  case  for  Taylor’s  model 
and  the  self-consistent  model.  However,  this  may  not  be  always  true  if  microcrack  interactions  arc 
considered.  The  reason  is  that  some  initially  open  stationary  microcracks  (ao)  in  the  very  close 
neighborhood  of  ^  =  7r/2  may  become  closed  due  to  interactions  with  neighboring  propagated  mi¬ 
crocracks  (a/).  Biaxial  tension  loadings  can  be  treated  in  a  similar  fashion  and  therefore  are  not 
repeated  here. 

II3.3.  ^iaxial  tension/compression  loading  tests 

One  has  to  determine  the  open  domain  (0*)  and  the  corresponding  axial  tensile  stress  by 
solving  the  system  of  equations  (26)  for  specified  values  of  6^  and  aff  in  order  to  compute  overall 
compliances.  A  good  initial  guess  for  6’  can  be  obtained  from  Taylor’s  model  by  solving  Eqs.  (6) 
and  (19).  With  this  initial  guess  one  can  obtain  the  trial  value  by  using  Eq.  (26)i. 
Subsequently,  an  improved  trial  value  6’^^^  can  be  rendered  by  using  (26)2,  and  so  on.  In  what 
follows,  we  adopt  the  following  propenies:  i/  =  0.2,  po  =  0.6  and  k  =  0.5. 

To  illustrate  the  performance  of  the  proposed  numerical  algorithm  summarized  in  Box  1,  let  us 
consider  the  onset  point  of  the  “cleavage  1  ”  microcrack  growth.  Table  2  shows  the  convergence 
behavior  of  this  case  for  u;o  =  0.1  and  afj’  =  -0.1.  It  is  observed  that  remarkable  convergence  is 
reached  within  only  two  to  three  iterations. 

The  normalized  axial  stress-strain  curves  for  four  different  biaxial  tension/compression  tests 
with  Wo  0. 1, 0.4  and  lateral  compression  afj  =  -0.1,  -0.4,  respectively,  are  similar  to  those  of 
the  uniaxial  tension  tests  in  Figure  3.  However,  the  behavior  of  axial  stress/  lateral  strain  curves 
depends  on  (or  the  ratio  |)  and  wq;  see  Figures  6  and  7.  The  magnitude  of  lateral 

(contractive)  strains  arc  higher  for  higher  lateral  compression  =  —0.4;  compare  Figures  6 
and  7.  On  the  other  hand,  the  loading  ratio  | ^22/^1?  I  affects  the  size  of  the  open  and  process 
domains.  TTic  sizes  and  orientations  of  these  open  and  process  domains,  in  turn,  affect  the  effects 
of  microcrack  interaction  on  lateral  strains.  In  particular,  for  higher  loading  ratio  1.  the 

net  effects  of  microcrack  interaction  enhance  the  Mode  I  (dilatational)  contributions  and  reduce 
the  Mode  11  (contractive)  contributions.  By  contrast,  when  the  loading  ratio  I  is  relatively 

small,  the  net  effects  of  microcrack  interaction  enhance  the  Mode  II  (contractive)  contributions; 
sec  Fig.  7.  It  is  worth  mentioning  that  if  ag  continues  to  increase  in  Fig.  7,  then  the  two  solid  lines 
in  Fig.  7  will  also  intersect  each  other  just  as  they  do  in  Fig.  6. 
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Finally,  Figure  8  displays  the  normalized  {S22)  and  (Sn)  compliances  vs.  axial  stresses  for 
ujQ  =  0.4  and  =  —0.1.  Due  to  closure  of  some  microcracks,  overall  compliances  are  always 
anisotropic.  It  is  also  noted  that  due  to  microcrack  interactions,  the  overall  compliance  matrix 
becomes  non-symmetric  (i.e.,  Sn  4  52i). 


n.4.  A  higher-order  microcrack  interaction  model 

The  statistical  micromechanical  damage  models  presented  so  far  are  based  on  the  concept  of 
pairwise  microcrack  interaction.  This  pairwise  microcrack  interaction  mechanism  essentially  cor¬ 
responds  to  a  second-order  damage  theory  (in  w^).  Within  the  context  of  the  ensemble-volume 
average  approach,  one  can  systematically  incorporate  many-microcrack  interaction  mechanisms 
into  the  proposed  framework.  In  essence,  one  needs  to  re-derive  the  ensemble-average  of  the  per¬ 
turbation  in  local  stress  field  due  to  n-microcrack  interactions  (n  >  3). 

To  illustrate  the  foregoing  statements,  let  us  consider  a  threc-microcrack  (third-order  in  u;^) 
\ 

interaction  mechanism.  Following  the  definitions  and  assumptitms  described  in  Part  I  of  this  study, 
one  can  recast  the  local  ensemble  stress  perturbation  in  Eq.  (12)  of  Part  I  as  follows: 

{f)+{t)  =  ({|})+({|})  (29) 

where  (T)  is  the  first-order  local  ensemble  stress  perturbation  due  to  pairwise  microcrack  interac¬ 
tions,  and  (T)  is  the  second-order  local  ensemble  stress  perturbation  due  to  the  third  wder  micro¬ 
crack  interactions.  In  particular,  (f )  can  be  expressed  as  [cf.  Eq.  (13)  of  Part  I] 

(T)  *  JjTjix,  a,  &; xj ,au9i  jx2, 02,  W(^2,  «2,  ^zjx,  a,  B; Xi,  01,61)  d02da2dx2  (30) 

Here,  (T)(x,  a,  B\  Xi ,  aj ,  9\  |x2, 02,  ^2)  is  the  second-order  stress  perturbation  of  a  microcrack  cen¬ 
tered  at  X  with  (a,  6),  given  a  microcrack  centered  at  xi  with  (oi ,  61),  averaged  over  the  subclass  of 
realizations  having  a  microcrack  centered  at  X2  with  (02, 02)-  Further,  /(X2, 02,  ^2!*,  a,  Xj ,  a, ,  ^i) 
is  the  conditional  probability  function  for  finding  a  microcrack  centered  at  X2  with  (02, 62)  given 
two  microcracks  fixed  at  x  with  (a,  6)  and  at  Xj  with  (ai,5i).  /(X2,a2,^2ix,a,^;xj,ai,^i)  can  be 
simplified  to  /(x,  a,  6)  by  the  assumptions  of  local  homogeneity  and  reasonable  randomness. 

The  solutions  of  T  for  a  system  of  three  (or  many)  arbitrarily  located  and  oriented  microcracks 
were  previously  investigated  by  Horii  and  Nemat-Nasser  (1985),  and  Kachanov  (1987).  In  partic¬ 
ular,  Kachanov  (1987)  pursued  extensive  numerical  computations  in  order  to  obtain  the  “transmis¬ 
sion  factors”  for  local  stresses.  Alternatively,  one  can  derive  approximate  analytical  soluticms  for 
T  by  following  the  procedures  presented  in  Sec.  2.3  of  Part  I. 

For  clarity,  let  us  express  T  and  (e*)(x)  as  follows  [cf.  Eqs.  (25)  and  (28)  of  Part  I]: 

(T)  =  (T-  -^f -f-t)  =  {Ko  +  mfia,e){K)  +  f\x)fHa,e){K'))  •  (31) 
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(e')(x)  =  { (S*’)  +  (S*2)  -K  (S-3)}  •  =  (S')  ■ 

(32) 

and  [see  Eqs.  (29>-(30)  of  Part  I] 

(S*'}  =  J^g.Kofia,e)deda 

(33) 

(S-2)  =  ja^Je  (K)/^o, « dfda 

(34) 

(S*^)  =  liL^/3(x)  g  •  (K')/'(a,  0)  d$da 

(35) 

In  the  above  equations,  K,  Ko  and  (K)  have  been  previously  defined  in  Eqs.  (24),  (26)  and  (27)  of 
Part  I,  respectively.  Similarly,  K'  (or  T)  can  also  be  constructed  as  follows.  One  starts  by  expanding 
Eq.  (15)  of  Part  I  into  six  linear  equations  with  j  ~  1,2, 3.  Then,  one  obtains  Eqs.  (16)-(17)  of 
Pan  I  with  the  understanding  that  j  =  1,2,3  and  the  permutations  1-2,  2-3,  S-1  are  involved  in 
Eq.  (17).  Eq.  ( 1 8)  of  Pan  I  is  modified  to  include  the  Jtth  {k  ^  j  and  k  ^  i)  microcrack’s  contribution 
to  stress  perturbations.  Subsequently,  Eq.  (19)  of  Pan  I  is  expanded  to  a  six-cquadon  system  with 
a  denoting  a  six  by  six  coefficient  matrix.  The  components  of  a  are  identical  to  those  in  Eq.  (20) 
of  Part  I  with  the  understanding  that  the  1-2-3  permutations  are  involved.  Therefore,  we  arrive 
at  explicit  formulas  similar  to  Eqs.  (21)-(24)  in  Pan  I.  Finally,  it  is  noted  that  the  computation  of 
(K')  involves  integration  over  the  probabilistic  domain  of  all  possible  positions  and  orientations 
of  two  active  neighboring  microcracks.  Following  the  standard  procedures  presented  in  Sec.  2.4 
of  Part  I,  it  can  be  shown  that  (S*^)  in  Eq.  (35)  introduces  third-order  terms  (in  u?^)  to  overall 
compliance.  Therefore,  a  third-order  statistical  micromechanical  model  can  be  constructed.  By 
repeating  the  foregoing  procedures,  we  can  formulate  a  complete  hierarchical  family  of  statistical 
micTomechanical  theories  of  arbitrary  order. 
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UJ.  Conclusions 


In  this  paper,  basic  two-dimensional  formulations  presented  in  Part  I  of  this  investigation  are 
further  extended  to  account  for  “cleavage  1”  microcrack  growth  proces^s  with  microcrack  inter¬ 
actions.  Stemming  from  a  micromechanical  viewpoint,  a  mixed  fracture  criterion  is  used  to  derive 
microcrack  kinetic  equations  for  interacting,  randomly  located  and  oriented  microcracks.  As  a  re¬ 
sult,  the  proposed  approach  is  capable  of  predicting  the  progressive  stress-strain  curves  for  biaxial 
loadings.  This  constitutive  predictive  capability  is  not  found  in  the  current  literature  on  stationary 
microcrack  interaction  models.  Both  tensile  and  mixed  tensile/compressive  loadings  are  considered 
in  Sec.  2,  including  systematic  analyses  regarding  detailed  integration  domains.  Further,  efficient 
numerical  integration  algorithms  are  presented  for  biaxial  loadings  in  Sec.  3.  A  number  of  model 
predictions  and  comparisons  are  also  performed.  It  is  shown  that  load-induced  anisotropy  and  non- 
symmetry  of  overall  compliances  can  be  predicted  by  the  proposed  micromechanical  models.  It  is 
emphasized  that  the  proposed  models  do  not  employ  any  fitted  “material  parameters”  and  do  not 
require  the  use  of  Monte  Carlo  simulations. 

The  proposed  pairwise  interaction  framework  is  subsequently  generalized  to  a  higher-order 
damage  model  to  account  for  higher  (third  and  above)  order  microcrack  interaction  effects  in  Sec.  4. 
However,  a  higher-order  micromechanical  damage  model  requires  much  more  complicated  anal¬ 
yses.  Experimental  validation  of  the  proposed  micromechanical  evolutionary  damage  models  is 
needed  to  further  assess  their  validity  for  practical  engineering  problems. 
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Box  1.  Numerical  int«*ffration  algorithm 


(1)  Specify  and  af] ,  and  guess  an  initial  bound  of  opening  domain  6*  (e.g.,  from  Taylor’s 
model  analysis). 


(2)  Integrate  (3)-(5)  over  6  and  6io.  Note  that 


% 


fa/,  if  d  in  process  domains; 
\  oq,  otherwise. 


(3)  Integrate  (3)-(5)  over  d^a  for  each  value  of  ^  +  ^lo.  Note  that 


fa/,  if  ^20  in  process  domains; 
I  ao,  otherwise. 


(4)  Integrate  Eqs.  (4)-(5)  from  r  =  (oi  +  02)  to  max{lQa\ ,  20a2). 

(5)  For  a  given  fracture  criterion,  find  the  corresponding  threshold  stress  (26) j.  Cal¬ 

culate  the  updated  opening  angle  bound  9*  by  (26)2  or  (27).  If  the  difference  between 
the  initial  guess  and  the  updated  ^’-valuc  is  not  acceptable,  then  go  back  to  Step  (2)  and 
re-iterate  until  proper  convergence  is  reached. 

(6)  Go  to  Step  (1)  and  specify  new  9^  and  values. 
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Table  1.  Convergence  behavior  (r^,,  = 


No.  Gauss  points 

io 

20 

30 

40 

80 


0.5862290322 

0.5864980549 

0.5864981079 

0.5864981080 

0.5864981080 


{^23) 

0.2405392992 

0.2406183532 

0.2406183630 

0.2406183630 

0.2406183630 


Table  2.  Convergence  behavior 

No.  of  iteration  (0 

0 

1.2645189576 

- 

1 

1.2592833120 

5.235645617  lE-03 

2 

1.2592966258 

1.3313781921E-05 

3 

1.2592965919 

3.3915545705E-08 

4 

1.2592965920 

8.63957794 19E- 11 

5 

1.2592965920 

2.196021 1427E- 13 

0.7.  Figure  captions 


Figure  1.  Schematic  sketch  for  the  definitions  of  process  domains  (0’’)  and  opening  domains 
(&•). 

Figure  2.  Schematic  sketch  for  various  integration  domains,  with  x  denoting  the  global 
reference  axis,  (x, ,  y,)  denoting  the  first  set  of  local  coordinates,  and  (x^,  y2)  denoting  the  second 
set  of  local  coordinates. 

Figure  3.  The  nonnalized  axial  stress-strain  curves  for  uniaxial  tension  tests.  The  dotted 
lines  and  solid  lines  are  predictions  from  Taylor’s  and  the  proposed  models,  respectively. 

Figure  4.  The  normalized  (522)  and  (Sjj)  compliances  vs.  axial  stress  curves  for  uniaxial 
tension  tests  with  wo  =  0. 1.  The  compliances  have  been  normalized  by  the  factor  (1  —  v^)IE.  The 
dotted  lines  and  solid  lines  are  predictions  from  Taylor’s  and  the  proposed  models,  respectively. 

Figure  5.  The  normalized  (S12}  and  (S21)  compliances  vs.  axial  stress  curves  for  uniaxial 
tension  tests  with  wo  =  0.1  and  wo  =  0.4.  The  dotted  line  is  the  prediction  from  Taylor’s  model 
for  both  Wo  =  0.1  and  wo  =  0.4.  The  dash-dot  lines  and  solid  lines  are  {S12)  and  (S21)  predictions, 
respectively,  from  the  proposed  model. 

Figure  6.  The  normalized  axial  stress  vs.  lateral  strain  curves  for  biaxial  tension/compression 
tests  with  5-^  =  -0.1.  The  dotted  lines  and  solid  lines  are  predictions  from  Taylor’s  ami  the 
proposed  models,  respectively. 

Figure  7.  The  normalized  axial  stress  vs.  lateral  strain  curves  for  biaxial  tension/compression 
tests  with  afi  =  -0.4.  The  dotted  lines  and  solid  lines  are  predictions  from  Taylor’s  and  the 
proposed  models,  respectively. 

Figure  8.  The  normalized  {522)  and  (5ii)  compliances  vs.  axial  stress  curves  for  biaxial 
tension/compression  tests  with  w©  =  0.4  and  =  -0.1.  The  dotted  lines  and  solid  lines  are 
predictions  from  Taylor’s  and  the  proposed  models,  respectively. 
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Figure  2.  Schematic  sketch  for  various  integration  domains,  with  x  denoting  the  global 
reference  axis,  (Xj,  yj)  denoting  the  first  set  of  local  coordinates,  and  (Xj,  y^)  denodng  the  second 
set  of  local  coordinates. 


normalized  {Sn)  and  (Sn)  compliances  vs.  axial  stress  curves  for  uniaxial 
=  0.1.  The  compliances  have  been  normalized  by  the  factor  (1  -  i/*)/E.  The 
litjcs  arc  predictions  from  Taylor's  and  the  proposed  models,  respectively. 


Figure  6.  The  normalized  axial  stress  vs.  lateral  strain  curves  for  biaxial  tension/compression 
tests  with  =  —0.1.  The  dotted  lines  and  solid  lines  are  predictions  from  Taylor’s  and  the  pro¬ 
posed  models,  respectively. 


Figure  8.  The  noimalized  {S^i)  and  {S\\)  compliances  vs.  axial  stress  curves  for  biaxial 
tension/compression  tests  with  wq  =  0.4  and  =  —0.1.  The  dotted  lines  and  solid  lines  arc 
predictions  from  Taylor’s  and  the  proposed  models,  respectively. 


PART  III 

A  Three-Dimensional  Statistical  Micromechanical 
Theory  for  Brittle  Solids  with  Interacting  Microcracks 


# 


nLO.  Abstract 

A  three-dimensional  statistical  micromechanical  theory  is  presented  to  investigate  effective 
elastic  moduli  of  brittle  solids  with  many  randomly  located,  penny-shaped  microcracks.  The  macro¬ 
scopic  constitutive  relations  are  statistically  and  micromechanically  derived  by  taking  the  ensemble 
average  over  all  possible  realizations  which  feature  the  same  statistical  distribution  of  microcracks. 
Approximate  analytical  solutions  of  a  two-microcrack  interaction  model  are  presented  to  account 
for  pairwise  microcrack  interaction  among  many  randomly  located,  aligned  microcracks.  There¬ 
fore,  the  ensemble-averaged  stress  penurbations  due  to  microcrack  interaction  can  be  constructed 
in  closed-form.  The  overall  effective  compliances  of  microcrack-weakened  brittle  solids  arc  de¬ 
rived  by  further  taking  the  volume  average  of  the  ensemble-averaged  stress-strain  relations  over  the 
entire  mesostructural  domain  of  a  representative  volume  element.  Some  numerical  examples  arc 
given  to  illustrate  the  behavior  of  the  proposed  method.  Comparisons  with  some  existing  methods 
are  also  appended.  Finally,  a  higher-order  ensemble-average  formulation  of  microcrack  interaction 
is  briefly  discussed.  The  proposed  framework  is  fundamentally  different  from  existing  “effective 
medium  methods”  which  do  not  depend  on  microcrack  locations  and  conflgurations.  It  is  empha¬ 
sized  that  no  Monte  Carlo  simulations  are  necessary  in  the  proposed  framework. 
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DLL  Introduction 


Damage  mechanics  has  emerged  as  an  important  branch  of  solid  mechanics  pertaining  to  mi¬ 
crocrack  nucleation,  gro^vih,  coalescence  and  interaction,  as  well  as  their  effects  on  overall  mechan¬ 
ical  responses  of  materials.  Basically,  there  are  two  types  of  models:  (a)  phenomenological  damage 
models  (e.g.,  Ju  (1989a, b),  Chow  and  Wang  (1987)),  and  (b)  micromechanical  damage  mo  lds.  We 
refer  to  Krajcinovic  (1989)  for  an  excellent  comprehensive  literature  review  on  damage  mechanics. 

The  majority  of  micromechanical  damage  theories  falls  into  the  category  of  “effective  medium 
methods”.  When  microcrack  interactions  are  totally  ignored,  the  resulting  micromechanical  mod¬ 
els  are  termed  “Taylor’s  models”;  see,  e.g.,  Krajcinovic  and  FancUa  (1986),  FaiwUa  and  Krajci¬ 
novic  (1988),  and  Ju  (1991b)  for  some  micromechanical  evolutionary  models.  The  “self-c(Misistent 
method”  (Hill  (1965))  was  applied  to  damaged  solids  by  Budiansky  and  O’Connell  (1976)  focus¬ 
ing  on  randomly  distributed,  weakly  interacting  microcracks.  The  self-consistent  method  was  fur¬ 
ther  explored  by  Horii  and  Nemat-Nasser  (1983),  Laws  et  al.  (1983),  Laws  and  Dvorak  (1987), 
Sumarac  and  Krajcinovic  (1987, 1989),  Krajcinovic  and  Sumarac  (1989),  Horii  and  Sahasakmon- 
tri  (1990),  Nemat-Nasser  and  Hori  (1990),  Ju  (1991a),  Ju  and  Lee  (1991),  and  Lee  and  Ju  (1991), 
etc..  In  addition,  a  three-phase  “generalized  self-consistent  model”  was  proposed  by  Cliristensen 
and  Lo  (1979),  and  Christensen  (1990).  The  “differential  scheme”  was  investigated  by  Roscoe 
(1952,  1973),  McLaughlin  (1977),  Laws  and  Dvorak  (1987),  Hashin  (1988),  Nemat-Nasser  and 
Hori  (1990),  and  so  on. 

On  the  other  hand,  based  on  variational  principles,  Hashin  and  Shtrikman  (1962,  1963)  pro¬ 
posed  upper  and  lower  bounds  for  elastic  composites  with  elastic  inclusions.  Extension  of  their 
work  to  penny-shaped  microcracks  is  possible,  see,  for  example,  Willis  (1977).  Moreover,  the 
“Mori-Tanaka  method”  was  studied  by  Mori  and  Tanaka  (1973),  Benveniste  (1986),  Zhao,  Tandon 
and  Weng  (1989),  and  Weng  (1990).  The  Mori-Tanaka  method  produces  identical  results  to  the 
Hashin-Shtrikman  bounds  (1962, 1963)  under  many  situations;  see  Weng  (1990)  for  details. 

It  is  noted  that  effective  medium  methods  are  only  valid  for  low  or  at  most  moderate  micro- 
crack  concentrations  since  they  do  not  depend  on  microcrack  locations  or  configurations.  When 
microcrack  concentrations  are  higher  and  microcrack  spacings  are  smaller,  strong  microcrack  in¬ 
teractions  become  significant  and  microcrack  locations  should  be  accounted  for.  Excellent  strong 
micTocrack  interaction  models  were  proposed  by  Gross  (1982),  Horii  and  Nemat-Nasser  (1985), 
and  Hori  and  Nemat-Nasser  (1987)  for  two-dimensional  deterministic  microcracks  (not  at  the  cot- 
stitutive  level).  More  comprehensive  studies  were  investigated  by  Kachanov  (1985),  Kachanov 


and  Montagut  (1986),  Kachanov  (1987),  Montagut  and  Kachanov  (1988),  Kachanov  and  Laures 
(1989),  and  Laures  and  Kachanov  (1991)  for  two-  and  three-dimensional  deterministic  arbitrary  mi¬ 
crocrack  arrays  (at  the  constitutive  level).  In  particular,  the  valuable  work  due  to  Kachanov  (1987) 
and  Kachanov  and  Laures  (1989)  are  ideally  suited  for  deterministic  microcracks.  However,  one 
needs  to  perform  Monte  Clarlo  simulations  associated  with  their  method  in  the  event  of  statistically 
distributed  miciocrack  arrays. 

In  this  paper,  we  attempt  to  establish  a  statistical  micromechanical  framework  to  predict  effec¬ 
tive  moduli  of  brittle  solids  with  many  interacting,  randomly  located  penny-shaped  microcracks. 
Randonf  microcrack  locations  and  relative  configurations  are  accommodated  through  probability 
density  functions  and  the  ensemble  average  approach.  It  is  noteworthy  that  Batchelor  (1970), 
Batchelor  and  Green  (1972),  Willis  and  Acton  (1976),  Hinch  (1977),  and  Chen  and  Acrivos  (1978a,b) 
employed  similar  methods  to  characterize  elastic  particulate  composites.  On  the  other  hand,  the 
proposed  method  is  quite  different  from  that  proposed  by  Hudson  (1980,  1981,  1986).  Hudson's 
method,  though  using  the  ensemble  average  approach,  is  based  on  a  second  ^rder  stiffness  tlKory 
and  therefore  leads  to  irrational  behavior  for  solids  with  moderate  or  high  miciocrack  concentra¬ 
tions  (Sayers  and  Kachanov,  1991).  The  proposed  approach  is  free  from  this  anomaly. 

This  paper  is  organized  as  follows.  In  Section  2,  the  ensemble  average  approach  to  derive  local 
damaged  stress-strain  relations  is  introduced.  Approximate  closed-form  solutions  are  subsequently 
presented  for  the  interaction  problem  of  two  aligned  but  arbitrarily  located  penny-shaped  micro- 
cracks.  The  overall  moduli  of  a  representative  volume  element  are  then  derived  by  the  volume 
averaging  process  in  Section  3.  In  Section  4,  applications  are  made  to  some  special  cases.  In 
particular,  in  the  event  of  dilute  microcracks,  the  present  approach  recovers  the  well-known  Taylor’s 
model  by  neglecting  interactions  among  microcracks.  Comparisons  with  some  existing  methods 
are  also  presented.  Finally,  a  higher-order  microcrack  interaction  formulation  within  the  proposed 
framework  is  discussed  in  Section  5. 
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IIL2.  An  ensemble  average  appruach  to  3-D  aligned  microcrack  interaction 

In  this  section,  we  start  by  studying  the  ensemble  averages  of  microcrack-perturbcd  stresses 
and  strains.  Approximate  closed-form  analytical  solutions  of  two  arbitrarily  located,  aligned  mi¬ 
crocracks  are  subsequently  given.  Finally,  some  test  problems  are  examined  to  assess  the  accuracy 
of  the  proposed  approximate  analytical  solutions. 

Following  the  cunent  literature,  it  is  assumed  that  the  volume-average  stress  a  is  approximately 

equal  to  the  far-field  stress  <t°°.  This  common  assumption  can  actually  be  removed  and  further 

improvements  are  possible.  These  and  related  issues  will  be  addressed  in  a  forthcoming  paper. 

% 

in.2.1.  Ensemble  averages  of  microcrack-perturbed  stresses  and  strains 

Due  to  the  existence  and  interaction  of  microcracks,  local  stresses  and  strains  in  the  matrix 
material  are  perturbed.  The  ensemble  average  approach  hinges  on  the  concept  that  local  stresses, 
strains  and  compliances  (or  stiffnesses)  at  a  typical  point  within  a  representative  volume  element 
(RVE)  of  a  microcrack-weakened  solid  can  be  obtained  by  averaging  over  the  ensemble  of  all  sta¬ 
tistical  realizations  of  randomly  distributed  microcracks.  Batchelor  (1970),  Batchelor  and  Green 
(1972),  and  Hinch  (1977)  applied  this  approach  to  the  study  of  fluid  suspensions  within  the  frame- 
woric  of  pairwise  (second  order)  interaction.  The  ensemble  average  approach  with  pairwise  in¬ 
teraction  was  later  applied  to  composite  materials  with  interacting  inclusions  (inhomogeneities)  by 
Willis  and  Acton  (1976),  as  well  as  Chen  and  Acrivos  (1978a,b).  Recently,  Chen  and  Ju  (1991)  and 
Ju  and  Chen  (1991)  proposed  two-dimensional  (second-order  and  higher-order)  micromechanical 
damage  theories  for  brittle  solids  with  interacting  slit  microcracks  by  employing  the  ensemble  av¬ 
erage  method  and  micromechanical  fracture  mechanics.  It  is  noted  that  local  displacements,  strains 
and  stresses  vary  with  positions  within  a  RVE. 

For  simplicity,  we  consider  a  two-phase  composite  composed  of  a  linear  elastic  matrix  and 
many  randomly  located  penny-shaped  microcracks.  The  local  strain  tensor  at  a  point  x  within  the 
RVE  takes  the  form: 

e(x)  =  S°  :  <r(x)  +  e*(x,  C)  (1) 

where  c  and  tr  denote  the  local  strain  and  stress,  respectively;  c*  is  the  perturbed  strain  due  to  the 
existence  and  interaction  of  microcracks;  and  C  denotes  the  set  of  all  possible  configurations  of 
microcracks.  It  is  emphasized  that  e*  is  zero  if  x  is  a  point  in  the  matrix  and  non-zero  if  x  is  a  point 
on  the  microcrack  surfaces.  Taking  the  ensemble  average  over  Eq.  (1),  we  arrive  at 

(£)(x)  =  S®  :  (<r)(x)  +  (e*)(x)  (2) 
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where  the  angle  brackets  (•}  signify  the  ensemble  average. 


^  Throughout  the  development  of  this  paper,  for  simplicity,  we  shall  assume  that  the  solid  is 

locally  homogeneous  (Hinch  (1977))  and  penny-shaped  microcracks  do  not  intersect  one  another. 
Local  homogeneity  implies  that  all  probability  density  functions  (PDF)  do  not  vary  under  small 
translation  on  a  macroscopic  length  scale.  It  can  be  shown  (Chen  and  Ju  (1991))  that,  in  the  case 
®  of  microcracks,  the  perturbed  strain  can  be  expressed  as 

(c-)(x)  =  fix)  ^  l(IIu]I  ®  n  -h  n  (2)  |Iu]l)(x'|x,  G)fiQ)  dQdS,  (3) 

% 

0  Here,  G  =  (u,n)  indicates  the  microcrack  radius  a  and  orientation  n.  Furthermore,  x'  denotes  a 

point  on  the  surfaces  (5,)  of  a  microcrack  centered  at  x;  fix)  is  the  PDF  for  a  microcrack  being 
centered  at  x;  HuJ  is  the  vector  of  microcrack  opening  displacements;  and  fiG)  is  the  PDF  for  a 
microcrack  with  a  geometry  G-  If  all  penny-shaped  microcracks  are  aligned  (parallel)  and  of  equal 
^  size,  then  there  is  no  variation  in  G.  This  could  correspond  to,  for  example,  microcracks  generated 

by  fiber  breaks  in  unidirectional  fiber  composites  (see,  e.g..  Laws  and  Dvorak  (1987)).  In  this  event, 
Eq.  (3)  can  be  simplified  as  follows 


le‘)(x)  =  f (x) 


f  I 


(4) 


It  is  well  known  that,  for  an  open  penny-shaped  microcrack  with  radius  a  embedded  in  an 
infinite  linear  elastic  isotropic  matrix,  the  microcrack  opening  displacements  at  x'  (at  a  distance  p 
from  the  center  of  the  microcrack)  are: 


8(1  -  u^) 
7rE(2  —  u) 


(5) 


where  E  and »/  =  the  Young’s  modulus  and  Poisson’s  ratio  of  the  virgin  matrix  material,  respectively. 
Moreover,  p,  s  and  t  are  the  ^-direction  normal,  the  i-direction  shear  and  the  y-direction  shear 
stresses  projected  on  the  microcrack  surface  in  its  local  coordinates;  see  Fig.  1  for  a  schematic  plot. 

If  all  (open)  microcracks  are  aligned  (parallel),  then  we  can  define  n  =  (0,0, 1)^.  By  sub¬ 
stituting  Eq.  (5)  into  (4)  and  caning  out  the  integration,  we  obtain  (assuming  equal  microcrack 
size) 


69 


where  (the  Voigt’s  notation) 


{  (et)  ("xx) 

IS  kJ 

ja‘\  =  J  ^^3/  ,  —  ,  , 

-  1  (el)  (  =  i  2(£:,)  f 

H  ^K) 

I  K)  ^  I  2{eU)  J 

In  addition,  g  is  the  transformation  matrix  and  T  is  the  local  stress  vector: 


0  0  O' 

0  0  0 

2-1/  0  0 

0  0  0 


■H 


In  the  event  of  distributed  (nonuniform)  microcrack  lengths  and  orientations,  similar  (though  more 
complicated)  expressions  can  be  constructed  accordingly. 

On  the  other  hand,  the  local  stress  vector  (T)  can  be  shown  to  be  (see  Eq.  (14)  in  Sec.  2.2): 


T  =  T°°  +  T  = 


■{•:)•{! 


Here,  T°®  denotes  the  unperturbed  local  stress  vector  due  to  remote  loading,  and  T  denotes  the 
local  stress  perturbation  due  to  three-dimensional  microcrack  interactions.  In  what  follows,  atten¬ 
tion  is  focused  on  pairwise  microcrack  interactions.  Higher-order  microcrack  interactions  will  be 
discussed  in  Sec.  5.  By  assuming  that  all  microcracks  are  aligned  with  a  chosen  global  coordinate 
system,  the  stress  T°°  due  to  far  field  loads  can  be  expressed  as 


=  Ko  •  r' 


where 


0  0  1  0  0  0 
Ko=  0  0  0  0  0  1 
0  0  0  0  1  0 


= 


The  ensemble  average  of  local  stress  perturbation,  on  the  other  hand,  takes  the  form  (assuming 
uniform  size  and  aligned  orientation): 


(t)(«)=  /(t)(: 


;xjxi)/(xi|x)  dxx 


70 


where  )  =  the  ensemble-average  stress  perturbation  for  a  microcrack  centered  at  x  over  the 

subclass  of  realizations  having  a  microcrack  centered  at  Xi;  and  /(xi  |x)  =  the  conditional  PDF  for 
finding  a  microcrack  centered  at  x\  given  a  microcrack  centered  at  x.  Further,  H  designates  the 
active  (qpen)  integration  domain  which  depends  on  the  loading  conditions. 

The  conditional  PDF  /(xi|x)  can  be  simplified  to  /(xi)  if  microcracks  do  not  intersect  and 
reasonable  randomness  holds  (Hinch  (1977)).  The  local  homogeneity  assumption  enables  us  to 
further  approximate  f{X\)  by  /(x).  Therefore.  Eq.  (12)  can  be  recast  as: 

(t)(x)  =  /(x)jf(t)(xlx,)dx,  (13) 

The  quantity  T(xlx])  corresponding  to  the  pairwise  microcrack  interaction  will  be  the  main  subject 
of  the  next  section. 

ni.2.2.  Approximate  explicit  solutions  for  pairwise  interaction  of  aligned  microcracks 

The  objective  of  this  section  is  to  construct  approximate  closed-form  explicit  expressions  for 
perturbed  stresses  T  due  to  the  two-microcrack  interaction  so  that  the  ensemble-average  formalism 
proposed  in  Sec.  2.1  can  be  realized.  It  is  actually  possible  to  derive  explicit  expressions  for  T  in 
the  interaction  problem  involving  two  arbitrarily  oriented  and  located  penny-shaped  microcracks. 
However,  general  solutions  of  the  iwo-microcrack  interaction  problem  are  rather  complicated  and 
no  reasonably  compact  closed-form  explicit  solutions  can  be  presented  within  normal  journal  page 
limit.  Therefore,  for  demonsuation  purpose,  attention  will  be  foucused  on  explicit  solutions  of 
two  randomly  located,  aligned  (parallel),  equal-sized,  penny-shaped  microcracks  embedded  in  an 
infinite,  linear  elastic  isotropic  matrix.  Nonetheless,  in  a  fonheoming  paper,  arbitrary  microcrack 
orientations,  locations  and  sizes  will  be  accommodated  through  numerical  ‘‘microcrack  interaction 
matrix”  and  the  ensemble  average  approach. 

The  “pseudo-traction”  method  is  adopted  here  to  derive  approximate  expressions  for  T  since 
exact  solutions  are  not  yet  available.  For  mathematical  simplicity,  only  the  first  term  of  Taylor’s 
expansion  of  the  local  stress  field  is  used  to  represent  the  average  stress  across  the  microcrack 
surface.  Higher-order  terms  in  polynomial  expansions  may  be  included  if  desired,  however,  at  the 
high  cost  of  a  much  larger  system  of  equations  and  much  more  complicated  analytical  expressions. 
Stemming  from  a  different  viewpoint,  the  more  accurate  “transmission  factor"-typc  formulation 
proposed  by  Kachanov  (1987)  is  well  suited  for  deterministic  numerical  computations.  Within  the 


framework  of  statistical  ep®'*rp.ble  average,  nevertheless,  closed-form  explicit  solutions  for  stress- 
interaction  are  preferred. 

Figure  1  shows  the  local  coordinate  systems  for  microcracks  1  and  2  of  radius  a.  The  z-axis 
is  chosen  as  the  direction  normal  to  the  microcrack  surface.  In  accord  with  the  pseudo-traction 
concept,  the  problem  of  two  interacting  microcracks  subjected  to  far  field  stresses  can  be  decom¬ 
posed  into  a  homogeneous  problem  and  two  sub-problems  (sec  also  Chen  and  Ju  (1991)).  In  the 
homogeneous  problem,  a  microcrack-free  solid  is  subjected  to  applied  stresses  at  far  field.  In  the 
sub-problem  j  (j  =  2),  an  infinitely  extended  solid  contains  only  one  penny-shaped  microcrack 

and  is  subjected  to  zero  remote  stress  at  infinity.  Since  stresses  vanish  on  microcrack  surfaces  (S^), 
the  following  boundary  conditions  must  be  satisfied  0  =  1. 2): 

-Pj+Pr+P;=0  ;  +  Sj  =  0  ;  -0+<f  +  i'i=O  (14) 


In  the  first  sub-problem  (containing  only  the  ‘microcrack  1  ’),  let  us  define  a  cylindrical  coordi¬ 
nate  system  with  the  center  of  the  ‘microcrack  1  ’  as  its  origin.  Accordingly,  the  center  location  of 
the  ‘microcrack  2’  can  be  characterized  by  (p,  4,  z);  see  Fig.  1.  For  convenience,  we  shall  introduce 
the  following  definitions: 

^  ^ yy  ♦  ^2  — -  ^xx  ^yy  (15) 


—  C zz  »  (7xz 


If  the  surfaces  of  the  ‘microcrack  1’  are  subjected  to  applied  normal  stress  then  perturbed 
stresses  at  the  ‘microcrack  2’  location  are  given  by  (Fabrikant  (1989,  p.  252-257)): 


cr,  = 


.  n  X  f  <^(^2  -  .  _l  a^i2a^  +  2z^  -  3p^)] 

\hj\  (/2  _  /2)3(/2  _  ^2)‘/2 


_/(i)  _ 


{l-2u) 


[aH6ll  -  21]  +  p^)  ~  5/1 
(/?  -  -  ah 


2p^  fa 


{(i  - 


.  _i  /  a  \  az^  [/J  +  a^(2a^  +  2z^  —  3p^)]  1 

UJ  (l|  -  P,)HP  -  I 


=  2p"*  \a\in  -  5p^) ..  i;l 

k(p2  - 


where 
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Furthermore,  when  the  surfaces  of  the  ‘mirrocrack  1  ’  are  loaded  by  the  combined  shear  stresses 
j.(i)  =  ^(1)  +  the  perturbed  stresses  at  the  ‘microcrack  2’  location  are  (Fabrikant  (1989,  p.  257- 
261)): 


=  2(fWe»^r»>e-‘^) 
^  7r(2  —  u) 

—  2c’^  f 


^  h(^-w  * - - / 


_ 

<7  2- 


7r(2  —  v) 


I  ^  ^  wl  -  /?)  hill  -lb  1^2  (r  +  T  e  )J  J 


„(i)  ^  2(fd)e»»  -h  rd)e-«^)  zhOl  -  [q^(4/^  -  5pb  -^  /?] 

7r(2  -  u)  I2HI  -  /^)3 


•  yj* - ar- 

,  2,1/2  ,  (fe“ 


—  1?)''^^  [/f  +  a^(2a^  +  2z^  —  3p^)] 

oFW 


where  fd)  is  the  complex  conjugate  of  rd^  The  total  perturbed  stresses  due  to  combined  normal 
loading  pd)  and  shear  loading  rd)  are  therefore  obtained: 

^<l)  _  ^  /T"d)  .  ^0)  _  ^'d)  .  ^//O) 

(Jt  SS  ^  .  "T  t  1  <Ti  35  ^  *  T*  o 

111.  2  2  2  ^20) 

<rd)  =  •  rd)  =  /’>  +  r"<’> 

It  can  be  shown  that  the  Canesian  stress  components  are 

</<■>  =  5  W'’  +  R«(4'’)]  ;  5  W"  -  RcC"?’)] 

;  4"  =  «'> 

<7d>  =  Im(rd))  ;  ad)  =  Re(Td)) 

where  “Re”  and  “Im”  are  the  real  and  imaginary  parts,  respectively,  of  a  complex  variable. 

Eq.  (21)  can  be  recapitulated  into  the  following  matrix  form: 


r'*'  = 


0) 

<7yj, 

(Tzz 

•  = 

^ry 

(^yz 

.  . 

Cl  *02  1  d) 

2  2  2 

fel— -62  Cl-^C2  01  -»02 

2  2  2 

bi  C3  dj 

^  £i.  ii 

2  2  2 

65  C5  <(5 

hh  C6  db  . 


Ifl' 


where  definitions  of  the  parameters  6,,  c,,  and  <f,  can  be  found  in  Appendix  I.  Symbolically,  we 
can  express  rd)  as  a  function  of  coordinates  and  normal  and  shear  stresses  as  follows: 


=  T{p,d>,z\puS\Ai) 
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Since  micrcx:racks  1  and  2  are  aligned,  the  perturbed  n''”T'.'!l  and  shear  stresses  along  the  sur¬ 
faces  at  the  ‘microcrack  2’  location  arc 


^  =  e,  ■  •  e,  = 

52  =  e,  •  •  Ci  =  ai'J 


(24) 


<2  =  •  r 


<».e 


where  e*.  e,,,  and  are  the  unit  base  vectors  in  the  Cartesian  (x,  y,  and  z)  coordinates. 

Similarly,  in  the  sub-problem  2  (containing  only  the  ‘microcrack  2’),  it  can  be  shown  that  the 
perturbed  normal  and  shear  stresses  along  the  surfaces  at  the  ‘microcrack  1  ’  location  arc 


Pi  =  ; 


where  symbolic  representations  similar  to  Eq.  (23)  and  (24)  have  been  employed: 

=  r(/>,  It  -4),  z\  p2,  -52,  <2) 

Eq.  (24)  together  with  (25)  then  leads  to 
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(27) 


For  convenience  in  the  following  derivations,  let  us  define 


T,_2  =  < 


'  Pi  ' 

fpri 

'  Pi  ' 

Si 

sV 

ti 

tr 

^00 

-  ;  t,_2  =  < 
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(28) 


and  a  =  the  6  X  6  "microcrack  interaction  matrix”  in  Eq.  (27).  It  is  noteworthy  that  our  a  matrix, 
in  fact,  corresponds  to  the  A  matrix  (the  “transmission  factor”  or  “crack  interaction  matrix”)  in 
Kachanov  (1987)  and  Kachanov  and  Laures  (1989).  With  these  notations  at  hand,  Eq.  (27)  can  be 
rewritten  as 

ti_2  =  aT,_2  (29) 


Since  Ti_2  =  Tflj  +Ti_2.  T)_2  can  be  solved  from  Eq.  (29): 


t,_2  =  K  •  Tit2  ;  where  K  =  a  ■  (I  -  a)"' 


(30) 


Further,  we  define  T  =  (pi,  si,  and  Ki  =  the  first  three  rows  of  K.  Hen''**,  from  Eq.  (30).  the 
perturbed  stresses  on  surfaces  of  a  microcrack  due  to  the  existence  of  the  second  microcrack  arc: 


f  =  Ki  T 


OO 

1-2 


(31) 


Finally,  since  the  two  microcracks  under  consideration  are  aligned,  it  can  be  easily  shown  that 


^OO  _ 

»  1-2  = 


0  0  1  0  0  0' 

0  0  0  0  0  1 

0  0  0  0  1  0 

0  0  1  0  0  0 

0  0  0  0  0  1 

0  0  0  0  1  0 


=  Kj  •  t” 


From  Eq.  (32)  and  (31),  we  conclude  that  the  average  penurbed  stress  vector  T  over  surfaces  of  a 
miciocrack  is  simply: 

f  =  K  •  ;  where  K  =  K,  ■  Kj  (33) 


Substitution  of  Eq.  (33)  into  (13)  then  renders  realizations  for  (T)  in  the  previous  section.  There- 
fwe,  the  ensemble  average  approach  is  completely  defined.  For  non-aligned  penny-shaped  micro¬ 
cracks,  more  complicated  derivations  will  be  involved. 


1II.2.3.  Some  test  problems  for  two-microcrack  interaction 


A  number  of  test  problems  are  considered  in  this  section  to  examine  the  performance  of  the  ap¬ 
proximate  analytical  solutions  presented  in  Sec.  2.2.  These  include  normal  and  shear  loadings  for 
two  aligned  coplanar  or  stacked  penny-shaped  microcracks.  It  is  not  our  intention,  however,  to  pro¬ 
pose  highly  accurate  analytical  solutions  to  compute  local  stresses  and  stress  intensity  factors  (SIFs) 
at  all  points  on  microcrack  surfaces  for  a  two-microcrack  interaction  problem.  Instead,  reasonably 
accurate  analytical  average  perturbed  stresses  (in  terms  of  elementary  functions)  over  microcrack 
surfaces  are  sought  in  order  to  exploit  the  ensemble  average  approach.  In  fact,  if  one  is  interested 
in  deterministic  3-D  microcrack  interaction,  excellent  numerical  method  has  been  proposed  by 
Kachanov  and  Laures  (1989).  It  is  noted  that  key  steps  in  Kachanov  and  Laures’  (1989)  method 
also  focus  on  the  computation  of  average  pairwise  stress  “transmission  factors”  and  average  trac¬ 
tions.  Once  average  tractions  become  known,  one  can  certainly  compute  projected  local  stresses 
and  SIFs  at  any  point  on  microcrack  surfaces.  In  general,  the  simple  analytical  solutions  presented 
in  Sec.  2.2  are  not  as  accurate  as  those  proposed  in  Kachanov  and  Laures  (1989).  Nonetheless,  the 
latter  relies  on  extensive  numerical  computations  of  transmission  factors  for  all  points  on  micro¬ 
crack  surfaces  and  is  therefore  not  employed  here. 


It  is  emphasized  that,  given  an  existing  microcrack,  the  location  of  the  second  microcrack  is 
random  within  the  ensemble  average  framework.  Although  stress  interactions  between  two  mod¬ 
erately  spaced  random  microcracks  are  not  very  strong,  their  cumulative  effects  are  important  due 
to  high  spatial  probabilities.  On  the  other  hand,  stress  interactions  between  two  closely  spaced  ran¬ 
dom  microcracks  are  strong  yet  their  spatial  probabilities  are  lower.  Therefore,  contributions  from 
both  closely  and  moderately  spaced  microcracks  should  be  accounted  for  in  the  ensemble  average 
approach.  The  ensemble-average  spacing  of  microcrack  arrays,  clearly,  depends  on  microcrack 
concentration. 

Cask  I:  Two  equal-sized  coplanar  microcracks  under  normal  loading.  Though  numerical  re¬ 
sults  for  SIFs  are  available  for  two  coplanar  microcracks  under  normal  and  shear  loadings  (Fab- 
rikant  (1987,  1989)),  exact  results  for  average  tractions  projected  over  microcrack  surfaces  are  not 
documented.  Nevertheless,  Kachanov  and  Laures  (1989)  show  that  their  results  for  SIFs  are  very 
close  to  those  of  Fabrikant  (1987, 1989)  for  two  equal-sized  coplanar  microcracks.  Therefore,  the 
average  tractions  computed  by  Kachanov  and  Laures  (1989)  should  be  quite  accurate  in  the  case 
of  coplanar  microcracks.  The  ratios  of  perturbed  vs.  far-field  average  normal  stresses  iplp°°)  ob¬ 
tained  by  the  proposed  method  are  listed  in  Table  1  for  various  l/2a  values  (the  Poisson  ratio  u 
=  0.25).  Here,  I  signifies  the  center-to-center  distance  between  two  microcracks  and  2c  is  the  mi¬ 
crocrack  size.  For  example,  I /2a  =  1.00025  means  that  the  smallest  distance  between  microcrack 
tips  is  only  0.00025a.  The  results  reported  in  Kachanov  and  Laures  (1989)  and  relative  differences 
between  the  two  results  are  also  given  in  Table  1  for  comparison  purpose.  From  Table  1,  it  is  clear 
that  interaction  renders  average  stress  amplification.  The  effect  of  microcrack  interaction  decays 
as  the  distance  between  two  microcracks  increases.  Moreover,  the  differences  between  our  simple 
calculations  and  those  of  Kachanov  and  Laures  (1989)  are  small. 

Case  II:  Two  equal-sized  coplanar  microcracks  under  shear  loading.  The  ratios  of  perturbed 
vs.  far-field  average  shear  stresses  (t/t°°)  obtained  by  the  present  method  are  listed  in  Table  2 
for  various  //2a  values  (the  Poisson  ratio  u  =  0.5).  From  Table  2,  it  is  seen  again  that  interaction 
renders  average  stress  amplification.  TTie  differences  between  our  simple  calculations  and  those  of 
Kachanov  and  Laures  (1989)  are  small. 

Case  ni:  Two  equal-sized  stacked  microcracks  under  normal  loading.  The  ratios  of  perttu-bed 
vs.  far-field  average  normal  stresses  {plp°°)  are  listed  in  Table  3  for  various  //2a  values  (the  Pois¬ 
son  ratio  V  =  0.25).  Obviously,  from  Table  3,  interaction  renders  strong  average  stress  shielding. 


The  degree  of  shielding  in  this  case  is  much  stronger  than  the  degree  of  amplification  in  Cases  I  or 
II.  Fabrikant  (1989)  did  not  provide  numerical  results  for  SIFs  nor  average  stresses  for  two  stacked 
microcracks  under  normal  or  shear  loadings.  Therefore,  exact  results  are  not  available.  For  com¬ 
parison,  the  differences  between  our  simple  calculations  and  those  of  Kachanov  and  Laures  (1989) 
are  listed  in  Table  3.  It  is  observed  that  the  difference  is  not  very  significant  if  //2a  is  greater  than 
0.5.  The  difference  increases  as  the  two  microcracks  move  closer,  due  to  sharp  variations  of  stress 
fields  in  the  close  neighborhood  of  interacting  microcracks. 

Case  IV:  Two  equal-sized  stacked  microcracks  under  shear  loading.  The  ratios  of  perturbed  vs. 
far-field'average  normal  stresses  (plp°°)  are  given  in  Table  4  for  various  Ijla  values  (the  Poisson 
ratio  V  -  0.2S).  From  Table  4,  interaction  renders  average  stress  shielding  when  the  distance  I  is 
small  and  very  weak  stress  amplification  when  Ijla  >  0.5.  The  differences  between  our  simple 
calculations  and  those  of  Kachanov  and  Laures  (1989)  are  small  for  //2a  greater  than  0.25. 

Since  pairwise  microcrack  interactions  are  random  (probabilistic)  within  the  framewOTk  of  the 
ensemble  average  approach,  the  errors  associated  with  the  present  approximate  analysis  should  be 
statistically  averaged  over  all  possible  realizations.  Therefore,  the  ensemble-average  error  of  the 
present  method  should  be  small  as  long  as  microcrack  concentration  is  not  extremely  high. 


IIL3.  Eniective  moduli  of  brittle  solids  with  interacting  microcracks 

With  the  ensemble  average  framework  and  analytical  pairwise  interaction  solutions  at  hand, 
we  are  now  ready  to  construct  the  ensemble-average  constitutive  equations  of  brittle  solids  with 
many  randomly  located,  aligned  interacting  microcracks.  By  substituting  Eq.  (10)  and  (33)  into  (9) 
and  taking  the  ensemble  average,  we  arrive  at 


(T)=(Ko-*-/(x)(K)).r“ 


(34) 


where 


J  Kr^  sin  V*  drdxpdO 


(35) 


In  the  foregoing  equations,  it  is  implicitly  assumed  that  all  microcracks  are  of  equal  size  and  of  same 
orientation.  In  addition,  the  spherical  coordinate  system  (r,  0, 6)  is  used  to  describe  the  random 
location  (xi)  of  the  second  microcrack  relative  to  the  first  random  microcrack  centered  at  x.  Note 
that  0  varies  from  0  to  tt  and  6  ranges  from  0  to  2n.  If  we  normalize  r  with  resp)ect  to  the  microcrack 
radius  a  (i.e.,  (  =  r/a),  then  Eq.  (35)  can  be  recast  as 


(K)  = 


s\n^didxl>de 


(36) 


Combining  Eq.  (6),  (34)  and  (36),  we  obtain  the  local  ensemble-averaged  damage-induced 
strain  (at  a  typical  point  x): 

(€*)(x)=(S*)(x)r~  (37) 

Here  the  ensemble-averaged,  damage-induced  local  compliance  has  two  components: 

(S*)(x)  =  (S“)(x)+(S*^)(x)  (38) 


where 


<S*')(x) 

{S*^)(x) 


16(1  -  j/2) 

3Ei2  -  u) 
16(1  - 1/2) 
3Ei2  -  u) 


/(x)a’g  •  Ko 
/2(x)a‘g  •  (K) 


(39) 


It  is  noteworthy  that  Eq.  (39)  actually  reveals  physical  nonlocal  effects  in  constitutive  equations 


through  microcrack  interaction  and  the  ensemble  averaging  process.  That  is,  the  stress-strain  laws 


at  a  material  point  x  within  a  RVE  depend  on  the  constitutive  laws  of  all  neighboring  points.  This 
is  a  physical  nonlocal  approach,  at  variance  with  ti.e  postulated  nonlocal  theories  due  to  Eringcn 


and  Edelen  (1972). 


To  obtain  volume  averaged  moduli  due  to  microcracks  within  a  RVE,  one  simply  applies  the 
volume-average  operator  to  Eq.  (38)-(39).  As  a  consequence,  we  have 


where 


TcTT  -  L«>  -  K  o’ 

'  '  3£(2-i/)*  “““  V 

7^  -  liO  -  /if,  „s  fv  A») ''« 
'  3£(2^*  '  '  V - 

Letbs  assume  that  there  are  N  microcracks  in  the  RVE;  i.e., 


/  /(x)  dx  =  N  (42) 

Jv 

Further,  consider  the  case  in  which  the  variance  of  the  PDF  (/ (x))  for  locations  of  microcracks  is 
small  (e.g.,  uniform  probability).  We  may  therefore  write 


Jv  V  V 


Substitution  of  Eq.  (42)  and  (43)  into  (41)  then  leads  to 

16(1 -r/^), 

^  ^  3£’(2-z/)^  *^‘*' 


where 


is  the  (volume-averaged)  microcrack  concentration  parameter.  We  emphasize  that  our  definition 
of  u  is  different  from  that  defined  by  Budiansky  and  O’Connell  (1976).  The  two  definitions  differ 
by  a  scalar  factor  47r/3.  It  is  noted  that  (S*’)  actually  corresponds  to  the  first-order  contribution 
due  to  non-interacting  microcracks;  i.e.,  the  simple  Taylor’s  model  is  recovered.  Moreover,  (S*^) 
represents  the  second-order  (in  w^)  contribution  due  to  pairwise  microcrack  interaction. 

Finally,  the  overall  (volume-ensemble  averaged)  effective  moduli  for  a  microcrack-weakened 
solid  is  obtained  by  adding  the  elastic  compliance  S°  to  (S*): 

(46) 
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In  panic'i’r',  in  the  case  of  3-D  linear  isotropic  elasticity,  S°  reads 

■  1  -j/  -1/  0  0  O' 

—V  1  —u  0  0  0 

co  _  ^  —V  —V  10  0  0 

^  £  0  0  0  2(1  + 1/)  0  0 

0  0  0  0  2(1+1/)  0 

0  0  0  0  0  2(1+1/) 


(47) 


It  is  emphasized  that  local  macroscopic  constitutive  laws  for  a  RVE  are  recovered  by  the  volume  av¬ 
eraging  process.  However,  if  one  stops  at  the  ensemble  averaging  level  without  further  performing 
the  volume  averaging  process,  the  constitutive  relations  are  nonlocal. 


# 
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I1I.4.  Some  numerical  examples 


# 


In  this  section,  some  numerical  examples  involving  aligned  microcracks  arc  presented  to  il¬ 
lustrate  the  proposed  ensemble-volume  average  approach  in  three-dimension.  In  the  first  example, 
effects  of  microcrack  interactions  arc  neglected  and  the  well-known  Taylor’s  model  is  recovered. 
Subsequently,  the  proposed  second-order  microcrack  interaction  model  is  implemented.  Finally, 
we  compare  the  results  of  the  prcsent  approach  with  some  existing  methods,  including  the  Tay¬ 
lor’s  model,  the  self-consistent  method  and  the  differential  scheme.  Closed  (nonactive)  microcrack 
contributions  to  compliances  arc  neglected  in  this  work. 

nL4.1.  bilute  non-interacting  aligned  microcracks 


Let  us  consider  the  case  in  which  effects  of  micrxxrrack  interactions  arc  totally  neglected.  This 
can  be  done  by  simply  dropping  the  term  (S*^)  in  Eq.  (46): 

^  =  (48) 


From  Eq.  (8),  (1 1),  and  (44),  we  obtain 


(S-i)  = 


16(1  - 
3£’(2  —  i/) 
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0  0  0  0  0  0 

0  0  2-1/  0  0  0 


0  0 
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0 

0 
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0  0  0 
0  2  0 
0  0  2 


(49) 


Combination  of  Eq.  (47)  and  (49)  then  yields 
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(51) 


(52) 
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It  is  observed  that  only  the  three  compliance  components  in  Eq.  (51)  arc  changed  due  to  the 
presence  of  randomly  located,  aligned,  penny-shaped  microcracks.  One  notes  that  Eqs.  (50)-(52) 
indeed  recover  the  well-known  Taylor’s  model.  Further,  the  effective  constitutive  equations  belong 
to  the  category  of  transverse  isotropy.  This  is  typical  of  fiber  breaks  in  unidirectionally  reinforced 
fiber  composites. 

nL4.2.  Aligned  interacting  penny-shaped  microcracks 

We  now  focus  on  the  proposed  ensemble-volume  averaged,  pairwise  interaction  model.  For 
a  given  microcrack  concentration  w,  the  overall  effective  moduli  in  Eq.  (46)  can  be  evaluated  by 
carrying  out  the  integration  in  Eq.  (36).  That  is,  one  needs  to  perform  integration  for  the  matrix  K 
over  the  active  microcrack  domain  (^,  j/',  6).  It  is  emphasized  that  the  values  of  K  decay  rapidly 
as  the  distance  between  a  pair  of  randomly  located  microcracks  increases.  Therefore,  compliance 
contributions  due  to  remote  integration  region  can  be  neglected. 

The  integration  in  Eq.  (36)  can  be  effectively  computed  by  the  Gauss  quadrature  scheme  with 
three  independent  variables:  i,  and  Q\  see  Table  5  for  the  convergence  behavior.  Our  numerical 
experiments  show  that  use  of  ^mar  =  20  (or  r,nax  =  20a)  is  quite  acceptable;  i.e.,  contributions 
from  the  ^  >  20  domain  can  be  neglected.  In  Table  5,  the  minimum  radius  of  integration  for  r 

A  A 

(or  0  is  assumed  to  be  r„„„  =  2a  (or  =  2).  The  only  nonzero  components  in  (K)  are  {K\j,), 
{Kta)  (■^^  35)  latter  two  are  equal).  It  is  observed  that  the  normal  component  ( A'13)  is  much 
greater  than  the  shear  components  { A'26)  and  ( A35).  Furthermore,  use  of  =  2a  implies  that  the 
minimum  allowable  distance  between  any  two  microcrack-centers  is  2a.  Therefore,  according  to 
the  “face-center  cubic”  calculation,  the  maximum  allowable  u;  is  l/4>/2  ~  0.1768  for  =  2a. 
Note  that  the  maximum  allowable  microcrack  density  would  be  0.741  if  a;  is  defined  as 
(Budiansky  and  O’Connell,  1976). 

To  achieve  higher  maximum  permissible  microcrack  densities,  various  minimum  radii  of  in¬ 
tegration  are  used  to  compute  (K);  see  Table  6  for  details.  The  value  of  ^mor  is  chosen  as 
20.0  in  Table  6.  We  observe  significant  increases  in  both  (A'n)  and  {Ku)  as  4m.n  decreases.  The 
maximum  microcrack  concentration  w  corresponding  to  the  =  1 .  la  case  is  found  to  be  1 .0627 
(or  4.451  if  the  definition  of  Budiansky  and  O’Connell  (1976)  is  employed).  It  is  emphasized  that 
Table  6  is  not  a  summary  of  convergence  behavior,  but  a  display  of  (K)  for  different  maximum 
allowable  microcrack  densities. 
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Eq.  (8)  and  (44),  we  obtain 
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(53) 


Adding  Eq.  (53)  to  (50),  we  thus  arrive  at  the  expression  for  overall  effective  compliances  with  the 
second-order  microcrack  interacdon.  Apart  firom  the  three  compliance  components  (S33),  (S55) 
and  (566).  it  is  seen  that  all  other  components  are  idendcal  to  the  elasdc  components.  In  particular, 
the  three  compliance  components  which  change  under  microcrack  interaction  take  the  following 
explicit  forms: 


(533)  =  53^3  + 


16(1  -  1/2) 
3E 


(u;  +  k\(jLj^) 


(Sjs)  =  (566)  =  555  +  — ~  (w  + 


(54) 


3£:(2  - 1/) 

Clearly,  the  relative  weight  of  the  second-order  compliance  terms  (in  to  the  first-order  terms 
depends  on  the  values  of  microcrack  densities.  For  dilute  (very  low)  w,  the  second-order  terms  arc 
negligible.  By  contrast,  the  compliance  contributions  are  significant  for  high  u  values. 


III.4.3.  Comparison  with  some  existing  methods 


The  effects  of  fiber  breaks  and  aligned  penny-shaped  microcracks  on  the  stiffness  of  unidirec¬ 
tional  fiber  composites  have  been  studied  extensively  in  the  literature.  We  refer  to  Laws  and  Dvorak 
(1987)  for  an  excellent  presentation  by  using  the  self-consistent  method  and  the  differential  scheme. 
It  is  noted  that,  in  Laws  and  Dvorak  (1987),  the  microcrack  density  is  defined  as  a  =  %Na?fV; 
i.e.,  their  a  is  equal  to  8-/  in  this  paper.  The  maximum  microcrack  concentration  considered  in 
Laws  and  Dvorak  (1987)  is  u?  =  a/8  =  0.125.  The  macroscopic  (overall)  material  behavior  is,  not 
surprisingly,  transversely  isotropic.  It  should  be  realized  that  the  Taylor’s  model,  the  self-consistent 
method  and  the  differential  scheme  all  belong  to  the  category  of  “effective  medium”  approaches. 
Namely,  theses  methods  do  not  depend  on  the  distributions  of  microcrack  locations  at  all. 


In  the  following  numerical  computations,  the  Young’s  modulus  E  is  taken  as  0.5  MPa  and  the 
shear  modulus  G  as  0.2  MPa  (1/  =  0.25)  for  the  virgin  matrix  material.  TTic  overall  (ensemble- 
volume  averaged)  longitudinal  normal  and  shear  compliances,  (533)  and  (555),  arc  plotted  against 
the  microcrack  concentration  parameters;  in  Figures  2  and  3,  respectively.  The  values  of  hi  and  ^2 
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in  Eq.  (54)  are  taken  as  1 .049792  and  0.355988,  respectively,  which  corresoond  to  (min  =  1-1  (i-c., 
the  maximum  allowable  density  ui  =  1.0627).  Note  that,  once  more,  u;  =  1  corresponds  to  o  =  8 
in  Laws  and  Dvorak  (1987).  For  comparison,  the  results  obtained  by  using  the  Taylor’s  model,  the 
self-consistent  method  and  the  differential  scheme  are  also  displayed  in  Figures  2  and  3.  We  would 
like  to  comment  that:  (a)  the  Taylor’s  model  is  really  suitable  for  dilute  microcrack  concentrations; 

(b)  the  self-consistent  method  and  the  differential  scheme  are  suitable  for  low  or  moderate  w;  and 

(c)  the  proposed  statistical  pairwise  microcrack  interaction  model  is  suitable  for  moderately  high 
(not  extremely  high)  w. 


Figtlres  4  and  5  show  the  effects  of  microcrack  density  u;  on  the  normalized  longitudinal 
Young’s  modulus  EifE  and  shear  modulus  Gl/G  for  four  different  models.  We  have  employed 
the  standard  notation;  El  =  ^/{Sn)  and  Gl  =  1/(555).  In  addition,  it  is  straightforward  to 
express  the  normalized  moduli  in  Figures  4  and  5  as  follows  for  the  Taylor’s  model: 


El  _  1  .  Gl  ^  1 

E  1+5u;  ’  G  1-h2.28571u; 


(55) 


and  for  the  prf "  nt  model: 

El  _  1  .  Gl  _  1 

E  ~  1  -^5a;(l  1.05^')  ’  G  ~  1  +2.2857M]  +0.356u,') 


(56) 
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ni.5.  A  higher-order  ensemble-average  formulation  of  microcrack  interaction 

So  far,  attention  has  been  focused  on  the  second-order  statistical  model  based  on  the  concept 
of  pairwise  microcrack  interaction.  In  this  section,  extension  to  accommodate  higher-oider  in¬ 
teraction  mechanism  is  considered  within  the  ensemble-volume  average  framework;  sec  also  Ju 
and  Chen  (1991)  for  two-dimensional  problems.  TTie  proposed  higher-order  formulation  hinges 
on  the  derivation  of  higher-order  (the  third  or  above)  corrections  on  local  stress  fields  (T)  due  to 
n-microcrack  interaction  (n  >  3);  see  Eq.  (9>-(l2). 

As  an  illustration,  let  us  consider  a  threc-microcrack  (the  third-order  in  interaction  mcch- 
anism  within  the  proposed  frame'vork.  For  simplicity,  we  assume  that  all  microcracks  are  aligned 
and  of  equal  size.  Following  the  definitions  and  assumptions  previously  described  in  Section  2.1, 
one  can  rephrase  the  local  ensemble  stress  p)erturbation  as  (cf.  Eq.  (9)): 

(t)  +  (t)  =  (|j|)  +  (|f|)  (57) 

where  (T)  is  the  second-order  local  ensemble  stress  perturbation  due  to  pairwise  microciack  inter¬ 
actions  (see  Eq.  (9)  and  (12)),  and  (T)  is  the  third-order  local  ensemble  stress  perturbation  due  to 
3-microcrack  interactions.  In  particular,  (’f)  can  be  expressed  as  (cf.  Eq.  (12)) 

(f)  =  ^(f)(x;xi|x2)/(x2|x;xi)  rfx:  (58) 

Here,  (T)  (x;  Xi  1x2)  is  the  third-order  ensemble-average  stress  perturbation  of  a  microcrack  centered 
at  X,  given  a  microcrack  centered  at  Xj,  over  a  subclass  of  realizations  which  have  a  microcrack 
centered  at  X2.  Further,  /(X2|x;  Xi )  is  the  conditional  probability  density  function  (PDF)  for  finding  a 
microcrack  centered  at  Xi  given  two  microcracks  fixed  at  x  and  at  Xj ,  respectively.  The  conditional 
PDF  /(xalx;  Xi)  can  be  further  simplified  to  /(x)  by  the  assumptions  of  local  homogeneity  and 
reasonable  randomness  (i.e.,  statistical  independence).  The  active  (open)  integration  domain  H 
depends  on  loading  conditions. 

For  clarity,  let  us  express  (T)  and  (e*)(x)  as  follows  (cf.  Eq.  (34),  (37)  and  (38),  assuming  local 
homogeneity  and  reasonable  randomness): 

(T)  =  (T-  f  t)  =  (Ko  +  /(x)(K)  f  (x){K'))  ■  r'-  (59) 

(e*>(x)  =  { {S*‘)(x)  +  (S*^)(x)  {S*5)(x))  ■  =  {S*)(x)  •  (60) 
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where  T  =  K'  •  r*'  and  (cf.  Eq.  (39'>'> 

(S*3)(x) 


16(1  -i/2) 


r<x)a’g  .  {K') 


(61) 


3£’(2-i/)- 

In  Eq.  (61),  we  have  defined  (K')  =  1/a^  (K').  We  recall  that  both  Ko  and  <C  are  expressed 
explicitly  in  closed-form  formulas.  Similarly,  K'  (or  T)  can  also  be  constructed  in  closed-form  as 
follows.  One  stares  by  expanding  Eq.  (14)  into  nine  linear  equations  with  j  =  1,2, 3.  Then,  one 
obtains  expressions  similar  to  Eq.  (22),  with  the  understanding  that  permutations  1-2, 2-3,  3-1  are 
involved.  Eq.  (24)  in  Sec.  2.2  is  modified  to  include  the  third  microcrack’s  contribution  to  stress 
perturbations.  Subsequently,  Eq.  (27)  is  expanded  to  a  9  by  9  system  with  a  denoting  a  9  by  9 
coefficient  matrix.  Therefore,  we  arrive  at  explicit  formulas  similar  to  Eq.  (30)-(33);  and  R',  t" 
can  be  expressed  in  closed-form  as  in  the  two-microcrack  interaction  problem.  Finally,  it  is  noted 
that  the  computation  of  (K')  involves  integration  over  the  domain  of  all  possible  positions  of  two 
active  neighboring  microcracks  around  a  specified  penny-shaped  microcrack  (cf.  Eq.  (36)). 

Following  the  same  procedure  presented  in  Sec.  3,  it  can  be  shown  that  (S*^)  in  Eq.  (61)  gives 
rise  to  the  third-order  terms  (in  u;^)  to  overall  compliances  due  to  the  third-order  microcrack  in¬ 
teractions.  Therefore,  a  third-order  statistical  micromcchanical  model  can  actually  be  constructed. 
By  repeating  the  foregoing  procedure,  we  can  formulate  a  complete  (though  very  complicated) 
hierarchical  family  of  statistical  microcrack  theories  of  arbitrary  orders. 
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III.6.  Conclusions 


# 


Based  on  the  concept  of  ensemble-volume  average  and  pairwise  microcrack  interaction,  an 
innovative  three-dimensional  statistical  micromechanical  theory  is  proposed  for  brittle  solids  vi/ith 
many  randomly  located  (but  aligned),  interacting  microcracks.  The  effects  of  interaction-induced 
stress  penurbations  and  random  locations  of  microcracks  are  manifested  by  (T) .  Approximate  ex¬ 
plicit  (analytical)  interaction  solutions  for  two  arbitrarily  located,  aligned  microcracks  are  given 
in  detail.  Therefore,  overall  effective  moduli  of  brittle  solids  with  interacting  microcracks  can  be 
fwmally  derived  on  the  ground  of  statistical  and  micromechanical  informations.  However,  explicit 
closed-fbrm  solutions  to  arbitrarily  located  and  arbitrarily  oriented  microcracks  are  very  involved, 
and  no  reasonably  compact  explicit  forms  can  be  presented  within  normal  page  limit.  In  a  forthcom¬ 
ing  paper,  nonetheless,  general  microcrack  geometry  will  be  fully  accommodated  through  numeri¬ 
cal  computations  of  the  “microcrack  interaciion  matrix”  a  and  the  ensemble  average  approach.  In 
addition,,  it  is  interesting  to  note  that  a  physically  “nonlocar  description  of  the  material  behavior 
can  be  obtained  during  the  ensemble  averaging  process. 

The  proposed  approach  is  fundamentally  different  from  existing  effective  medium  methods 
which  do  not  depend  on  locations  and  configurations  of  microcracks.  Further,  the  proposed  micro- 
crack  interaction  framework  does  not  require  the  use  of  Monte  Carlo  simulations.  Some  numerical 
examples  are  presented  to  illustrate  the  behavior  of  the  proposed  model.  The  resulting  predictions 
are  compared  with  some  existing  methods.  Finally,  a  higher-order  microcrack  interaction  formu¬ 
lation  is  briefly  summarized. 

The  proposed  method  provides  a  simple  framework  to  accommodate  statistical,  micromechan¬ 
ical,  and  interaction  aspects  of  distributed  microcrack  arrays.  Applications  may  be  made,  for  ex¬ 
ample,  to  aligned  matrix  cracks,  fiber  breaks,  and  interlaminar  delaminations  of  brittle  composite 
materials.  Further  research  will  be  needed  in  the  future  to  assess  the  applicability  of  the  proposed 
method  to  practical  engineering  problems. 
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Table  1.  Two  equal-sized  coplanar  microcracks  under  normal  loading  (v  =  0.25) 


Ilia 

Present 

Kachanov 

Diffcrcncc(%) 

1.00025 

1.0354 

1.0837 

4.46 

1.005 

1.0348 

1.0779 

4.00 

1.05 

1.0296 

1.0529 

2.21 

1.1 

1.0251 

1.0398 

1.41 

1.15 

1.0215 

1.0315 

0.97 

1.25 

1.0161 

1.0214 

0.52 

1.5 

1.0088 

1.0104 

0.16 

2.0 

1.0035 

1.0038 

0.03 

2.5 

1.0018 

1.0019 

0.01 
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Table  2.  Two  equal-sized  coplanar  microcracks  under  shear  loading  {y  =  0.5) 


Ijla 

Present 

Kachanov 

Differcncc(%) 

1.005 

1.0681 

1.1017 

3.05 

1.05 

1.0580 

1.0703 

1.15 

1.25 

1.0317 

1.0292 

0.24 

1.5 

1.0174 

1.0144 

0.30 

1.75 

1.0106 

1.0084 

0.22 

2.0 

1.0070 

1.0054 

0.16 

2.5 

1.0035 

1.0026 

0.09 

3.5 

1.0013 

1.0006 

0.07 
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Tables.  Two  equal- sized  slacked  microcracks  under  normal  loading  (v  =  0.25) 


l/2a 

Present 

Kachanov 

Diffcrcnce(%) 

0.05 

0.5004 

0.5583 

10.37 

0.25 

0.5383 

0.6689 

19.52 

0.35 

0.5836 

0.7158 

18.47 

0.5 

0.6667 

0.7777 

14.27 

0.75 

0.7928 

0.8562 

7.40 

1.0 

0.8754 

0.9073 

3.52 

1.5 

0.9505 

0.9588 

0.87 
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Table  4.  Two  equal- sized  stacked  microcracks  under  shear  loading  (i/  =  0.25) 


//2a 

Present 

Kachanov 

Difference(%) 

0.05 

0.5549 

0.6613 

16.09 

0.25 

0.8214 

0.8886 

7.56 

0.5 

1.0002 

0.9837 

1.68 

0.75 

1.0233 

1.0053 

1.79 

1.0 

1.0180 

1.0084 

0.95 

1.25 

1.0120 

1.0072 

0.48 

1.5 

1.0080 

1.0055 

0.25 

2.0 

1.0039 

1.0031 

0.08 

2.5 

1.0021 

1.0026 

0.05 
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Table  5.  Convergence  behavior  of  (K)  vs.  ^mai  (given  ^m.n  =  2) 


^max 

(A'i3> 

{Ku)  =  {A'35: 

4.0 

0.150807 

0.059503 

6.0 

0.168554 

0.066656 

8.0 

0.173029 

0.068467 

10.0 

0.174644 

0.069121 

20.0 

0.176139 

0.069728 

40.0 

0.176327 

0.069804 

80.0 

0.176351 

0.069828 

160.0 

0.176354 

0.069829 

320.0 

0.176354 

0.069829 

640.0 

0.176354 

0.069829 
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Table  6.  Numerical  integration  ot  (K)  for  different  ^rn,n  (given 


,  “  ^mtnl ^ 

(A'26)  =  (A'35 

2.0 

0.176139 

0.069728 

1.75 

0.305726 

0.105974 

1.5 

0.601154 

0.195689 

1.25 

0.850189 

0.281202 

1.125 

0.975302 

0.327491 

1.1 

1.049792 

0.355988 

=  20) 
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111.8.  Appendix  1:  Parameters  for  Eq.  (22) 


The  parameters  6,,  c,  and  tf,  in  Eq.  (22)  can  be  shown  to  be: 


2  r 

bj  -  —  (1  +  2i/)l'i  + - 

rr  52 


2  1  , 

b2  = —aff?  (1  —  2t^)ki  + — kj  cos  2^ 
L  52  . 

6,  =  1  L  -  Jt, 

^  L  S2  J 

2  r  z^  ] 

64  =  —agl  (1  —  2i')ki  +  — kj  sin2ci 

T  02 


65  =  — zg2gjkiSin4> 

TT 

2 

66  =  — zg2g2k4COs4> 

TT 


Cl  =  fs  cos  4> 

Cl  =  ft  cos  +  /s  cos  ‘i4>  +  /6(C0S  </>  +  cos  3<^) 
C3  =  /4  sin  <i>  +  fs  sin  3<^  +  /6(sin  3<f>  +  sin  <^) 
C4  =  —2/6  cos 
cj  =  /?  +  h  cos  2(;6 
C6  =  /g  sin  20 


di  =  /3  sin 

d2  =  —  ft  sin  (f>  +  fs  sin  34>  +  /6(sin  34>  —  sin  <f>) 
di  =  U  cos  d>-  fs  cos  3<f>  +  /6(cos  <!>  -  cos  3<^) 
^4  =  -2/6  sin  (;!• 
ds  =  /g  sin  20 

<^6  =  /7-/8COs20 
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where 


il 

95 


/2  =  7 

95 


h  - 

/4  = 


8 


?r(2  -  v) 
8(1  -  V) 


7r(2  -  I/) 

=  K*" 


[-2(1  +  v)agjfi  +  zg295^4] 

053/1 


/8  = 


ir{2  -  u) 


glit/af2  +  zg\h) 


9^  =  (o^  - 

52  =  m  - 


a 
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“  ‘2  M 


k\  -  0/2  -  sin  *(5^4) 

/f +  aW  +  222_  3p2) 


^2  = 


^/5 


A'3  = 


a2(6/|-2/?  +  p2)_5/« 


55 


*4  = 


_£2(4^|_V^ 


55 


(65) 


(66) 


(67) 
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ni.9.  Figure  captions 


Figure  1. 

Figure  2. 

Figure  3. 

Figure  4. 
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Figure  5. 


Coordinate  systems  for  two  microcracks. 

Comparison  of  overall  longitudinal  normal  compliance  (533). 
Comparison  of  overall  longitudinal  shear  compliance  {Sss). 
Comparison  of  normalized  longitudinal  Young’s  modulus. 
Comparison  of  normalized  longitudinal  shear  modulus. 
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PART  IV 

An  Improved  Two-Dimensional  Micromechanical  Theory  for 
Brittle  Solids  with  Randomly  Located  Interacting  Microcracks 


rV.O.  Abstract 

This  paper  presents  applications  of  accurate  orthogonal  /unction  approximation  methods  to 

the  two-dimensional  problems  of  brittle  solids  containing  randomly  located  and  interacting  mi- 
% 

crocracks  within  the  framework  of  rigorous  micromechanics  and  (probabilistic)  ensemble  average 
approach.  The  random  two-crack  interaction  problems  axe  solved  by  the  highly  accurate  Legen¬ 
dre  and  Tchebycheff  OTihogor\z!i  polynomials  to  any  desired  orders.  The  complex  stress  potential 
method  is  subsequently  employed  to  micromechanically  derive  microcrack  opening  displacements 
^  under  complex  loadings  due  to  microcrack  interaction  effects  —  including  concentrated  loadings, 

arbitrary  loadings  and  polynomial  loadings.  Improved  local  ensemble-averaged  and  overall  effec¬ 
tive  elastic  compliances  due  to  microcracks  and  their  interactions  are  systematically  constructed 
by  using  the  pairwise  microcrack  interaction  mechanism  and  the  ensemble  average  approach.  A 
•  number  of  interesting  analytical-numerical  examples  containing  different  random  microcrack  con¬ 

figurations  are  also  presented  to  illustrate  the  capabilities  of  the  proposed  framework. 
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IV.l.  Introduction 


The  effects  of  microcrack  interactions  on  overall  elastic  moduli  of  brittle  solids  have  been  in¬ 
vestigated  by  many  researchers  in  the  solid  mechanics  community.  Among  others,  Ju  and  Chen 
(1991a,b)  and  Ju  and  Tseng  (1992a)  proposed  two-dimensional  and  three-dimensional  statistical 
micromechanical  damage  models  to  estimate  effective  elastic  moduli  of  brittle  solids  containing 
many  randomly  dispersed  and  interacting  microcracks.  In  particular,  statistical  aspects  and  micro¬ 
crack  interactions  are  explicitly  accounted  for.  Based  on  a  simplified  stress  interaction  fonnulation, 
analytical  closed-form  solutions  of  microcrack  interactions  are  derived  by  Ju  and  Chen  (1991a)  and 
Ju  and  Tseng  (1992a)  at  the  expense  of  some  accuracy  in  local  stresses.  Gross  (1982)  and  Horii 
and  Nemat-Nasser  (1985)  proposed  some  asymptotic  models  for  microcrack  interaction  problems. 
Embedded  in  their  models  is  the  restriction  that  microcracks  are  not  too  close  to  one  another.  On 
the  other  hand,  accurate  numerical,  deterministic  stress  interaction  approximations  were  introduced 
by  Kachanov  (1987)  and  Benveniste  et  al.  (1989).  Explicit  analytical  forms  arc  not  available  in 
their  models;  however,  higher  accuracy  in  deterministic  local  stresses  is  achieved  and  stronger 
microcrack  interactions  are  made  possible. 

In  this  study,  numerical  orthogonal  function  approximations  are  employed  to  improve  the  accu¬ 
racy  of  local  stresses  of  microcrack  interactions  and  therefore  the  accuracy  of  the  overall  ensemble 
averaged  effective  elastic  moduli.  In  Section  2  of  this  paper,  orthogonal  function  approximations 
are  applied  to  solve  the  arbitrary  (probabilistic)  local  two-microcrack  interaction  problems.  In 
particular,  Legendre  and  Tchebycheff  polynomial  families  are  implemented  to  assess  their  effec¬ 
tiveness.  Other  approximation  methods  are  discussed  and  compared  with  the  orthogonal  function 
approximation  method  in  Section  4.  Applications  of  complex  stress  potentials  on  microcrack  inter¬ 
action  problems  are  discussed  in  Section  4.  Stress  potentials  for  concentrated  loading  conditions  are 
integrated  for  arbitrary  loading  cases.  Microcrack  opening  displacements  and  strains  for  polyno¬ 
mial  loadings  of  any  arbitrary  order  are  also  presented.  In  Section  5,  overall  effective  elastic  moduli 
arc  derived  by  taking  the  ensemble-volume  averages  in  the  probability  space  over  all  possible  mi¬ 
crocrack  realizations.  Finally,  some  numerical  examples  are  presented  to  illustrate  the  proposed 
theory  in  Section  6. 


TV  2.  Orthogonal  function  approximations  of  two-crack  interaction  problems 
rv.2.1.  Two-crack  interaction  problems 

The  problem  of  two  interacting  microcracks  embedded  in  an  elastic  solid  loaded  by  far  field 
stress  <r°°  has  been  studied  by  various  methods  in  the  literature.  In  particular,  a-°°  =  (cr^., 
in  two-dimension;  see  Figure  1  for  a  schematic  plot  of  coordinate  systems. 

Let  p  and  q  denote  the  local  normal  and  shear  tractions  on  the  surface  of  a  microcrack,  respec¬ 
tively.  system  of  equations  for  the  tractions  on  two  arbitrary,  interacting  microcracks  can  be  set 
up  by  decomposing  the  original  problem  into  three  subproblems.  The  boundary  conditions  that 
miciocrack  surfaces  are  traction-free  lead  to 


=  pT ;  gi(a-i)  =  9r 
Pii^i)  =  pT  '■>  ^2(^2)  =  9r  ^(^2) 


(1) 


where  X]  and  12  are  local  horizontal  coordinates  along  the  microcrack  lines.  Further,  subscripts  1 
and  2  denote  microcrack  1  and  microcrack  2,  respectively.  Local  tractions  due  to  far  field  loading 
only  are  p°®  and  q’^.  Traction  perturbation  due  to  inter-crack  iraeractions  only  are  denoted  by 
P2^  92)-  Assuming  that  a\  and  02  are  the  half-lengths  of  the  two  microcracks,  respectively, 
Eq.  (1)  are  valid  for  -ay  <  xy  <  ay  and  -02  <  12  <  °2-  Following  Kachanov  (1987),  Benveniste 
et  al.  (1989)  and  Ju  and  Chen  (1991a),  it  is  assuif  .ed  further  that  microcracks  do  not  intersect  each 
other  arbitrarily. 


Following  Benveniste  et  al.  (1989),  we  let  Pnix),  n  =  0, 1 . 2, ....  N  denote  a  family  of  orthogo¬ 
nal  functions  defined  in  {- 1 , 1  ] .  The  tractions  pi  (i  1 ),  91  (1 1 ),  P2('r2),  and  92(12)  can  be  approximated 
by  the  base  function  set  as: 


P2M  =  ;  <nM  -  '!;p. 


Accordingly,  the  perturbed  tractions,  pi(x|),  91(^1),  ^(12),  and  %(x2)  are  linear  combinations  of 


perturbations  on  one  crack  -.vhen  the  other  crack  is  loaded  by  the  base  functions: 


N 


N 


n=0  n=0 

9i(ii) = E  + E  9? 

n^O  n=0 

km = EprF^<">(i,) + ^  ?rF="”>(x2) 

n=0  n=0 

92(1!)  =  +  E  9rFi'‘"'(i2) 


ns^ 


nj=0 


where 


=  normal  traction  on  crack  a  due  to  normal  loading  crack  jd 

=  normal  traction  on  crack  a  due  to  shear  loading  P„  ^ 

=  shear  traction  on  crack  a  due  to  normal  loading  Pn  crack  jd 

=  shear  traction  on  crack  q  due  to  shear  loading  Pn  (  “  1  on  crack  0 


Substituting  Eq.  (2)  and  Eq.  (3)  into  Eq.  (1),  we  obtain 


^  kfJ  ^  =  pr + Ef5^» + E  9?C“’(-'i) 

E 9r + E  + E  9;^"  >(x,) 

Ef?^"  (r)  +E'’r'tr’‘^^>^E9rC”W 

n=0  \02/ 

E  9?^«  f  ?)  '  9r  +  E  +  x;  9?F,T''’<^j) 


(3) 


(4) 


(5) 


The  property  of  orthogonality  of  P„  is  then  applied  to  arrive  at  the  following  system  of  cqua- 


tions: 


forn=0.1,2,.,.,A' 

(^”  fe)  (^)).’"  =  G“  ""  (^))/E 

92  forn=0,l,2,...,N 

("■  (I)  (1)).,^” = 

?r  forn=0,l,2,...,A- 

('■-fi)'-(S)).«'-(<.-'-(f;)).,*§{':*‘“'-©)/> 

^(fi'rW.P.  9?  fort.  =0,1, 2 . A- 


in  which  (f,g)u..  denotes  the  weighted  norm  of  /  and  g. 


All  the  norms  in  Eq.  (6)  arc  simple  constants  after  carrying  out  the  integrations.  The  system 
of  4(A'  +  1)  equations  is  sufficient  to  solve  for  the  4(A'  +  1)  unknowns,  p",  9",  pj,  and  g^.  In  the 
case  of  N  =  4,  there  are  20  equations  in  terms  of  20  unknowns.  For  more  accurate  solutions,  the 
number  N  can  be  increased. 


The  total  tractions  accounting  for  microcrack  interactions  are  obtained  by  substituting  the 
solved  unknowns  back  into  Eq.  (2).  However,  another  choice  exists  as  we  look  at  Eq.  (3).  That  is, 
in  addition  to  the  left-hand  side  of  Eq.  (2),  the  right-hand  side  provides  another  (better)  approxima¬ 
tion  [Benvenistc  et  al.  (1989)].  The  difference  between  these  two  choices  is  that,  with  the  former 
choice,  tractions  on  each  crack  are  solved  directly  from  the  system  of  equations  while  for  the  latter 
tractions  on  one  crack  arc  the  sum  of  tractions  due  to  the  far  field  applied  loadings  and  the  pertur¬ 
bations  induced  by  the  interaction  with  the  other  crack.  This  stress  perturbation  is  the  amount  of 
tractions  on  one  crack  while  the  other  crack  is  loaded  by  the  total  tractions.  In  another  words,  the 
total  tractions  solved  from  Eq.  (6)  are  projected  onto  the  other  crack.  This  process  is  embedded 
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in  Eq.  (3).  Stress  projections  are  equivalent  to  the  cvaluat’O^*  of  those  inrtucnce  functions.  F^, 

rP 

Ca0  pa0  ar>Apo0 
^P<t  '  ^<IP  ’  ^qq  • 

IV.2.2.  Legendre  and  Tchebycheff  polynomials 

Among  the  orthogonal  function  sets,  Legendre  and  Tchebycheff  polynomials  arc  widely  used 
in  various  applications.  Both  polynomials  are  tested  in  this  study.  The  definitions  of  Legendre 
Polynomials  are 

Loix)  =  0,  L„(x)  =  :i^[(x2-in  O) 

2'*n!  ax" 

with  the  following  orthogonal  property 


lo, 


forn=j 


otherwise 


The  weighting  function  for  Legendre  polynomials  is  wl(x)  -  1.  The  system  of  equations  for  two- 
crack  interaction  with  the  application  of  Legendre  polynomials  is  as  follows 

- 1  (^))  ri- 1  (^)) 

=  2a,pf“(5„o,  forn  =0,1,2,. ..,Af 

- 1  (f ))  p;- 1  (fl)) 

=  2ai9j"<5„o,  for  n  =0, 1,2,...,  A' 

-  t  in  (f ))  K  -  t  in  (f )) 

-lazpfSno,  forn  =0, 1,2, ...,  A 

-i  (fi))  Pi  -  t  (F-<n,,.,,,i,  (^))  ,n 


=  2a29^<5„o,  forn  =  0, 1,2,...,A 


where  Sij  =  0  for  i  =  j  and  0,  otherwise. 


On  the  other  hand,  the  triangle  family  of  orthogonal  polynomials,  Tchebycheff  polynomials,  is 


defined  as 


with  the  orthogonal  propeny 


7’o(x)=l,  r,(x)  =  x 

Tn*l(x)  =  2ir„(i)  -  T„_i(x) 


{TT,  for7=j=0 
f,  fort  =;■  7^0 
0,  otherwise 
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The  weighting  function  for  Tebycheff  family  is  wt(x)  =  (1  -  If  t^’'  family  is  applied,  we 

have  the  following  equations  for  the  two-crack  interaction  problem; 


TTQi 


2  -  ^nO 


TTOi 

2^n0 


naO 

=  TToipf’^nOifor  n  =0, 1,2,  ...,N 


92 


WJ* 


WT  nsO 

=  5raigj"6„o»for  n  =  0, 1,2,  ...,N 


=  naipfSno.foT  n  =  0, 1,2,  ...,N 


7ra2 


2-6, 


•nO 


-92 


=  fffi29r^n0,  for  n  =  0, 1,2, ..., 


(12) 


Due  to  the  fact  that  base  functions  of  Legendre  andTchebychcff  approximations  are  polynomi¬ 
als,  if  no  stress  projection  is  performed,  the  tractions  obtained  by  solving  either  Eq.  (9)  or  Eq.  (12) 
and  substituting  the  coefficients  into  Eq.  (2)  can  be  rearranged  by  collecting  the  terms  of  the  same 
power  of  I. 


Pi(i,)=pr+E’"r(^)'  ;  9.(x,)=?r+E2?(r) 

n=0  '  *  '  »»=0  ■' 

P2(X2)=Pr+X]^2  i  92(j^2)  =  9r  +  5^22 

«=n  n=0  ^ 


(13) 


Tractions  with  the  operations  of  stress  projection  are  similarly  obtained  by  the  right-hand  side  of 
Eq.  (5). 

N  N 

pi(xi) = pr + E  + E 


n=0 

N 


naO 

N 


qiixi)  =  qr  +  Yi  ^ 

n=0 

AT  N 

P2(X2) = ft” + E  ^r^S‘"’<x2) + E  "’(^2) 

n=0  naO 

N  N 

®(X2)  =  ft”  +  ?’r-^i',''"’(x2) + XI  e?;r;,“"’(i2) 


(14) 


n=0 


n=0 
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in  which  influence  functions  are  similar  to  those  defined  in  Eq.  (4) 


_  normal  traction  on  crack  a  due  to  normal  loading  ^  crack  0 

_  normal  traction  on  crack  a  due  to  shear  loading  ^  crack 

=  shear  traction  on  crack  a  due  to  normal  loading  ^  crack  0 

=  shear  traction  on  crack  o  due  to  shear  loading  crack  0 


(15) 


IV.23.  Comparisons  of  Legendre  and  Tchebycheff  approximations 

Both  Legendre  and  Tchebycheff  approximations  employ  polynomials  as  the  base  functions. 
They  are  both  orthogonal  and  have  the  same  symmetry  property: 


Z,„(a-)  =  (-l)"P„(-x)  ;  r„(x)  =  (-l)"r„(-x)  (16) 

The  leading  orders  of  the  n*^  functions,  I„(x)  and  T„ix),  are  both  equal  to  n.  TTie  shapes  of  base 
functions  of  the  same  order  are  quite  similar.  Weighting  functions  are,  peibaps,  the  most  noticeable 
difference.  For  the  Legendre  family,  the  weighting  function  is  simply  1.  The  weighting  function 
for  the  Tchebycheff  family  is  singular  at  the  two  crack  tips. 


To  compare  the  effectiveness  of  two  approximation  function  sets  on  the  problem  of  two-crack 
interaction,  let  us  consider  the  stress  intensity  factors  (SIF)  at  crack  tips.  For  Mode  I  loadings, 
•  SIFs  at  both  crack  tips  can  be  calculated  by  integrating  over  the  crack  line  on  the  product  of  total 

traction  p(x)  and  the  SIF  for  a  unit  concentrated  normal  loading  [Tada  et  al.  (1973)): 


flit 
a  ^  t 


p{t)dt 


(17) 


For  the  cases  in  which  p(x)  are  simple  polynomials,  analytical  SEF  can  be  obtained.  Assuming  that 
pix)  =  (4)",  it  can  be  shown  that  [by  carrying  out  the  integration  in  Eq.  (17)) 


Kf‘"Uc”v^  (18) 

where  c"  is  a  simple  constant  which  can  be  easily  calculated  by  symbolic  mathematical  softwares 
like  MAPLE,  MATHEMATICA,  or  MACSYMA  for  any  number  of  n.  Manuals  for  these  softwares 
are  provided  in  the  reference.  Table  1  shows  c”  for  n  ranging  from  0  to  8.  The  authors  like  to  point 
out  that  the  applications  of  symbolic  computation  softwares  in  various  research  fields  are  becoming 
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Table  1.  Constants  for  SIFs  of  polynomial  loadings 


a 

0 

1 

2 

3 

4 

5 

6 

7 

■ 

1 

la 

B 

n|«o 

-H 

3 

8 

IBI 

1^1 

popular  as  the  capability  of  computers  improved  by  the  modem  technology.  loaldmidis  (1990) 
applied  these  power  tools  to  the  crack  problems  in  fracture  mechanics.  More  references  about  the 
so-called  semi-analytical/numerical  (SAN)  method  can  be  found  in  loakimidis  (1992)  in  which 
the  numerical  and  symbolic  computations  are  combined  to  solve  the  singular  integral  equations 
commcmly  found  in  fracture  mechanics. 

These  simple  results  are  applicable  to  the  problems  of  solving  crack  dp  SIFs  for  two-crack 
interaction  without  stress  projection.  Combining  Eq.  (13),  Eq.  (17),  and  Eq.  (18),  SIFs  of  the  two 
micTOcracks  can  be  expressed  by  simple  summations  on  the  coefficients  solved  from  the  system  of 
equations 

N  N 

=  +  ’  Kf/^>  =  K^(l  +  ^c"'P2")  (19) 

n=0  n=0 

where  the  SIFs  for  non-interacting  case  are  defined  by 

Kioj  =  pf*  y/iFa  ;  =  pf  y/iFa  (20) 


Table  2.  Comparison  of  normalized  SIFs  for  two  interacting  collinear  cracks 


Inner  tip 

Outer  tip 

k 

Exact 

Legendre 

Tchebycheff 

Exact 

Legendre 

Tchebycheff 

0.9 

1.0004 

1.0004 

1.0004 

1.0003 

1.0003 

1.0003 

0.75 

1.0028 

1.0028 

1.0028 

1.0024 

1.0024 

1.0024 

0.5 

1.0176 

1.0176 

1.0176 

1.0125 

1.0125 

1.0125 

0.25 

1.0804 

1.0804 

1.0804 

1.0409 

1.0410 

1.0409 

0.2 

1.1125 

1.1124 

1.1125 

1.0517 

1.0517 

1.0517 

0.1 

1.2551 

1.2539 

1.2551 

1.0863 

1.0870 

1.0863 

0.05 

1.4729 

1.4650 

1.4733 

1.1198 

1.1229 

1.1198 

0.02 

1.9046 

1.8593 

1.9132 

1.1589 

1.1727 

1.1595 

0.01 

2.3716 

2.2475 

2.4201 

1.1841 

1.2149 

1.1864 

0.001 

5.3947 

4.0501 

10.6993 

1.2443 

1.4029 

1.3654 

0.0001 

13.3456 

6.0830 

-4.7501 

1.2808 

1.6190 

1.0973 

10-6 

93.0293 

7.9860 

-1.2083 

1.3212 

1.8234 

1.1646 
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Tables  2,  3,  and  4  show  the  comparisons  of  crack  tip  SIFs  for  the  three  cases  depicted  in 
Figure  2.  Since  the  two  microcracks  are  parallel  to  each  other  in  all  the  three  cases,  SIFs  of  both 
microcracks  are  identical.  The  non-interacting  SIF,  K^,,  is  used  to  normalize  all  SIFs  in  the  three 
tables.  Normalized  SIFs  in  Table  2  are  always  greater  than  one.  This  is  in  agreement  with  the 
physical  observation  since  two  collinear  cracks  intuitively  enhance  the  tractions  for  each  other.  On 
the  contrary,  two  stacked  cracks  act  against  each  other  from  opening  the  crack  lines  and  therefore 
SIFs  given  in  Table  3  are  all  less  than  one. 


Table  3.  Comparison  of  normalized  SIFs  for  two  stacked  cracks 


rja 

Inner/Outer  tip 

Legendre 

Tchebycheff 

100 

1.000 

1.000 

10 

0.986 

0.986 

5 

0.951 

0.951 

2 

0.843 

0.843 

1 

0.774 

0.788 

0.4 

0.772 

0.721 

0.2 

0.677 

0.696 

0.1 

0.628 

0.674 

0.05 

0.588 

0.654 

0.01 

0.535 

0.620 

0.005 

0.524 

0.610 

0.001 

0.506 

0.596 

Note  that  SIFs  at  the  two  crack  tips  for  stacked  cases  are  identical.  Both  Legendre  and  Tcheby- 
cheff  family  lead  to  very  close  results  in  the  stacked  cases  even  when  the  distance  between  the  two 
crack  is  very  small.  For  the  third  (aligned)  case  in  Table  4,  significant  differences  are  observed  only 
at  the  inner  tips  when  the  two  cracks  are  very  close  to  each  other.  Exact  solutions  given  by  Erdogan 
(1962)  are  included  in  Table  2  for  comparison  for  the  first  case  in  which  two  cracks  are  collinear. 
Table  1  shows  that  if  the  cracks  are  very  close,  Legendre  family  still  works  well  while  Tcheby- 
cheff  family  results  in  negative  SIFs  at  the  inner  tips  which  are  actually  unreasonable.  Both  sets 
yield  nearly  identical  results  if  two  cracks  are  far  away  from  each  other.  It  is  interesting  to  observe 
that  except  those  unreasonable  results,  TchebychefF  method  tends  to  over-estimate  the  interaction 
effect.  Consequently,  SIFs  at  the  inner  tips  calculated  through  Tchebycheff  method  are  larger  than 
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the  exact  values.  On  the  contrary,  the  interaction  effects  are  under-estimated  by  the  Legendre  ap¬ 
proximation  and  the  inner-tip  SIFs  are  less  than  the  exact  ones.  Table  2  also  shows  that  if  the  tv,o 
collinear  crack  are  not  too  close,  k  >  0.01  for  example,  Tchebycheff  method  appears  to  be  a  bener 
approximation  to  the  exact  solution  than  the  Legendre  method  though  the  approximation  by  the 
later  is  actually  quite  close  to  the  exact  result.  However,  in  the  cases  which  two  collinear  cracks 
are  extremely  close,  Legendre  method  seems  to  be  a  better  choice  since  the  Tchebycheff  method 
behaves  unreasonably.  One  possible  reason  for  the  poor  approximation  of  Tchebycheff  family  at 
extremely  close  distance  is  the  singularities  at  both  crack  tips.  In  the  collinear  cases,  tractions  at  the 
inner  crack  tips  rise  sharply  as  the  distance  between  the  two  cracks  decreases.  This  fast  growth  is 
amplified  by  the  weighting  function  which  is  singular  at  crack  tips.  The  comparisons  above  show 
that  Legendre  family  is  a  more  robust  choice  for  the  crack  problems  especially  when  cracks  are 
very  close  and  collinear. 


Table  4.  Comparison  of  normalized  SIFs  for  two  aligned  cracks 


r/a 

Inner  tip 

Outer  tip 

Legendre 

Tchebycheff 

Legendre 

Tchebycheff 

100 

1.000 

1.000 

1.000 

1.000 

10 

1.004 

1.004 

1.005 

1.005 

5 

1.010 

1.010 

1.020 

1.020 

2 

0.847 

0.847 

1.078 

1.077 

1 

0.481 

0.478 

1.087 

1.091 

0.4 

0.234 

0.235 

1.050 

1.044 

0.2 

0.180 

0.183 

1.016 

1.011 

0.1 

0.152 

0.168 

0.986 

0.990 

0.05 

0.135 

0.162 

0.964 

0.977 

0.01 

0.115 

0.158 

0.936 

0.962 

0.005 

0.113 

0.160 

0.931 

0.959 

0.001 

0.126 

0.162 

0.932 

0.955 

Table  5  compares  the  crack-tip  SIFs  calculated  by  fourth  order  Legendre  approximation  method 
with  and  without  the  stress  projection.  The  two  interacting  cracks  are  collinear  as  in  the  first  case 
of  Figure  2.  The  numerical  experiments  show  that  the  stress  projection  is  not  important  if  the  two 
cracks  are  not  very  close.  In  Table  5,  SIFs  with  and  without  stress  projection  are  quite  close  for  k 
greater  or  equal  to  0.01 .  However,  when  the  two  cracks  arc  extremely  close,  the  operation  of  stress 


projection  appears  to  be  necessary  if  the  accuracy  is  concerned.  The  approximated  SIFs  are  much 
improved  toward  the  exact  solution  for  the  cases  o(  k  -  O.CX) 1,0.0001,  and  10~*  in  which  cracks 
are  extremely  close. 


Table  5.  Comparison  of  normalized  SIFs  for  two  interacting  coUinear  cracks 
by  the  fourth  order  Legendre  approximation  with  and  without  stress  projection. 


k 

Inner  tip 

Outer  tip 

Exact 

Solution 

No 

Projection 

Exact 

Projection 

Exact 

Solution 

No 

Projection 

Exact 

Projection 

0.9 

1.0004 

1.0004 

1.0004 

1.0003 

1.0003 

1,0003 

1.0028 

1.0028 

1.0028 

1.0024 

1.0024 

1.0024 

0.5 

1.0176 

1.0176 

1.0176 

1.0125 

1.0125 

1.0125 

0.25 

1.0804 

1.0804 

1.0804 

1.0409 

1.0410 

1.0409 

0.2 

1.1125 

1.1124 

1.1125 

1.0517 

1.0517 

1.0517 

0.1 

1.2551 

1.2539 

1.2551 

1.0863 

1.0870 

1.0863 

0.05 

1.4729 

1.4650 

1.4733 

1.1198 

1.1229 

1.1198 

0.02 

1.9046 

1.8593 

1.9093 

1.1589 

1.1727 

1.1595 

0.01 

2.3716 

2.2475 

2.3899 

1.1841 

1.2149 

1.1858 

0.001 

4.0501 

5.5333 

1.2443 

1.4029 

1.2614 

0.0001 

13.3456 

6.0830 

11.6430 

1.2808 

1.6190 

1.3244 

lo-*^ 

93.0293 

7.9860 

26.3239 

1.3212 

1.8234 

1.3785 
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IV.3.  Other  approximation  methods 


IV.3.1.  The  zeroth-order  approximations 


Several  approximation  techniques  have  been  proposed  in  the  literature  to  solve  Eq.  (1).  For 
example,  Ju  and  Chen  (1991),for  simplicity,  expanded  the  traction  functions  in  terms  of  polyno¬ 
mials  with  respect  to  the  center  of  microcracks,  and  (for  computational  simplicity  only)  took  only 
the  zeroth-order  terms  which  were  constants  for  specified  microcrack  locations  and  orientations. 
Therefore,  Eq.  (1)  can  be  recast  as 


Pi  =Pr+Pi 
qi=qT‘+qi 
P2  =  pT+P2 


(21) 


qi=qT  ■‘■ft 

By  employing  Sneddon  and  Lowengrub  (1969)  solutions,  Ju  and  Chen  (1991)  obtained  approxi¬ 
mate,  analytical  solutions 


Pi ) 

r  0 

0 

Ol 

02" 

r  pi ' 

qi  1 

P2  1 

>  = 

0 

«5 

0 

03 

0 

04 

0 

< 

ft 

P2 

ft  i 

.07 

08 

0 

0  . 

1  92 

Definitions  for  a,  ’s  can  be  found  in  the  Appendix  of  Ju  and  Chen  (1991). 


(22) 


For  compactness,  let  us  define  Ti_2  =  (pi,ft,P2,92)^.  —  (pr»9r»Pri9r)^»*‘”dTi-2  == 

(Pi>  ftj  ft)^‘  Therefore,  Eq.  (21)  and  Eq.  (22)  arc  rewritten  in  matrix  and  vectorial  forms 

Ti_2  =  a  •  Ti_2 

where  a  denotes  the  4  by  4  matrix  in  Eq.  (22).  The  perturbed  tractions  due  to  two-microcrack 
interaction  are  solved  from  Eq.  (23) 


f,_2  =  aa-a)-‘-Tr_2 


(24) 


Definition  of  Ki_2  =  a  •  (I  —  a)-‘  and  transformation  of  the  far  field  loading  into  the  local 
microcrack  coordinates  then  render 


where 


Ti_2  =  Ki_2  •  Gi_2  • 


Gi_2  = 


sin^  9 
—  4  sin  29 

sin*(^  +  <t>) 
.-|sin2(6>4-^) 


co^9 
5  sin  29 

cosH9  +  <f>) 
5  sin  2{9  +  <l>) 


—  sin  29 
cos  29 

-  sin  2(9  +  <f>) 
cos  2(9  +  (f>)  . 


(25) 


(26) 
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Moreover,  we  define  Ti  =  (/).  qxf  and  f  2  =  (p2,  Therefore,  Eq.  (25)  can  be  rephrased  for 
two  microcracks  as  follows 

t,  =K,-G,_2  0-~  (27) 

t2  =  K2  •  Gi_2  •  (28) 

where  Ki  contains  the  first  two  rows  of  K 1-2  and  K2  the  last  two  rows. 

In  addition  to  the  aforementioned  zeroth-order  polynomial  expansion  at  microcrack  centers, 
Kachanov  (1985,  1987)  proposed  a  more  accurate  approximation  which  also  led  to  a  system  of 
equation;  in  terms  of  some  constants.  Eq.  (1)  are  first  averaged  over  the  microcrack  lines.  In 
essence,  traction  perturbations  on  one  microcrack  due  to  the  interaction  with  another  microcrack 
loaded  by  an  averaged  traction  are  averaged  over  the  microcrack  line.  Accordingly,  the  following 
system  of  equations  can  be  obtained  (Kachanov,  1985) 

Pi  =  pT  ^12  Pi  ^12  92 

91=9“  +  p2  +  92  ^29) 

P2  =  Pf  +  A21  Pi  +  ^21  9i 
92  =  92“  +  A2?Pi+A5Ui 

For  a  given  deterministic  microcrack  configuration,  transmission  factors  Afj  must  be  computed 
by  numerical  integrations.  With  pre-computed  transmission  factors,  Eq.  (29)  become  a  system  of 
simultaneous  equations  in  terms  of  constants  only,  which  can  be  solved  for  the  unknown  average 
tractions.  Note  that  the  matrix  operation  in  Eq.  (24)  serves  the  same  purpose  as  the  solutions  of 
simultaneous  equations  in  Eq.  (29).  Certainly,  Eq.  (29)  are  more  accurate  than  (24).  However, 
much  less  computational  efforts  are  required  for  the  analytical  zeroth-order  polynomial  expansion 
since  no  transmission  factors  involving  numerical  integrations  are  needed.  Further,  Eq.  (24)  are 
expressed  in  close-form.  The  efficiency  of  analytical  expressions  over  numerical  computations 
might  not  be  very  significant  for  a  very  few  deterministic  microcracks.  Nevertheless,  for  ensemble- 
averaged  effective  compliance  computations  of  statistical  microcrack  arrays  (involving  hundreds 
of  random  microcracks),  the  computational  efficiency  and  simplicity  of  the  zeroth-order  polynomial 
expansion  are  more  favorable. 

rv.3.2.  Other  higher  order  approximations 

There  are  a  number  of  methods  in  the  literature  which  involve  higher-order  approximations  of 
microcrack  interaction  problems.  Emanating  from  Eq.  (1),  the  traction  funedons  were  expanded 


into  Fourier  series  in  Horii  and  Nemat-Nasser  ( 1 985),  Moreover,  Gross  ( 1 982)  approximated  trac* 
tion  functions  by  Tchebycheff  polynomials.  Both  methods  involve  choosing  a  parameter  which  rep¬ 
resents  the  ratio  of  the  distances  between  microcrack  centers  to  microcrack  lengths.  Higher-order 
terms  involving  a  chosen  parameter  are  asymptotically  dropped  such  that  the  number  of  terms  cor¬ 
responding  to  expanded  tractions  becomes  finite.  A  system  of  equations  containing  the  unknown 
coefficients  of  the  expansion  and  known  constants  are  then  constructed.  Approximate  tractions  are 
computed  by  substituting  those  coefficients  solved  from  simultaneous  equations  back  into  the  trac¬ 
tion  expansions.  Significant  errors  are  observed  for  the  foregoing  two  methods  when  microcracks 
are  closely  located. 

IV.33.  Comparisons  of  SIFs  for  zeroth  order  approximation 

For  the  case  of  two  collinear  interacting  cracks  in  Figure  2,  the  SIFs  solved  by  applying  the 
zeroth  order  approximation  at  center  points  are  listed  in  Table  6  and  Table  7  for  the  inner  and 
outer  tips,  respectively.  Results  by  Kachanov  and  Gross  are  also  compared  with  the  exact  values 
given  by  Erdogan  (1962).  Kachanov’s  method  which  is  a  special  case  of  Benveniste’s  method  ap¬ 
pears  to  be  the  best  approximation  in  both  tables.  If  no  projection  is  performed  in  the  center-point 
approximation,  the  results  are  no  better  than  Gross’.  However,  with  a  simple  stress  projection, 
the  center-point  method  becomes  better  than  the  fourth  order  asymptotic  method  of  Gross,  Actu¬ 
ally,  center-point  zeroth  order  approximation  with  stress  projections  is  comparable  to  the  method 
proposed  by  Kachanov. 

Although  center-point  zeroth  order  approximation  method  is  improved  by  the  operation  of 
stress  projection,  the  orthogonal  function  approximation  method  employed  in  this  paper  is  shown 
to  be  better  when  Table  6  and  Table  7  are  compared  with  Table  2.  Another  advantage  for  the 
orthogonal  function  approximation  method  is  that  higher  accuracy  may  be  achieved  by  simply 
increasing  the  number  of  order  used  in  the  numerical  implementation.  The  entire  formulation  is 
highly  systematic  so  that  the  computer  implementation  for  any  order  of  approximation  is  easy  and 
straightforward. 


121 


Table  6.  Comparison  of  normalized  inner-tip  SIFs  for  two  interacting  collinear  cracks 


k 

Exact 

Kachanov 

Gross 

Zeroth  Order  Approximation 

No  Projection 

Exact  Projection 

0.9 

1.0004 

1.0004 

1.0004 

1.0004 

0.75 

1.0028 

1.0028 

1.0026 

1.0028 

0.5 

1.0176 

1.0173 

1.0144 

1.0176 

0.25 

1.0804 

1.0696 

1.0507 

1.0797 

0.2 

1.1125 

1.112 

1.0911 

1.0646 

1.1109 

0.1 

1.2551 

1.251 

1.1564 

1.1061 

1.2423 

0.05 

1.4729 

1.452 

1.2062 

1.1380 

1.4200 

0.02 

1.9046 

1.809 

1.2441 

1.1630 

1.7039 

0.01 

2.3716 

2.134” 

1.2583 

1.1727 

1.9410 

0.001 

5.3947 

3.401* 

1.2719 

1.1819 

2.7825 

0.0001 

13.3456 

4.731* 

1.2733 

1.1829 

3.6462 

10-^ 

93.0293 

7.309* 

1.2734 

1.1830 

5.3798 

*  from  Benveniste  et  al.  (1989)  **  2.138  from  Benveniste  et  al.  (1989) 


IV.4.  Applications  of  complex  stress  potentials  to  crack  interaction  problems 


IV.4.1.  Concentrated  loadings 


Consider  the  problem  of  a  single  crack  with  half-length  a  loaded  at  x  =  6  by  normal  and 
shear  forces,  P  and  Q,  as  shown  in  Figure  3.  Complex  stress  potentials  which  solve  this  problem 
completely  are  readily  given  in  Tada  et  al.(1973): 


/  Z,(z)  1  ^  1  fP) 

I  Ziiiz) }  TT  \  Q  }  (z  —  b)\/z^  —  a? 


(30) 


(2  —  b)\/  z^  — 

where  z  =  x  +  iy  and  i  =  v/^.  For  normal  loading  (Mode  I)  the  stresses  and  displacements  are 

(7xx  -  ReZi  —  yImZ'j 

ayy  =  ReZi-¥ylmZ]  (31) 

Txy  =  -yReZ'j 


and 


IGu  =  (1  —  2v)RtZi  —  yImZi 
2Gv  =  2(1  —  v)lmZi  —  yR^Zi 


(32) 


I 


Table  7.  Comparison  of  normalized  outer-dp  SIFs  for  two  interacting  collincar  cracks 


Exact 


1.0003 


1.0024 


1.0125 


1.0409 


1.0517 


1.0863 


1.1198 


1.1589 


1.1841 


1.2443 


1.2808 


1.3212 


Kachanov 


1.052 


1.086 


1.118 


1.154 


1.175 


1.214* 


1.227* 


1.233* 


Gross 

Zeroth  Order  Approximatic 

No  Projection 

Exact  Projection 

1.0003 

1.0004 

1.0003 

1.0024 

1.0026 

1.0024 

1.0126 

1.0144 

1.0125 

1.0426 

1.0507 

1.0406 

1.0540 

1.0646 

1.0510 

1.0880 

1.1061 

1.0828 

1. 1136 

’.1380 

1.1096 

1.1332 

1.1630 

1.1339 

1.1406 

1.1727 

1.1451 

1.1476 

1.1819 

1.1591 

1.1484 

1.1829  ’ 

1.1613 

1.1484 

1.1830  ! 

1.1616 

*  from  Benveniste  et  al.  (1989) 


In  the  case  of  shear  loading  (Mode  11), 

Cxj.  -  llmZii  +  yReZ'ij 
Cyy  —  yReZji 

Tj.y  =  ReZii  -  ylmZ’ij 

and  _ 

2Gu  =  2(1  —  u)lmZn  +  yReZn 

2Gv  =  —(1  —  2i/)ReZii  —  ylmZn 

where 

[  ry  f  \  j  j  rrt  d^//(2) 


Zii  =  J  Z„{2)  dz  and  Z'lj  = 

Shear  modulus  and  Poisson’s  ratio  are  denoted  by  G  and  u,  respectively  in  Eq.  (32)  and  Eq.  (34). 
The  displacements  in  x  and  y  direction  are  u  and  w,  respectively. 


When  a  microcrack  is  loaded  by  only  normal  pressure,  the  problem  becomes  symmetric  with  re¬ 
spect  to  the  x-axis.  Relative  horizontal  opening  displacements  vanish  and  relative  vertical  opening 


displacements  are  simply  twice  the  vertical  displacements  along  the  microcrack  line.  By  contrast, 
in  the  case  of  shear  loading,  the  problem  becomes  anti-symmetric  respect  to  the  microcrack 
line.  Therefore,  the  relative  vertical  opening  displacements  vanish  and  the  relative  horizontal  open¬ 
ing  displacements  are  simply  two  times  the  horizontal  displacement  along  the  microcrack  line.  See 
Figure  4  for  a  schematic  representation. 

Accordingly,  the  effects  of  normal  and  shear  loadings  to  vertical  and  horizontal  opening  dis¬ 
placements  are  uncoupled.  Stress  potentials  arc  calculated  by  Eq.  (30).  The  displacement  com¬ 
ponents  along  the  microcrack  line  are  then  obtained  through  Eq.  (32)  or  Eq.  (34).  Therefore,  for 
normal  loading,  we  have 

llul]  =  0  and  |Iu]i  =  2t;|,=o  (?'>) 

and  for  shear  loading 

M  =  0  and  =  (36) 

Assuming  that  the  total  volume  of  the  solid  containing  the  microcrack  to  be  V’,  with  the  ap¬ 
plication  of  divergent  theorem,  the  a\  ,;rage  microcrack-induced  strain  can  be  written  as  [Horii  and 
Nemat-Nasser  (1983)] 

7=  =  ^  f  {m®n  +  n®m)dSi  (37) 

Js, 

where  n  denotes  the  normal  of  the  crack  lines  S/  and  vector  HuJ  =  (Hul,  II vj)^  the  crack  opening 
displacements.  Substituting  Eq.  (30)  into  Eq.  (32)  and  Eq.  (34)  and  by  the  nature  of  symmetry  and 
anti-symmetry,  respectively,  of  Mode  /  and  Mode  II  loadings  discussed  above  the  induced  strains 
are 

in  which  the  Young’s  modulus  is  denoted  by  E. 

IV.4.2.  Arbitrary  loadings 

For  any  arbitrary  forms  of  loadings  P(i)  or  Q(i),  the  concept  of  the  famous  Green’s  function  is 
applied  to  derive  to  the  stress  potentials  and  the  strains.  The  stress  potentials  given  by  Eq.  (30)  and 
the  corresponding  strains  given  by  Eq.  (38)  arc  for  the  case  of  a  crack  subjected  to  the  concentrated 
loading  applied  at  a  point  on  the  two-dimensional  crack  line.  When  a  normal  or  shear  loading  in 
the  form  of  an  arbitrary  function  apllied  on  the  crack  line,  the  two  equations  serve  as  the  Green’s 


function  for  this  problem.  Strains  arc  therefore  related  to  the  loadings  by  integrating  the  product  of 
loading  functions  and  Eq.  (38)  over  the  entire  crack  line 


4(1  -  1/2) 
VE 


dx 


Following  the  same  logic,  the  corresponding  stress  potentials  can  be  derived  from  Eq.  (30) 


(39) 


/  \  =  I 

/*“  y/aJ  —  x^ 

fPMX 

\  Zll(z)  J  j 

Lo  (z  —  xjv/z^  —  a* 

1  <?(!)/ 

dx 


(40) 


and  the  cenresponding  stresses  and  displacements  are  obtained  by  the  substitution  of  stress  poten¬ 
tials  calculated  fix)m  Eq.  (40)  into  Eqs.  (31)  to  (34). 

% 


rv.4  Polynomial  loadings 


Consider  the  special  cases  where  loadings  are  simple  polynomials.  Assuming  that  P(x)  = 
P n(f )"  and  Q{x)  ~  the  integration  in  Eq.  (39)  is  an  Euler  integration  of  the  first  kind  by  a 

change  of  variable.  With  the  help  of  the  Beta  function,  strains  are  written  as 


(41) 


where  the  Beta  function  is  defined  in  terms  of  the  Gamma  function  r(x): 

B{p,q)  =  f  x^-\\  ~  I)’-'  dx  =  (42) 

Jo  r(p  +  q) 

Similarly,  stress  potentials  for  polynomial  loadings  can  be  derived.  Listed  below  are  stress  poten¬ 
tials  for  the  polynomial  loadings  up  to  the  fourth  order: 


= 

Zj%)  =  Zf{z)  = 
Zf\z)  =  Zf^iz)  = 
Zf\z)  =  Zf{z)^ 
ZW(z)  =  zW(r)  = 


2  -  \/?  —  <P 

v/z^  _  q2 

-  z\/z^  -  €?■  -  dJjl 

ay/  z^  —  o2 

z^  —  z^y/z^  —  a2  —  za^jl 
aJs/z"^  —  aJ 

Z*  -  z^Vz^  -  dJ  -  z’^a^jl  -  fl*/8 
a^v/z^  —  a} 

z*  —  z*\/z^  —  —  z^aJfl  —  za*/% 

a*y/z^  -  aJ 


(43) 


Microcrack  opening  displacements  for  a  microcrack  subjected  to  polynomial  normal  and  shear 
loadings 


(44) 
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can  thus  be  derived  by  method  outlined  in  this  section; 


where  A„(x)  for  n  =  0, 1 , 2, 3, 4  for  plane  strain  are 

.  - : 


= 

Ai(x)  = 
A2(x)  = 


>l3(l)  = 


At(x)  = 


E 

4(1  — 

E 

4(1  -  t/2) 
E 

4(1  -  t/2) 
E 

4(1  —  u^) 
E 


3^  — 


\6  Sa^/ 

\  8a  4a^  J 

f-^ - r  /  3  \ 

V40  10a2  5a*  J 


The  factor  (1  -  v^)  should  be  removed  from  Eq.  (46)  for  plane  stress  problems.  It  is  noted  that  micro- 
crack  opening  displacements  corresponding  to  even-order  polynomial  loadings  arc  even  functions. 
On  the  other  hand,  the  odd-order  polynomial  loadings  result  in  odd  microcrack  opening  displace¬ 
ment  functions.  That  is, 

^n(-a:)  =  (-l)"^„(x)  (47) 

for  non-negative  integers  k.  Derivations  of  Eq.  (43)  and  Eq.  (46)  are  carried  out  with  the  assis¬ 
tance  of  symbolic  mathematical  softwares  MAPLE  and  MATHEMATICA.  Similar  results  may  be 
obtained  for  any  number  of  n. 

Now,  assume  that  there  are  N  cracks  which  arc  all  aligned  such  that  their  normal  directions 
coincide  the  j/-axis.  Further,  we  assume  that  tensile  and  shearing  forces  are  applied  homogeneously 
over  the  solids.  By  neglecting  the  interaction  effects,  the  volume  average  compliances  due  to  the 
presence  of  microcracks  can  be  derived  by  letting  n  =  0  in  Eq.  (41)  and  summing  over  N  cracks; 

jVo^4(l-t/^)  3  1 

Po  Qo  V  E 
1  - 


Note  that  the  well-known  Taylor’s  model  is  recovered  by  Eq.  (48). 
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rv.5.  Local  and  overall  effective  elastic  moduli 


IV.5.1.  Local  ensemble  averaged  moduli 

Assuming  that  the  distribution  of  microcracks  satisfies  reasonable  randtnnness  and  local  ho¬ 
mogeneity  [Hinch  (1977)]  and  that  no  microcracks  intersect  one  another,  statistical  tractions  at  a 
local  material  point  are  obtained  by  considering  the  effect  of  two-crack  interaction  and  averaging 
the  effect  over  all  possible  microcrack  configurations.  See  Ju  and  Tseng  (1992b)  for  a  more  de¬ 
tailed  discussion.  Let  T(i)  =  (p(x),  q(x))^  be  the  vector  of  local  tractions.  Ensemble  averaged 
tractions  at  a  local  point  are  the  sum  of  far  field  loading-induced  component  and  the  component 
due  to  two-crack  interacdon; 


(T)(xi,ai,  ^i)  =  {Ko(xi,ai,5i)  +  /(xiini5^i)(K)(Xi,ai,5i)}  • 


(49) 


where  the  transformarion  matrix  is  defined  by 


Ko(xi,ai,di) 

and  (in  the  probability  space) 


sin^ 

—5  sin  2^1 


cos^  61  —  sin  26\ 

5  sin  2^1  cos  2^1 


(50) 


(K)(xi,ai,0i)  =  j  KixuaiJuXz, 02,62)  <i^2<ia2d02  (51) 

in  which  K  denotes  a  two  by  three  matrix  whose  components  arc  solved  from  the  system  of  equa¬ 
tions  described  in  Section  2,  Each  column  of  K  represents  a  specific  loading  condition.  The  first, 
second,  and  third  column  correspond  to  the  cases  in  which  the  solid  is  loaded  only  by  ff",  and 

T~,  respectively.  The  two  components  in  each  column  are  the  normal  and  shear  tractions  solved 
from  the  system  of  equations.  Local  tractions  are  once  again  statistically  averaged  over  all  possi¬ 
ble  crack  configurations  at  that  material  point.  With  the  help  of  Eq.  (39),  ensemble  averaged  local 
strain  takes  the  form 


(e-)(x)  =  /  J  j"  .  (T)/(a,«)  dt  da  de  (52) 

where  g  is  local-global  transformation  matrix.  Note  that  the  ensemble  operator  is  perfwmed  in  die 
probability  space,  not  in  the  physical  space. 


Local  statistical  compliances  are  then  obtained  by  substituting  Eq.  (49)  into  Eq.  (52).  The 
microcrack-induced  compliances  are  separable  into  two  parts: 

.2(1  -  t/2) 

III  y  o"  —  T‘g-  a;  ai  aa  aa 

(53) 


(S'’)(x)  =  /(x)^^l-^  ^  ^  v/a2  -  Pg  ■  Ko/(a,  6)  di  da  d6 


E 

(S*2}(x)  =  £  v/^^g  ■  6)  dt  da  d6 


rv.5.3.  Overall  ensemble  averaged  effective  elastic  moduli 


We  assume  that  statistical  homogeneity  holds.  Therefore,  the  volume  average  is  equal  to  the 
ensemble  average.  Accordingly,  the  ensemble-volume  averaged  microcrack-induced  compliances 
take  the  form  ^ 

J A  Je  J-\  ^54^ 

/  /  r' vT^g .  (K)  rf{ 

^  ^  J A  J&  J-\ 

where  denotes  the  microcracks  density  parameter.  Note  that  (S*')  lepresents  the  first-order 

effects  due  to  the  presence  of  microcracks.  Microcrack  interaction  effects  are  represented  by  (S*^) 
which  is  of  the  second-order  in  u.  It  is  emphasized  that  our  second-order  compliance  formulation 
is  completely  different  from  the  dubious  second-order  stiffness  formulation  proposed  by  Hudson. 


Following  the  Voigt’s  notation  for  strains,  the  ensemble-volume  averaged  overall  constitutive 
equation  is 

(55) 

where  e  =  .  The  overall  ensemble- volume  averaged  compliance  consists  of  three 

components: 


(S)  =  S°  -t-  (S*»)  +  (S-2) 

in  which,  since  the  solid  is  assumed  to  be  isotropic  linear  "ilastic,  the  elastic  compliance  is 

1  —  1/  —V  0 

—  V  \  —  u  0 

0  0  2 


(56) 


(57) 


Remark  5.2. 1 :  The  effective  properties  normally  relate  the  average  strains  with  the  average  average 
stresses  or  mathematically,  (e)  =  (S)  •  {<t).  However,  in  Eq.  55,  the  relations  between  the  average 
strains  and  the  far  field  applied  stresses  arc  established.  This  defines  a  set  of  slightly  different 
effective  properties.  In  case  that  the  average  stresses  arc  very  close  to  the  far  filed  applied  stresses 
{tr°°  ~  (er)),  the  two  definitions  are  equivalent.  | 


IV.6.  Some  numerical  examples 


# 


# 


For  simplicity,  all  microcracks  are  assumed  to  be  of  equal  size  and  uniformly  distributed  in 
locations.  Legendre  polynomials  of  the  fourth  order  arc  adopted  as  base  functions.  The  number  N 
in  Eq.  (9)  is  therefore  set  to  4.  Twenty  unknowns  are  solved  from  twenty  equations  for  each  two- 
crack  interaction  problem.  Two  categories  of  problems  with  different  orientation  distributions  are 
considered.  In  the  first  category,  all  microcracks  are  assumed  to  be  alinged  (parallel)  yet  randomly 
located.  The  compliance  matrix  due  to  the  presence  of  microcracks  is 


_  \-iP-  0  0  0 

{S*l)  =  — — 2xu;  0  1  0 
^  0  0  1 


(58) 


which  is  in  consistent  with  the  Taylor’s  model  given  by  Eq.  (48).  On  the  other  hand,  only  two 
components  in  the  3  by  3  microcrack-interaction  induced  compliance  matrix  are  non-zero: 


(S-2)  = 


l-r/2 


-TTU; 


0 

ri 

0 


0 

0 

rz 


(59) 


The  COTstants  rj  and  r2  are  given  in  Table  8  for  two  cases  with  and  without  stress  projections. 
To  ensure  that  no  microcracks  intersect  each  other,  a  small  elliptical  probabilistic  region  surround¬ 
ing  a  microcrack  is  excluded  when  performing  the  numerical  Gauss  integration.  Major  axis  of  the 
elliptical  probabilistic  region  coincides  the  crack  line.  The  lengths  of  major  axes  for  all  computa¬ 
tions  are  twice  the  crack  size.  The  parameter  h  is  defined  as  the  ratio  of  minor  axis  to  the  crack 
size.  See  Figure  5  for  a  schematic  illustration.  A  smaller  k  value  implies  that  stronger  interacting 
microcracks  are  allowed  to  appear  in  the  neighborhood  of  a  given  microcrack. 


Table  8.  Results  for  aligned  microcracks  by  using  the  fourth  order  Legendre  polynomials. 


■ 

r2 

No  Projection 

Exact  Projection 

No  Projection 

Exact  Projection 

2.0 

1.253 

1.253 

0.489 

0.489 

1.5 

0.721 

0.721 

0.572 

0.572 

1.0 

0.248 

0.248 

0.618 

0.618 

0.5 

-0.006 

-0.006 

0.528 

0.528 

Randomly  oriented  and  located  microcracks  are  considered  in  the  second  category  of  examples. 


Much  more  computational  efforts  are  required  in  these  cases. 

_  1_,,2  0  0] 

(S'>)  = 0  1  0 
^  0  0  2 


(60) 


129 


Table  11.  Results  for  microcracks  with  random  orientation  by  using  the  center  point. 


B 

■8l 

52 

53 

No 

Projection 

Exact 

Projection 

No 

Projection 

Exact 

Projection 

No 

Projection 

Exact 

Projection 

2.0 

3.685 

3.497 

0.405 

0.364 

6.559 

6.265 

1.75 

3.443 

3.928 

0.236 

0.430 

6.411 

6.995 

1.5 

3.279 

4.487 

0.509 

6.432 

7.956 

There  are  three  non-zero  independent  components  in  the  microcrack-interaction  induced  compli¬ 


ance: 


(S*2) 


1-t/^ 

irE 


U} 


S\ 

•S2 

0 


S2 

-81 

0 


O' 

0 

^3. 


(61) 


Since  the  orientations  of  the  microcracks  are  random  and  the  original  material  is  isotropic,  the 
overall  compliance  should  be  isotropic.  In  this  case,  the  three  components  of  the  overall  compliance 
matrix  Eqs.  (57),  (60),  and  (61)  are  all  isotropic.  Therefor,  the  property  of  isotropy  is  preserved 
in  the  ensemble-volume  averaging  process.  Numerical  results  of  sj,  sj.  and  sj  for  cases  with  and 
without  the  operation  of  stress  projection  are  shown  in  Table  9. 


Table  9.  Results  for  microcracks  with  random  orientation  by  using  the  fourth  order  Legendre. 


k 

5l 

52 

53 

No 

Projection 

Exact 

Projection 

No 

Projection 

Exact 

Projection 

No 

Projection 

Exact 

Projection 

2.0 

1.036 

1.036 

0.014 

0.014 

2.046 

2.046 

1.75 

-0.364 

-0.364 

-0.257 

-0.257 

1.243 

1.243 

1.5 

-1.098 

-1.098 

-0.573 

-0.573 

0.408 

0.408 

As  a  comparison.  Tables  10  and  11  display  the  results  for  the  two  categories  by  using  only  the 
crack-center  points  approximations  (see  Ju  and  Chen  (1991a)]. 


Table  10.  Results  for  aligned  microcracks  by  using  the  center  point  interaction. 


■ 

ri 

>”2 

No  Projection 

Exact  Projection 

No  Projection 

Exact  Projection 

2.0 

1.173 

1.116 

0.481 

0.447 

1.0 

0.364 

0.022 

0.686 

0.542 

0.5 

0.640 

-0.398 

0.688 

0.372 
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From  Tables  8-11,  we  observe  that  there  is  no  difference  on  the  overall  compliances  with 
or  without  the  stress  projection  procedure  when  the  fourth  order  Legendre  family  is  employed  to 
approximate  the  interactions.  It  suggests  that  the  fourth-order  Legendre  approximation  itself  is 
accurate  enough  so  that  a  further  stress  projection  dose  not  affect  the  results.  On  the  other  hand, 
center  point  interaction  approximation  appears  to  be  less  accurate.  Different  results  arc  obtained  for 
the  two  categories  with  or  without  a  stress  projection.  It  is  observed  that  the  column  ri  in  Table  8 
exhibits  a  decrease  as  h  becomes  smaller.  Since  microcracks  are  aligned  in  this  case,  a  decrease 
in  the  value  allows  more  interaction  effects  from  the  random  second  crack  which  is  likely  to  be 
approxinjately  stacked  on  top  of  the  given  first  crack.  As  mentioned  in  Section  2,  stacked  cracks 
tend  to  interact  against  each  other.  A  shielding  effect  which  decreases  the  value  of  rj  should  be 
expected.  Therefore,  the  trend  in  the  ri  column  in  Table  8  is  rational. 
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rv.7.  Conciusions 


Accurate  microcrack  interaction  method  is  incorporated  into  the  framework  of  statistical  mi¬ 
cromechanical  model  to  estimate  the  effective  elastic  moduli  of  brittle  solids.  A  linear  system  of 
equations  is  obtained  through  the  numerical  function  approximation  method  by  using  the  family 
of  orthogonal  base  functions.  The  proposed  method  is  capable  of  solving  the  problems  of  very 
strongly  interacting  microcracks.  Higher  accuracy  is  achieved  by  increasing  the  order  used  in  the 
approximation  and  solving  the  resulting  larger  system  of  equations.  The  system  of  equ-dons  for 
many-microcrack  interaction  problems  can  be  constructed  in  a  similar  manner. 

% 

Complex  stress  potentials  are  shown  to  be  a  useful  tool  for  the  crack-interaction  and  compli¬ 
ance  problems.  Based  on  the  single  crack  stress  potentials  corresponding  to  concentrated  loads, 
relevant  quantities  such  as  stresses,  strains,  displacements  and  complainces  are  derived  for  any 
arbitrary  loadings.  Simple  and  useful  results  are  obtained  for  the  special  cases  of  polynomial  load¬ 
ings.  Extensive  numerical  computations  are  required  for  the  overall  ensemble-volume  averaged 
elastic  moduli.  Numerical  examples  show  that  the  fourth  order  Legendre  approximation  is  accu¬ 
rate  enough  so  that  further  stress  projection  is  not  necessary.  However,  the  stress  projection  should 
be  applied  to  improve  the  accuracy  if  center-point  interaction  approximation  is  used. 
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rv.9.  Figure  captions 


Figure  1. 
Figure  2. 
Figure  3. 
Figure  4. 
Figure  5. 
Figure  6. 


Coordinate  systems  for  two  microcracks. 

Three  cases  of  two-crack  interaction  configiuations. 

A  single  crack  subjected  to  a  pair  of  concentrated  nonnal  and  shear  loads. 
Symmetry  and  anti-symmetry  for  nonnal  and  shear  loadings. 

Active  integration  domain. 

Four  cases  of  microcrack  configurations. 
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Coordinate  systems  for  two  microcracks. 
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Case  1.  Two  collinear  cracks 
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Case  2.  Two  stacked  cracks 


Case  3.  Two  offset  aligned  cracks 


Figure  2.  Three  cases  of  two^rack  interaction  conSguradons. 
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figure  3.  A  single  crack  subjected  to  a  pair  of  concentrated  normal  and  shear  loads. 
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Symmetry 


Anti-symmetry 
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figure  4.  Symmetry  and  anti-symmetry  for  normal  and  shear  loadings. 


Case  4 


Figure  6.  Four  cases  of  microcrack  configurations. 
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